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ABSTRACT 


A numerical model is developed which calculates the atmospheric 
corrections to infrared radiometric measurements due to absorption 
and emission by water vapor, carbon dioxide, and ozone. The corrections 
due to aerosols are not accounted for. The transmissions functions for 
water vapor and carbon dioxide are those given by Davis and Viezee (1964). 
The transmission function for ozone is that given by Moskalenko (1969). 

The model requires as input the vertical distribution of temperature and 
water vapor as determined by a standard radiosonde. The vertical distri- 
bution of carbon dioxide is assumed to be constant. The vertical distribu- 
tion of ozone is an average of observed values. The model also requires as 
input the spectral response function of the radiometer and the nadir angle 
at which the measurements were made. A listing of the FORTRAN program is 
given with details for its use and examples of input and output litsings. 
Calculations for four model atmospheres are presented to illustrate the 
nature of the atmospheric correction for a variety of atmospheric condi- 
tions and for two different radiometers. 


v 
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INTRODUCTION 

The infrared radiometers used on aircraft and satellites to determine 
the surface temperature of the earth sense radiation in the most transpa- 
rent parts of the absorption spectrum of a cloudless, dustless atmosphere. 
Presently, the most used part of the spectrum is the 800 - 1200 cm“^ wave- 
number (8-14iim) atmospheric "window". This window is somewhat dirty due 
mainly to absorption by water vapor, carbon dioxide, ozone and aerosols 
{haze, smog, pollution, etc.). Most clouds are so highly absorbing 
(McDonald, 1960) as to make infeasible the detection of radiation from the 
earth's surface through the clouds. When viewing clouds the radiometer 
essentially measures cloud top temperatures. For a comprehensive treat- 
ment of the absorption characteristics of the atmosphere, the reader is 
referred to Goody (1964). 

In order to ascertain surface temperature from the radiance measured 
by the radiometer, the absorption and emission by water vapor, carbon 
dioxide, ozone and aerosols must be accounted for. This report describes 
the second version (Boudreau, 1972) of a computerized model which is being 
developed for calculating the atmospheric corrections to radiometric measure- 
ments made in the 800 - 1200 cm - ^ window. The atmospheric correction is 
defined as the difference between the actual temperature and the temperature 
as determined remotely by radiometer. This model is written ’n FORTRAN IV 
and calculates the atmospheric correction due to water vapor, carbon dioxide, 
and ozone only. The transmission functions for water vapor and carbon 
dioxide are taken to be those given by Davis and Viezee (1964). The trans- 
mission function for ozone is that given by Moskalenko (1968, 1971). 
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The vertical distribution of temperature, water vapor, carbon dioxide 
and ozone must be specified for the model from the surface to the altitude 
for which the correction is to be calculated. The variation in the verti- 
cal distribution of carbon dioxide is small and assumed to be constant at 
330 ppm. An average vertical distribution of ozone for 20°N given by 
Hering and Borden (1964) and shown in Fig. 1 is used since ozonesonde 
observations are not readily available. The model accepts as input the 
vertical distribution of temperature and water vapor as determined by a 
standard radiosonde. The model assumes that relative humidity is 10% when 
the humidity is low enough to cause "motorboating" to be reported in the 
radiosonde. Above the highest altitude for which measurements have been 
made by radiosonde and to a pressure altitude of 1 mb, the model uses the 
temperature structure from average soundings for 30°N which are given by 
Valley (1965). For the months of May through October and November through 
April, the 30°N soundings for the months of July and January, respectively, 
are used. The distributions for 30°N are used in the model since the model 
will be applied to remotely sensed data taken around that latitude, but it 
m^y be modified easily to incorporate other distributions. 

The water vapor profiles used for the upper atmosphere are those 
proposed by Clark (1973) based on the work of Goldman, et al_., (1973), 
Mastenbrook (1971), McKinnon and Morewood (1970), Murcray et_ ah , (1969), 
and Scholtz et al_. , (1970). Two water vapor profiles are used: a sunnier 
profile for the months of May through October and a winter profile for the 
months of November through April. These profiles are shown in Fig. 2. 

Since measurements of the vertical distribution of aerosol concentra- 
tion and composition are not usually available and since the transmission 
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Fig. 2. 


Proposed average water vapor profiles for the 
upper atmosphere for summer and winter (from 
Clark, 1973). 
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functions for aerosols are not well known, the model at this stage of 
development does not include the atmospheric correction due to aerosols. 
The present model was developed in part to study the effects of aerosols 
on remotely sensed infrared data. The difference between the actual 
atmospheric correction determined from a remote sensing experiment and 
the calculated atmospheric correction given by the model described here 
for the same experimental condition would be due mainly to aerosols. 

Hence, the gross aerosol correction can be studied empirically using 
this model . 

The atmospheric correction is specific to each radiometer, being a 
function of the radiometer's spectral response function, which describes 
the variation in the radiometer's response across the spectral band of 
800 - 1200 cnH . Therefore, the model also requires as input the spectral 
response function of the radiometer for which the atmospheric correction 
is being calculated. 
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DEVELOPMENT 

J 


Radi ometri call y determined temperature . If no atmosphere were 
present, an infrared radiometer at altitude z would detect a radiance, 

D 0 , expressed by 

v 2 

D 0 = J Y(v)S{v,e,T s )dv , (1) 

v i 

where v is wavenumber, y(v) is the radiometer's response function, v x and 
v 2 are the limits within which y(v)>o, S(v,T s ) is the radiant intensity 
leaving the earth's surface, T s is the surface temperature of the earth, 
and e is the nadir angle at which the radiometer is veiwing the earth's 
surface. Eq. (1) and the equations for radiance which follow implicity 
express the radiance integrated over the solid angle viewed by the radio- 
meter. 

Surface temperature may be obtained from measurements of D 0 as 
follows. Eq. (1) may be written as 

Do ~ y$(vq,6,T s )av, ( 2 ) 


where v 0 is defined such that 


S ( Vq 1T5 ) " / S(v »0»T S )dv/{v2-vi ) 9 


"1 


V2 


Vl 


~ f y ( v )s( v > 9 » t s ) dv / J S(v,0,T s )dv, 


VI 


VI 


( 3 ) 


( 4 ) 


and 
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Av = v - v, . (5) 

2 1 

Assuming that the earth is a black-body radiator, we take 

S(v,e,T $ ) = B(v,T s ) = av 0 3 / [exp (bv 0 /T s ) -1 ] » (6) 

where a = 8.9349 x 10“ 13 cal crr^sec^ster* 1 and b = 1.4385 cm deg. The 
use of (6) in (2) allows us to solve for T s » 

T s = bv 0 /ln(( Y Ava 3 0 /D 0 )+l ) . (7) 

The surface temperature determined by this method is referred to as 
an equivalent blackbody temperature or a brightness temperature since it 
is not the actual surface temperature unless the surface radiates as a 
blackbody. 

Sensing through a horizontally homogeneous atmosphere at a nadir 
angle e, a radiometer at height, w r , detects a radiance, D(w r ) , given by 
v2 l"w r 

D( w r) ~ f y(v) / B(v,w)^(v,(w r -w)f(e))dw 

vi o 


+ B(v,T s )t(v ,w r f(e) ) dv. 


( 8 ) 


where x is the transmissivity of the atmosphere, B(v,w) is the Planckian 
function given by (6), f ( e) a function which accounts for the increased 


atmospheric path length when viewing at nadir angle e (if a flat earth is 
assumed, f(e) = sece), and height, z, is expressed in terms of pnecipita- 
ble water, w. The relation between height and precipi table water is given 



o 
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where p and p' are the density of water vapor and liquid water, respectively. 
The Planckian function, B(v,w), can be expressed as a function of w because 
atmospheric temperature, T=T(w). A derivation of (8) is given by Boudreau 
(1968). 

Spherical earth geometry . Due to its spherical shape, the earth appears 
as a disk below the aircraft or satellite. The edge of the disk is located 
at a nadir angle, e 1 , which with the help of Fig. 3 is seen to be given by 

6* = arc sin (R/R+H) (10) 

where R = 6370 km, the radius of the earth, and H is the altitude of the 
sensor. Hence, we only need concern ourselves with data taken at e < 9'. 

For 0 < 9', the path length, ds, at z is given by Wark, et al_., (1963) as 

ds = [(R+z)m 2 dz]/[(R^z) 2 m J 2 -m 0 2 (R+z 0 ) 2 ] l/2 , (11) 

where n z and n Q are the indices of refraction of air at z and the surface, 
respectively, and z 0 is the minimum height of the refracted ray. Now, n z 
and n 0 are obtained from the formula 

n = 1 + CP/T (12) 

where P and T are air pressure (mb) and temperature (K) and 

C = 77.526 X 10 -6 K/mb. (13) 

% is obtained from Eqs. (2) and (4) of Wark et al_. , (1963), i.e., 


z Q = ((R+z)sine)/n 0 -R. 


(14) 
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Transmissivity functions . The transmissivity, t w , of water vapor 
for wave number intervals of 25 cm -1 in the 800 to 1200 cm -1 window is 
given by Davis and Viezee (1964) as 

b(v) 

t w (Av,w) = exp {-[L(v)wP] }. (15) 

Occasionally it is required to calculate atmospheric corrections for an 
instrument that has a response that extends slightly beyond the 800 to 
1200 cm -1 window. Outside the 800 - 1200 cm -1 window, Davis and Viezee 
(1964) give the following transmissivity functions. For water vapor 


transmissivity at v < 800 cm -1 , 

t w (Av,w) * exp{-AL( v)w(l+3. 1 7AL (v)w/P)”° • 5 ) , (16) 

in which 

A = [0.76+(0.58+0.48P 2 ) 0,5 ] p _0 • 1 (T/T 0 ) b(v ^ . (17) 

For water vapor transmissivity at v > 1200 cm' 1 , 

t w (Av,w) - exp{-BL(v)w(l+4.9BL(v)w/P)' 0 * 5 }, (18) 

in which 

B = [1 . 18+( 1 . 38+0 .48P 2 ) 0 » 5 ]p-° . 15 . (19) 

For carbon dioxide transmissivity at 550 <_ v <_ 800 cm" 1 , 

t c (Av,w c ) = exp[-0.4P°* e {£ck(v)w c +l]°‘ 5 -l }], (20) 

where 

c = (0.87+P°- e )P a (T/T o ) c < v ) (21) 
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and 


a = 1 .2 - 0.15 log P. 


( 22 ) 


In (15) through (22), P = p/1013.2, where p is pressure in millibars, 

T 0 = 273. 16K, w c is the optical path in atmos-cm of carbon dioxide, and 
values of L(v), k(v), b(v), and c(v) for increments of Av=25 cm -1 are listed 
in DECK 3, BLOCK DATA PROGRAM of the FORTRAN program given in the Appendix. 

The transmissivity, t 0 (6v,w 0 ) of ozone for irregular wave number 
intervals of < 10 cm -1 in the 602-847, 965-1175, 2000-2180 cm -1 bands is 
given by Moskalenko (1969, 1971) 

T„(«v.wo) = exp[-M'')t( 0 m(v) P n ( v )]. (23) 


where w Q is the optical depth of ozone and values of e(v), m(v) and n(v) 
for 6v <_ 10cm -1 are listed in DECK 3, BLOCK DATA PROGRAM of the program 
given in the Appendix. 

The ozone transmissivity is specified for 6v <_ 10cm" 1 and is matched 
to the 25cm -1 intervals of the water vapor transmission as follows. The 
ozone transmissivity, t 0 (av,w 0 ), for the 25* 1 cm intervals is taken to be 
Planckian - weighted transmissivity of the <$v intervals. 


t o( Av ’ w q) 



B(v,T)T 0 (<5v,w 0 )dv 


V 


(I 



B(v ,T)d 


v , 


I 

V 


( 24 ) 


where Av = v" -v 1 , the 25cm -1 interval of water vapor. The relation- 
ship between Av and 6v|< - v^+i - is shown in Fig. 4. Eq. (24) is 
calculated as (ommitting the w 0 symbolism for symplicity) 
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— 8v, 



K 8 v n“H 


» 


♦ 


V 


1 



V 


n+i 


v 


1 1 



FIG. 4. The relationship between the irregularly 
spaced ozone intervals, 8v k , and the larger 
(2 5cm" 1 ) water vapor intervals ,Av . 
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t 0 (Av) - ^CB-|(T)t(6v-| )(v 2 -v* ) + B2(T)t(6v2)(v3~V2) 

+ ••• + B m (T)T{6v m )(v"-v n] )] , (25) 

in which 

M = B-|(T)(v 2 -v') + B2(T)(v3-v 2) + + Bn (T)(v"-vn ) , (26) 

Bi(T) = B((v'+v 2 )/2,T), (27) 

B k (T) = B{ (v k +v k+ i )/2 ,T) , (28) 

B n (T) = B((v n +v")/2,T). (29) 


Since the highest ozone concentration are found in the stratosphere, an 
average T = 240°K is used in calculating (25) - (29). 

Water vapor, carbon dioxide, and ozone radiate mutually in some bands 
in the spectrum. For these bands, the transmissivity is taken to be the 
product of the individual transmissivities, i.e., 

T = t w (Av,w)t c (Av,W c )t 0 (Av,W 0 ) . (30) 

outside the mutual bands, t is taken to be t w . 

Numerical model . We solve a numerical analog to (8) since is can 
not in general be solved in closed form. In finite difference notation, 
(8) may be approximated as 

n I" m-1 

D ( w m) B(vj,Tj ) {t(vj ,(w m -Wi+i )f(e))- T (vj ,(w m -w-j ) 

j=l L i=l 


sec)} + B(vj »Ts ) x ( v j >w m f ( 0 ) )| Av j , 


(31) 
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where the Avj are constrained by the transmissivity functions to be 25cm’’ 1 , 
vj is the value of v at the midpoint of the 25cm -1 interval, and 

T. = (T i+1 + T-j )/2 ♦ (32) 


In the computation of the T(vj,W| < ) from the Davis and Viezee functions and 
the Moskalenko function, in place of P and T, we use an effective pressure, 
Pe. and temperature, T e , given by 

Wm . Wm m-1 



f dwi w^-iry] ™ t (X wl -X t )/2, 
J m k 

wu *=k 


(33) 


where X can be pressure or temperature, w is the optical depth of the 
atmospheric constituent for which t(av,w) is being calculated [e.g. w=w Q 
for t 0 (av,w 0 )J, and 


m i ~ W £ + l " W £* 


(34) 


The concentration of carbon dioxide is assumed to be constant; therefore, 
the effective pressure of carbon dioxide can be calculated from (33) as 
fol lows: 



(35) 


through the use of (9). The use of the ratio, q c , for carbon dioxide 


q c = Pc/Pa* (36) 

where p c and p a are the densities of carbon dioxide and air, respectively, 
and the hydrostatic equation. 
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dP 

lz = -Pa9» (37) 

in which g is the acceleration of gravity (98 cm sec' 2 ), in (35) obtains 



Now, q c = constant, therefore (38) may be integrated to obtain 

P e = <Pk + Pm>/ 2 - (39) 

The radiance D(w m ) detected by the radiometer is converted to an 
equivalent blackbody temperature, T r , through use of the calibration 
relati on 

v 2 n 

D c = y y(v)B(v,T r )dv = ^^Y(vj)B(vj,T r )Av j. (40) 

vi j=l 

Eq . (40) is used to calculate a table of versus T r over the range of 
T r that is expected to be encountered when using the radiometer. Then we 
interpolate in this table to find T r that corresponds to the value of 
D(w m ) obtained from (31). The atmospheric correction, aT s , is given by 

aT s - T s - T r , (47 ) 

where T s is assumed to be the air temperature at the surface as reported 
by the radiosonde observation. 

Optical depths . Precipi table water, w, is computed from temperature, 
dew point temperature, and pressure as follows. The hydrostatic equation, 
(37), is used to write (9) as 
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w 



dp/p 'g 


Po 


L 



q 44/p' g, 


4 = 1 


(42) 


where q is the average specific humidity in the layer Ap. 

In (42), q is obtained from the relation (Haltiner and Martin, 1957) 

q = 0.622e/(p-0.378e) > (43) 


in which vapor pressure, e, is determined from dew point temperature 
using the Goff-Gratch relation specified in List (1971). 

The optical depth of carbon dioxide in atm-cm is computed from 
z 

w c = p-. f P C < 44 > 

0 


where is the density of carbon dioxide at STP. The use of (36) and 
(37) in (44) gives 


P 



which, since q c = O.SxlO -3 , may be integrated to obtain 


(45) 


w c = 0.260 (p 0 - p) 


(46) 


where p is in millibars. 

The optical depth of ozone in atm-cm is computed from 
2 P 


W 0 = -f- L f »£.dp, 

P n * po 9 J Pa 


P l 

o 


(47) 
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by (37). In (47), p 0 is the density of ozone and p' 0 is the density of 
ozone at STP. Letting 



(47) is calculated from 


L 



*=1 


where q 0 is the average value for the layer Ap. 


(48) 


(49) 
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SAMPLE CALCULATIONS 

For purposes of illustrating the atmospheric correction model and the 
nature of atmospheric corrections, calculations are presented for four 
typical distributions of water vapor and temperature given by Valley (1965) 
and for two radiometers: the NOAA-2 SR and the RS-18. The response func- 

tions for the NOAA-2 SR and the RS-18 are shown in Fig. 9. The vertical 
distributions to a pressure altitude of 100 mb of temperature and dew 
point for the four model atmospheres are shown in Figs. 10 - 13. The 
atmospheric corrections which correspond to the four model atmosphere are 
shown for the N0AA-2 SR in Fig. 14 and for the RS-18 in Figure 15. 

As can be seen in Fig. 9, the response of the NOAA-2 SR avoids the 
ozone and carbon dioxide bands and is therefore only sensitive to absorp- 
tion and emission by water vapor. Fig. 14 shows that for the four model 
atmospheres the atmospheric correction is positive and increases with 

height. In the upper levels (> 300 mb) the correction increases more 
slowly than at lower levels because specific humidity becomes very low. 

As the amount of water vapor approaches zero with increasing height, the 
atmospheric correction approaches a constant value with height. 

The spectral response of the RS-18 radiometer is much broader (see 
Fig. 9) than that of the NOAA-2 SR; it's response extends into the carbon 
dioxide band and contains the ozone band. Fig. 15 shows that the atmos- 
pheric correction for the RS-18 does not approach zero with increasing 
height above 300 mb but due to emission and absorption by carbon dioxide 
and ozone continues to increase to approximately 10 mb and then decreases 
with increasing altitude. The decrease in correction above 10 mb is due 
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Wavelength (urn) 


Fig. 9 . Response Functions for the RS-18 and NOAA-2 SR. 






20 



Fig. 10. Vertical Distribution of Temperature (Solid Line) and dew 
point (dashed line) for Model #1 - Tropical Storm. 
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Fig. 11. Vertical Distribution of Temperature (Solid Line) and Dew 
Point (Dashed Line) for Model #2 - Sub-tropical Summer. 
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Fig. 12. Vertical Distribution of Temperature (Solid Line) and Dew 
Point (Dashed Line) for Model #3 - Mid-latitude Summer. 
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Fig. 13. Vertical Distribution of Temperature (Solid Line) and Dew 
Point (Dashed Line) for Model #4 - Mid-Latitude Winter. 
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vs. Pressure Altitude for the RS-18 Radiometer 

















to the increase of temperature with height in the stratosphere (~100 - 
1 mb) such that at about 10 mb the temperatures are greater than the 
effective radiating temperature of the atmosphere below. The maximum 
concentration of ozone occurs around 10 mb and the increased emission due 
to ozone at these warmer temperatures causes the radiance to increase 
with height above 10 mb which results in decreasing the atmospheric cor- 
rection. The RS-18 is an airborne scanning radiometer designed for ob- 
taining low altitude infrared imagery. It is obviously ill-suited for 
satellite application due to the effects of ozone and carbon dioxide. 

The temperature lapse rates do not vary significantly in the model 
atmospheres. As such, the magnitude of the atmospheric correction is 
primarily a function of the amount of water vapor. The atmospheric 
correction for the NOAA-2 SR varies (see Fig. 14) from 9.8C for the very 
humid atmospheric model 1 (Tropical Storm) to 1.7C for the relatively dry 
atmospheric model 4 (Mid-latitude Winter). Similarly, the atmospheric 
correction for the RS-18 being flown at a pressure altitude of 700 mb 
(~ 10,000 ft.) varies from 4.4C for the humid Tropical Storm to 1.8C for 
the drier Mid-latitude Winter model. 

The above calculations were made for zero nadir angle. The effect of 
viewing at a nadir angle greater than zero would be to displace the curves 
in Figs. 14 and 15 to the right, i.e. toward larger atmospheric correction. 
This larger correction is due to the longer atmospheric path encountered 
when the surface is viewed at a nadir angle greater than zero. The longer 
atmospheric path increases the magnitude of the atmospheric correction. 

Atmospheric corrections for the Eglin, AFB, Florida, 1115 GMT, 21 
November 1971 radiosonde (shown in Fig. 16) illustrate the effects of a 
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Fig. 16. Eglin AFB, Florida, 1115 GMT 21 November 1971 Radiosonde. 
Temperature (Solid Line) Dew Point (Dashed Line). 
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temperature inversion at the surface on atmospheric correction. The tempe- 
rature increases approximately 5 degrees C from the surface to 1500 feet. 
The emission from this layer of air which is warmer than the surface causes 
the atmospheric correction as shown in Fig. 1 7 to be negative to an alti- 
tude of 14,000 ft. Above 14,000 ft. the correction becomes positive. 



Altitude ( 10 3 feet 



Atmospheric Correction vs Altitude for Eglin AFB, Florida, 1115 
GMT, 21 November 1971, Radiosonde. PRT-5 at Zero Nadir Angle. 
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SUMMARY AND RECOMMENDATIONS 

A numerical model has been developed which calculates the atmospheric 
corrections due to water vapor, carbon dioxide and ozone. Future versions 
of this model should include the corrections due to aerosols so that the 
model will have more general applicability. When more general versions of 
this model are developed, they should be subjected to the test of compari- 
son with observation. Of course the present model can be used to calculate 
atmospheric corrections when aerosol concentration is low. Furthermore, the 
difference between the actual correction determined from a remote sensing 
experiment and the calculated correction given by the radiation model for 
the same experimental condition would be due mainly to aerosols. Hence, 
the gross aerosol correction can be studied empirically using the model. 
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FORTRAN PROGRAM OF MODEL 
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General , In order to facilitate the comparison of the model as de- 
scribed by formulae in the text and the FORTRAN program which follows, 
a listing of FORTRAN symbols and their corresponding text symbols is given 

in Table 1 . 


Table 1. FORTRAN symbols and their corresponding text symbols. 



FORTRAN TEXT 

A 


DEFINITION 

A temporary variable used in transmission 
computations 


ABSOR 

*"v 

An absorption coefficient for transmission 
through water vapor in the wave number 
range 25 to 2150 cm"* 

ABSORC 

k v 

An absorption coefficient for transmission 
through carbon dioxide in the wave number 



range 550 to 800 cm"’ 

ABSORO 


An absorption coefficient for transmission 
through ozone in the wave number range 
575 - 2150 cm -1 


ANGLE 

0 

Observation angle from nadir to target 

ATCENT 

W 

c 

A table of optical depths of carbon dioxide 
in atmosphere centimeters 

ATMOS 

B('-j-Tl) 

The blackbody function for wave number 
and temperature T. 

B 


A temporary variable used in transmission 
computations 

BB0DY 

B( Uj ,T.) 

The blackbody function for wave number vj 
and temperature T. 

BETA 

3 

A temporary variable used in transmission 
computations 

C 


A temporary variable used in transmission 
computations 
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Table 1. 

continued 


FORTRAN 

TEXT 

DEFINITION 

CPRESS 


Effective pressure of carbon dioxide 

D 


A temporary variable used in transmission 
computations 

DEGREE 


A table of temperatures whose midpoint is 
the surface target temperature 

DELTRA 


The change in transmissivity between layers 

DETECT 

D < w m> 

o 

The detector intensity, cal cm sec 

DEWAVG 


Average dewpoint of an atmospheric layer 

DEWPT 


Dew point temperature 

DLTEMP 

aT 

s 

Temperature correction from surface to 
sensor 

DRTEMP 

T r 

Simulated temperature observed by sensor 

EXPRESS 

p e 

Effective pressure of water vapor 

ETEMPS 

T e 

Effective temperature of water vapor 

EXPOMO 

m(v) 

An absorption exponent for transmission 
through ozone in the wave number range 
600 to 2150 crrf‘ 

EXPON 

b 

V 

An absorption exponent for transmission 
through water vapor in the wave number 
range 800 to 1200 cm"* 

EXPONC 

c 

V 

An absorption exponent for transmission 
through carbon dioxide in the wave num- 
ber range 550 to 800 cm"* 

EXPONO 

n(v) 

An absorption exponent for transmission 
through ozone in the wave number range 
600 to 2150 cm"* 

FUNCTI 

Y 

The percent response of a sensor to 
energy in a given wave number interval 

HEADER 


A table for instrument names 

HEAD1 

HEAD2 

HEAD3 


Temporary variables for these instrument 
names 
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Table 1. continued 



FORTRAN TEXT DEFINITION 


HEATER D 

c 


A blackbody intensity table simulating 
instrument response to a blackbody 


HEIGHT 

z 

A table of radiosonde observation altitudes 

HIGH 


Average height of an atmospheric layer 

HUM 

q 

A table of specific humidities for each 
atmospheric layer 

HUMID 


Average water vapor mixing ratio of an 
atmospheric layer 

I 

i 

Atmospheric layer index 

INTEN 


Sum used to determine integral detector 
intensity 

INTER 

v j 

A table of wave number mid-points 

INTERV 

AVj 

A wave number interval 

J 

i 

The wave number interval index 

LEVELS 


The number of levels at which radiosonde 
data is given or is interpolated 

LIM 

m-1 

The maximum number of layers to the level 
of observation 

LIMIT 

m 

The number of levels of data available 

ONE 


1.0 

OTCENT 

w o 

Optical depth of ozone in atmosphere 
centimeters 

OZONE 

% 

A table of ratios of the density of ozone 
to its density at STP for each atmospheric 
layer 

PART 

q 

A specific humidity for determining precip- 
i table water 


Table 1. continued 


FORTRAN 

TEXT 

DEFINITION 

PARTI 
PART2 
PART 3 
PART4 
PART5 
PART6 


Temporary variables used to compute vapor 
pressure 

PRATIO 

P 

The ratio of the average pressure within 
an atmospheric layer to a standard of 
1013.2 mb 

PRESS 

P 

A table of pressures 

RESPON 

Y(v) 

A table of response functions defining 
an instrument 

SFCBB 


The contribution from a surface target 
image 

SLOPEA 

SLOPEB 


Temporary variables used to perform 
linear interpolation on blackbody inten- 
sity functions 

STEPSZ 


The desired step size between atmospheric 
levels 

TEMPER 

T. 

1 

A table of atmospheric temperatures 

TEMPOR 

T i 

The average temperature of an atmospheric 
layer 

THETA 


Angle in radians 

TRANGE 


The width in degrees of the temperature 
calibration table 

TRANS 

T 

The transmissivity of the atmosphere 

TRANSC 

t c (Av, w c ) 

The transmissivity at a level due to 
carbon dioxide absorption 

TRANS M 

tCvjW^) 

The transmissivity of the atmosphere from 
the surface to level k 

TRANSO 

t 0 (5v,w 0 ) 

The transmissivity at a level due to ozone 
absorption 

TRANS ER 

t w (Av,w) 

The transmissivity at a level due to water 


vapor absorption 
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Table 1 . continued 


FORTRAN 

TEXT 

DEFINITION 

TRATIO 

T ' T o 

The ratio of a level's average temperature 
to a standard temperature of 273. 16K 


373.16/T 

The ratio of the steam point to a saturation 
temperature when used to compute vapor pres- 
sure 

VAPER 

e 

The vapor pressure used to compute precip- 
i table water 

WATER 

w 

A table of precipitable centimeters of 
water vapor 

WATERC 


The amount of carbon dioxide in atmos -cen- 
timeters for a given layer at a given angle 

WATERO 


The amount of ozone in atmosphere centimeters 
for a given layer at a given angle 

WATERS 


The amount of water vapor in precipitable 
centimeters for a given layer at a given 
angle 

WAVENO 

V 

The wave number mid-point 

WHOLE 


A temporary sum used to determine integral 
precipitable water and upwelling radiation 
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Program Description . The description of the FORTRAN program which 
follows is taken from Hal bach (1973). The program is written for use on 
a UN I VAC 1108 EXECUTIVE system and has various types of computation schemes 
and output displays which are appropriate to the users needs. The program 
control flow charts (Figs. 7 and 8) and lead card descriptions (Tables 2-7 
and Figs. 5 and 6) give a description of these options. The program requires 
25,000 locations in core. The user should be familiar with the radiative 
transfer equation and the logic of the program which is outlined in the 
program control flow charts (Figs. 7 and 8) in order to execute the options 
available in the program. The formats for input of data to the program 
are given in FORMAT statements in the program listing. 
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Table 2. An Example of Lead Card Arrangement 


Field Col uirois Format Name Identification 

Card No . 1 , Name CONTROL 

1 1-10 110 LIMIT An integer number telling how many 

non-control cards follow or defining 

a value used for control. 

2 11-17 A6 LABEL Six alphanumeric characters defining 

the type of data or action which 
follows : 

'TEMPER' Temperature range 

'ANGLE' Viewing angle from nadir 

'STEP' Pressure step size 

'SPHERE* Apply sphere function 

'FLAT' Do not apply sphere function 

'MAXCOR' S-C4020 minimum temperature 

' MINCOR' S-C4020 minimum temperature 

'MESSAG' One line message 

' RESPON ' Response function data 

'TRACE' Transmission function data 

'WATER' Radiosonde data 

'EXECUT' Use previous radiosonde 

'LIMB' Make limb display 

' DEFAUL ' Load mission radiosonde data 

3 23-28 A6 HEAD1* Six alphanumeric characters used as 

supplemental control data. 

'METERS' Change meters to feet 


*HEAD1 , HEAD2, and HEAD3 are sometimes used as display labels. 


Table 2. (Continued) 
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4 

29-34 

A6 

HEAD2 

5 

35-40 

A6 

HEAD3 


Card 

No. 2, Name 

WATER 


1 

2-8 

F7.1 

PRESS (I) 

2 

11-16 

F6.1 

DEWPT 

3 

21-27 

F7.1 

TEMPER 

(I) 

4 

31-40 

FI 0.1 

HEIGHT 

(I) 

5 

41-50 

E10.5 

OZONE ( i) 

Card 

No. 3, Name 

RESP0N 


1 

1-7 

F7.1 

INTERV 

2 

8-14 

F7.4 

FUNCT1 

Card 

No. 4, Name 

TEMPER 


1 

1-7 

F7.1 

TRANGE 

Card 

No. 5, Name 

MESSAG 


1 

2-61 

10A6 

MESSAG 

(I) 

Card 

No. 6, Name 

TRACE 


1 

1-10 

E10.5 

WEIGHT 

(I) 

2 

11-20 

E10.5 

ABSORO 


(I) 


Six alphanumeric characters used 
as supplemental control data. 

'KELVIN' Do not change temperatures 

Six alphanumeric characters used as 
supplemental control data. 

'OZONE' Load ozone data from cards 

Atmospheric pressure in millibars 
Dew point at pressure level 
Temperature at pressure level 

Height at pressure level 

Ozone at pressure level 

Wave-number midpoint interval 
Percent response of instrument system 
in this interval 

Temperature range chosen for calibra- 
tion table 

Lead message line for displays 

Index and weighting function for 
trace gas transmission coefficients. 
Coefficient for transmission function. 


Table 2. (Continued) 
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3 

21-30 

E10.5 

EXP0M0 

(I) 

4 

31-40 

E10.5 

EXP0N0 

(I) 

Card No. 

7, Name 

DEFAUL 


1 

2-8 

F7.1 

PRESUR 

2 

11-20 

E10.5 

OZONES 

3 

23-28 

F6.1 

TEMP0R 

4 

31-40 

E10.5 

HUMID 


First exponent for transmission function 
Second exponent for transmission function 

Atmospheric pressure in millibars 
Ozone at pressure level 
Temperature at pressure level 
Specific humidity at pressure level 
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DECK SETUP FOH PROGRAM 


Radiation Model 


PAGE NO. -1— OF- 



j 8ASG,T PROGAM, T,YYYYY 

r " ' " ' 

7 

8|{DG PROGRAMMER ERL BLDG. NO. PHONE NO. 


8RUN I, (P)(/OPT)J (RUNID , ACC OUNT . PRO J- ID . RUN-TIME . PAGES ) 


(Front of dick) 


Fig. 5. Example of a Lead Card Setup 


1 ‘ * 5. f< 030 
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Table 3. Card Input and Printer Output Using Card A (Fig. 5). 


Input 


1 

e 

IM'i.OzrctL 

AltkACt 

.23119*02 . C-'.WW-IT? 

.95000 +00. 

. 1 00-00 + VO 

3 

.24 223*02 .699’./*- 92 

.95000.00 

.100-00 +0-0 

4 

. 2 4 30 ! * 02 .1 300 - 0 1 

. 950.0.0*00 

. 10.000 +1-0. 

5 

.24316*02 .320’ .0-01 

.90.000m 0.0 

. 1 5000 +0.0 

e 

.24156*02 , 50000-01 

.90000+00 

. J 50.00+ O', i 

7 

.25162* 02 . 500' <0 -0 1 

.90000*00 

, 1 50.00+00 

6 

.25322* 02 . COOOO - 0 1 

. 66009 +t-0 

. 1 70-0.0. +0.0. 

9 

.25316* 02 . 90000 -0 1 

.8500/. : 0-0 

. 16000 *00i 

10 

♦ 2 5 1 97 * 02 , 1 ?' *99+ 00 

.6 500-0 *00 

.16010+00 

11 

.£6142 *02 . 1 2‘ '■90 * 00 

.65000+00 

, 1 600'.'. ♦ 00 

12 

. £63 63*02 . 1 5000* 00 

.64000*00. 

.20000+00 

13 

. 263 57 * 02 . 1 0000 ♦ 00 

. 620.00* 0*0- 

.220 VO ♦CO 

14 

. 26 1 30* 02 , 1 0000 ♦00 

.62000*00 

, 22000+00. 

15 

.27224*02 . 2217-0*00 

.62000+00 

.22V.OO+OOi 

16 

.27362* 02 . 3017.70 * OO 

,6217.70+0/0 

.25000+00 

17 

.27277* 02 . 36000* 170 

. 617.7.70 ♦ OO 

.25000+00 

16 

.27 1 37* 02 .230*00 * OO 

.63000+00 

.221/7. it tV.i 

19 

.261 03*02 . 230170*00 

.63000+00 

. 22000+00 

20 

. 26303 * 02 . 200* 70 * *-70 

.64000 + 170 

.20000 + lO x 

£1 

. £8 1 00+ 02 . 2*v_ 7,70 ♦ *70 

.62000 + 170 

.29/VO9++0O 

£2 

. £6434 **0£ . 4 400*0* OO 

,60000+00 

.2CO0O+0O 

£3 

. 29204 ♦ 0? . 4417, >0* '.70 

.60000*00 

.2Ci'oo++;o 

£4 

. 292*02 ♦ 02 . 30000* 170 

.60000 +00 

. 2217,7,. * 00. 

25 

. 294 3 7* 0*2 . 20*000 ♦ 00 

.60000*170 

.22000*00. 

2b 

.29157*02 . 2017.70 *170 

.60000*00 

.2200,0+00 

£7 

. 302 £ 5*02 . 1 417.70* 170 

.64000+00. 

.24000+00 

£6 

.30*479* 02 .11 17.70*170 

. 6 CO 00+00. 

.24000+00 

£9 

.3023 5*02 .11 (7.70 *170 

. t'CO-'OO+OO 

.24100:00 

30 

.31245*02 .9*7.7.70-01 

.6a77.i*lV0 

.22000*00 

31 

.31 5L'0 *02 . 70-04/ /-Or! 

.69000 + 00 

.200.00+001 

32 

. 3 1 235*02 . 700470-01 

. 69000.+ 00. 

.20000+00 

33 

. 3224 5 * 02 . 50-900-91 

.90000+00 

, 16000*00 

34 

, 32 5 59 * 02 . 300470-01 

.95000+001 

.15000+00 

35 

. 3 £ 1 9 5 * 02 . 304700 -0-1 

.9 59+99 *09 

.1 5000*00. 

36 

.33327+02 .1 017.7. 1-0 1 

.95000*00 

, 1 1 004 1 + OO 

37 

.33557*02 *70il700-02 

.96009*09 

. 1 OOOO + OO 

36 

. 33 1 1 7* 0£ . 701700-02 

. 96090+ OO 

. 1 oooo+ o/.-« 

39 

. 38600*02 . 35 300 -01 

.10000+0-1 

.00000 

40 

.30197*02 .35 3170- 0 1 

.10000*+ 0.1 

4CV.7.OO 

41 

.38194*02 .890-0/0-01 

- 79500+00. 

.0.0000 

4£ 

. 39206* 02 . 8917.70-0 1 

.79500+ 00. 

.00000 

43 

. 39403*02 . 3400/0*170 

.91200+00 

,99000-01 

44 

.39391*02 .61 7170* 99 

.65400+00 

. 1 5000+00. 

45 

. 40409 ♦ 02 • 2030/0* 01 

.91900*00 

.23500+00 

46 

.4039 7+ 02 . 4030/0 *0 1 

.91500+00 

.31000+00 

47 

,40104 + 02 * 5 3900. 1 0 1 

. 925'- ’.U 

. '.V* 


. . r *+ /s/*» . 4' 4 

nir iv: i r/’. 


■»*.• 

4i 4 i r.o * c 7r.i'< t 

o+ 6 *+■.<; 

,£660* .♦ OO. 

50 

.41391 + 02 . 33600 + 0 1 

.65100+00. 

.26100+00 

51 

.42410*02 , 70300 + 0* 1 

. 90400+00. 

.34999*90 

52 

.42397*02 . 1 5 50-' '» + 0- J 

. 7690-0 + 00 

.20500 + 00 

53 

.42194 + 02 .164 O', i + 00 

.76200+00 

.86500-01 

54 

. 4320 7*02 . 1 6 4 00- * 00 

.70200+ OO 

.66500-01 

55 

.43403*02 .97000-01 

, 72900 ♦ on 

.35000-01 

56 

. 4 30 1 9 ♦ 02 . 9 7 00-0-0 1 

.72900+0-0 

. 3f0.no- -0«i 


Output 




INTERVAL 

A 

e 

C 

HEIGHT 

375 . 0 - ec - o.o 

.9060 

.9500 

, 1009 

.1199 

€ 00 . 0 - 625.0 

. 00/60 

. 950*0 

, 1990 

,2230 

600.0— 625,0 

.0139 

.9500 

.1009 

.3919 

600 . 0 - 625.0 

.0329 

.9000 

. 1 50/9 

.3180 

600 . 0 - 625 . U 

.0500 

,9000 

,1590 

.1580 

625 . 0 - 650.0 

.05 09 

. 99*90 

.1500 

. 162*0 

625 , 0 - 650,0 

.0609 

. 680*9 

. 171 X 1 

.3229 

625 . 0 - 650.0 

.0990 

.8500 

.1800 

.3189 

€ 25 . 0 - 650.0 

.1200 

.8590 

.1600 

.1979 

650 . 0 - 675.0 

.1200 

.8590 

.1800 

*1420 

650 . 0 - 675.0 

. 1 500 

.8400 

.2900 

.3630 

650 . 0 - 675.0 

.1600 

.8200 

.2299 

.3570 

650 . 0 - 675.0 

.1600 

.8290 

.2290 

.1380 

675 . 0 - 700.0 

. 220+0 

.8299 

.2299 

.2240 

€ 75 . 0 - 7 T - C .0 

.3000 

,8200 

.2500 

.3620 

675 . 0 - 700,0 

.3600 

.eooo 

.2500 

.2770 

675 . 0 - 700.0 

.2300 

.8300 

.2200 

.1379 

700 . 0 - 725.0 

.2300 

. 830/0 

.2299 

,1639 

700 . 0 - 725.0 

.2000 

.8400 

.2090 

.3930 

700 . 0 - 725.0 

.2800 

.8200 

.2000 

.1000 

700 . 0 - 725.0 

.4400 

.8090 

.2690 

.4340 

725 . 0 - 750.0 

,4400 

, 80*00 

.2690 

.2940 

725 . 0 - 750.0 

.3000 

,8000 

.2299 

.202 a 

725 . 0 - 750.0 

.2000 

.8009 

.2200 

.4370 

725 . 0 - 750.0 

.2000 

.8000 

.2200 

.1570 

730 . 0 - 775.0 

.1400 

.8400 

.2400 

.2859 

750 . 0 - 775.0 

.uoa 

.8609 

.2409 

.4799 

730 . 0 - 775.0 

♦ 1 too 

. 86**0 

.2400 

.2350 

775 . 0 - 600,0 

.0900 

.8600 

.2200 

.24 50 

775 . 0 - eoo.o 

.0700 

.8900 

.2000 

.5209 

775 . 0 - 600.0 

.0700 

.6900 

.2000 

.2350 

800.0- 625.0 

,0500 

. 9000 

.1800 

.2450 

600 . 0 - 625.0 

.0300 

.9500 

. 1 590 

.5590 

600 . 0 - 625.0 

.0300 

.9509 

,1590 

. 1950 + 

825 . 0 - 630.0 

,91 09 

. 950.9 

.1100 

.3279 

625 . 0 - 650.0 

.9079 

.9699 

.1900 

.5579 

625 . 0 - 650.0 

.0979 

,9609 

.1900 

.1179 

350 . 0 - 975.0 

.93 53 

1 . 900*0 

. 0*999 

, 61*90 

950 . 0 - 975.0 

.0353 

1 . 0*099 

.0009 

.1979 

950 , 0 - 975.0 

.9899 

, 7959 

.0900 

.1949 

975 . 0 - 1000.0 

♦9899 

,7959 

.9900 

.2969 

975 . 0 - 1000.9 

.3409 

.9120 

,9990 

.4939 

975 . 0 - 1000.0 

.8170 

.8549 

.1500 

.3919 

ir + 00 , 0 - 1025.0 

2 . 930*9 

.9199 

.2350 

.4999 

100 - 0 . 0 - 1023.0 

4 .0300 

.9159 

.3100 

.3979 

i ■> V , « , 0 .- 10 * 25 . 0 

3 . 390/0 


. 3 3 90 

* 194 ’-l 

1025 n-t-'.MV n 

* *eiiv. 

Coro 


. r'ro 

to ? 5 . 0 - 1050.0 

4 . 67*,/ 1 

.91 AO 

.2669 


1025 . 9 - 1050.0 

3.3690 

.8519 

.2619 

.3919 

1950 . 9 - 1075.9 

7. 9390 

.9949 

.3400 

.4199 

1050 . 9 - 1075.0 

1.5 590 

,7899 

.2959 

.3979 

1959 . 0 - 1975.9 

. 1 840 

. 7629 

.9865 

.1949 

1075 . 9 - 1100.9 

.1849 

.7629 

.0865 

.2079 

1075 . 0 - 1100.9 

.0976 

.7290 

.9359 

.4039 

1075 . 9 - 1109.9 

.0976 

.7299 

.0350 

.0199 



oxer 

ossc* 

onv* 

Q9PC* 

0 w 
0602' 
0161* 
0S6C* 
on** 

09 PC* 
GW 
0602* 
oxer 

GS6C* 

GUV* 

0161* 

0C6I* 

GZ6C* 

ontv 
006 1* 
oco** 
OZCiZ* 
fl£6f 
W65* 
ooxv* 
l*.U6* 


001 * 3 * 

oocp* 

0090 * 

□• 0412 - 0*5212 

0022 * 

0009 * 

0055 * 

O' OS 12 - 0" 5212 

0012 * 

00*2 * 

0095 ’ 

0 ’ 05 12 - 0*4212 

0002 * 

OC 9 Z* 

009 * * 

0 * 5212 - 0*0012 

0002 * 

0009 * 

0095 * 

0 * 5212 * 0*0012 

0002 * 

0 O 9 Z* 

0005 * 

0 * 4212 - 0*0012 

0002 * 

DG 9 Z* 

0005 * 

U'OOIZ— 0 ' 4 Z 02 

0002 * 

0092 * 

OOP** 

0 * 0012 ^ 0 * 5 Z 02 

0002 * 

0002 * 

0092 * 

0 * 0012 - 0 * SZ 02 

0002 ’ 

0092 * 

GGM* 

D" 5 Z 02 -O*O 502 

0002 * 

DOOP* 

0990 * 

0 * 5 / 02 - 0*0402 

0091 * 

0009 * 

0550 * 

o* 5 zo 2 -a"osn 2 

OOP I* 

0009 * 

0550 * 

0 * 0502 - 0*5202 

0002 * 

0009 * 

02 MT 

0 *DSOZ - 0*5202 

0002 " 

0029 “ 

O 2 C 0 * 

0 * 0402 - 0*5202 

0002 * 

0029 * 

02 CO* 

0 * 5202 - 0*0002 

OOfiO* 

GOGS* 

2 * 20 * 

O’SZI 1 - 0*04 It 

0000 * 

0009 * 

zcco* 

0 ‘ 5 ZIt-O*O 5 tt 

sw 

091 2 * 

2 /VU* 

0 " 4 Zt 1 - 0*04 II 

06 W 

0529 * 

Z 990 * 

0 * 0411 - 0*5211 

tawi * 

0 * 69 * 

non* 

0 * 0411 - 0*5211 

cwr 

0606 * 

GIH* 

0 * 0511 - 0*4211 

G*sr 

0606 * 

GIVI* 

0 * 5211 - 0*0011 

IVS 2 I* 

0299 * 

0901 " 

0 * 6211 - 0*0011 

two’ 

0062 * 

ocvz* 

0 * 5211 - 0*0011 

0140 * 

GZ 92 * 

9 Z 9 G* 

0 * 0011 - 0 * 4/01 


ix+rz’zziS* 

LXl +00022 ’ 

00 + 000 x 2 * 
vo + 00002 * 
00+00002 ‘ 
00* 00002' 
lu * '. 4/. A. 1 2 
00*00002* 
00+00002* 
00*00002* 
00+ V0002 ' 
00*00091* 
iVi+OGGP i * 
■X+0Ci002* 

10 * ♦ o r .v.o z * 
00 * 0000 * 2 * 
00000* 
1*0000 * 
ig-oosvc* 
XG-0006>‘ 
00+00201 * 
00+00V5 l * 

00+ 00 VS X' 

oo+oocz 1 * 

10*00015* 


L'Ol +-VZ1GC 9 * 
OOjHVXOi*’ 
•*OiH‘OOi>i>" 

00 < UO09Z ‘ 

oo+vzvoop* 
00+000 3Z* 
GO : GO 09/ * 
00 » ov*o 9 z * 

VZ>OOOPZ“ 

* w+ooopz * 

00+00009* 
',‘0* +00009 ’ 

oo+oooop* 

00+00009* 
00+00029* 
■7j +00029’ 

vzv+oooos* 

OO+OVOO** 

oo+ GO 3 I t * 
co* + iv.wzfl * 

00+00*69* 

001 +00606 4 
GO +00606 4 
GO +00299* 
•VO +0006* 
0G+0GZ9Z* 


lOrO’*COOJ9* 

tOi+OOOCC* 

CZl+VOOSS" 

iw+oglw* 

oo+o.opc* 

*.Oj+;."OOOS * 
VO* COCO'S* 
VO+VOO’3>* 
VO+ 00092* 
OO+OGOn* 
10—1*0099 * 
10-0005 5 ■ 
10-00054* 
10-0002 V* 
xo-ooozt* 
UrW.tZt* 
1 0-002*2* 
ig-oozcc* 
tG-ooz/v* 
10-00*99* 

oo+ooo i i 4 

vo#*roiv i' 

r*'.**c6t> i ' 

00+0G901* 

Wi+V-GCVZ* 

io-og ur 


20+16159* 

20+56CS1T 

2G+-VIV59* 

20+995*9* 

ZG+WW 

20+603*9* 

ZO-IGlCP* 

zo+secc** 

ZG+VIVC9* 

20+99529* 

20+60230“ 
20+I6U9* 
20+SGCI9* 
zo+?i>X9* 
2G+IG109* 
20+56 19 V* 
20+Z659V* 
20+01 V9V" 
20+0654 V 
2G+5GW 

2G ♦/t'i2iO 

ZOtCGlW* 
2G + * 600 
20+OTVW* 
2 G + US 5 V* 


20 

Id 

op 

6Z 

PZ 

ZZ 

<*z 

cz 

rz 

tz 

2Z 

tz 

HZ 

69 

*9 

Z9 

99 

59 

V9 

59 

29 

19 

V .9 

64 

PS 

zs 


(panu. 1 q.u 03 ) •£ aiqei 


£i? 
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Table 4. Card Input and Printer Output Using Card B (Fig. 5) 


Input 


Output 



urt.nu 

-5 



1 


?5ilF.5fCtJS f*T~3 

WAVENUMBER Ml CPC 

z 

637.5 

.02 

637.5 

,0209 

3 

662 , 5 

.10 

662.5 

.1009 

4 

C87.5 

.25 

667.3 

.2509 

3 

712.5 

.45 

712.3 

.4509 

e 

737.5 

.63 

737.3 

. 6390 

7 

762.5 

.74 

762,3 

.7499 

e 

767 . 5 

.7* 

767.3 

.7699 

9 

612.5 

.*3 

612.5 

.6390 

10 

637.5 

.91 

637.3 

.9100 

11 

*62.5 

.92 

662.3 

.9299 

12 

667.5 

.92 

667.3 

.9299 

13 

912.5 

.93 

912.5 

.9309 

14 

937.5 

.96 

937.3 

.9699 

15 

962.5 

.99 

962.3 

.9909 

16 

9*7.5 

1 .00 

967.5 

1 .9900 

17 

lfil 2 . 5 

.99 

1012.3 

.9900 

16 

1037.5 

,S6 

1037.5 

.96 09 

19 

1062.5 

.91 

1062.3 

.9199 

20 

10*7.5 

.6G 

1037.5 

.6600 

21 

1112.5 

.62 

1112.5 

.6200 

22 

1137.5 

• 79 

1137.5 

.7990 

23 

1162.5 

.76 

1162.5 

.7600 

24 

11*7.5 

.71 

1187.5 

.7109 

25 

1212.5 

.56 

1212.3 

,5600 

26 

1237.5 

.22 

1237.3 

.2290 
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Table 5. Card Input and Printer Output Using Card C (Fig. 5) 


Input 


vrt.Hxtn 


1 


MCSSAG 



? 

m:ccl number i t iK*:r 

1 C Ai STORM 

AHCk R HTML *11954; 

3 

19 WATER 

MEirRSKtLVlN 

4 

1 Ol 3. Z 

300.2 

303.7 

• 

5 

930.0 

297.6 

301 .0 

577. 

6 

9tHl.LT 

295.3 

296.3 

1057. 

7 

eso.r. 

293.1 

256.1 

1 560 . 

B 

6' 

290.3 

293.4 

2066 . 

9 

730.0 

267.3 

290.7 

2645. 

ll'i 

TOO .« 

264.7 

267.6 

3233. 

11 

650.0 

261 .2 

264.4 

3656. 

ia 

600 ,T» 

277.6 

261.1 

4525. 

13 

550 .0 

274.0 

277.5 

5239 . 

14 

500.0 

271.3 

272.9 

6010. 

15 

<50,r. 

265.2 

266.0 

6647. 

16 

400.0 

259.6 

262.3 

7763. 

17 

350.0 

253.2 

255.9 

6776. 

16 

30*0.0 

244.9 

246.2 

9916. 

19 

250.0 

234.2 

236.5 

11215. 

20 

200.0 

221 .2 

226.5 

12734, 

21 

150.0 

373.2 

211.0 

14577. 

22 

100.0 

373.2 

169.2 

16950. 


Output 

H3CEL NVXJER l f lRCf>KAL STCftH AFTER RIEHL«1954> 


PRESSURE 

CEW POINT 

TEMPERATURE HEIGHT 

FRECIP 

ATHjS- 

CENT 

EFFECTIVE 

CONSTANT 

EFFECTIVE 

1013.2 

Z7.0 

30.5 

FEET 

WATER 

COtf 

GZCNE 

PRESSURE 

CONCENTRATION 

TEMPERATURE 

950.0 

24.4 

27.8 

1693.0 

1.3699 

16.4 

.001 

961 .6 

961.6 

302,3 

900.0 

22.1 

25.3 

3467.6 

2.3621 

29.4 

.003 

957,8 

956.6 

301,3 

850.0 

19.9 

22,9 

5116.0 

3,2758 

42.4 

.004 

934,7 

931.6 

300.2 

800.0 

17.1 

20.2 

6650.3 

4.1 046 

55.4 

.005 

912.6 

906.6 

299.1 

750.0 

14.1 

17.5 

6677.7 

4.6363 

66.4 

,006 

891.7 

681.6 

298 .0 

700.0 

H.5 

14.6 

10606.6 

5.4916 

81.4 

.007 

871 .9 

856.6 

297.0 

630 .O 

8.0 

11.2 

12657.3 

6.0645 

94.4 

.006 

653,3 

831.6 

295.9 

600.0 

4.6 

7.9 

14845.6 

6.5532 

107.4 

.009 

636.3 

606.6 

294,9 

550.0 

.6 

4.3 

17168.1 

6.9657 

120.4 

.010 

620.8 

781.6 

294.0 

500 .O 

-1.9 

-.3 

19717.6 

7.3229 

133.4 

,011 

606.4 

756.6 

293.1 

4 30.0 

-*.o 

-5.2 

22463.6 

7.6U/4 

146.4 

.013 

✓94 .O 

731.6 

292.3 

4UO.G 

-13.4 

-10.9 

25466.9 

7.M08 

159.4 

.014 

V64.4 

✓06.6 

291 .6 

350,0 

-20.0 

-17.3 

26796.9 

7.9519 

172.4 

.015 

777.1 

681.6 

291.0 

300.0 

-26.3 

-25.0 

32532.4 

8.0374 

165.4 

.016 

772.3 

656.6 

290.6 

250.0 

-39.0 

-34.7 

36794.2 

8.0769 

198.4 

.017 

769. 7 

631.6 

290.3 

200.0 

-52.0 

-4G.7 

41777.7 

6.0936 

211.4 

.016 

768.7 

606.6 

290.2 

130.0 

100.0 

-62.2 

47624.2 

8.0941 

224.4 

.019 

768.7 

561.6 

290.2 

100.0 

100.0 

-84.0 

55609.6 

8.0941 

237.4 

-020 

768.7 

556.6 

290.2 



Table 6. Card Input and Printer Output Using Card P (Fig. 5) 


Input Output 

Otf MJLTEC STRATOSFHtRC PRCFILCS 


1 

2 

iffi . dlf AULT /A 

34CCFAUL 

112.4 .36910-06 

203.1 

.13000-02 

PRESSURE OZCWt TEMP 

^DUNGING 

K) GM/RGM K 

112.4 ,36910-06 203 . 1 

SPECIFIC 
HUM! ClTt 
CH/KCM 
.131730-02 

5 

95.2 

. 50000-06 

205.2 

.14100-02 

95.2 

.30000-06 

205.2 

.14100-02 

4 

60. a 

.90000-06 

207 , 4 

.15300-02 

80.8 

.900043-06 

207.4 

.15300-02 

5 

68.7 

.16000-05 

209.6 

.16500-02 

68.7 

.160470-05 

209 . 6 

.16500-02 

6 

56.5 

.25000-05 

211.7 

.17900-02 

56.5 

.25000-05 

211.7 

.17900-02 

7 

49.6 

. 40000-03 

213.9 

. 19400-02 

49.8 

.417300-05 

213.9 

.19400-02 

e 

42.6 

.50000-05 

215.9 

.21LV.VV02 

42.6 

.30000-05 

215.9 

.211730-02 

9 

36,4 

. 55000-05 

217.9 

.22600-02 

36.4 

.5 5' 730-05 

217.9 

.226173-02 

10 

31 .2 

.70000-05 

219.9 

.24700-02 

31 .2 

.700473-05 

219.9 

.Z4nV3-02 

11 

26.6 

.60000-05 

211.6 

.26600-02 

26.8 

.8000/3-05 

211.6 

.26800-02 

12 

23.0 

. 1 01V.V3-04 

223.8 

.29000-02 

23.0 

.10000-04 

223.8 

.291730-02 

13 

19.6 

, 1 00' VO -04 

225,8 

.314'30-02 

19.6 

.100473-04 

225,8 

.31400-02 

14 

17.0 

.10000-04 

227.8 

.341 IV. i-02 

17.0 

. 1007-73-04 

227.8 

.341173-02 

13 

14.7 

.10500-04 

229.8 

.36900-02 

14.7 

.105073-04 

229.8 

.36900-02 

16 

12.7 

.105173-04 

231 .8 

.417//30-02 

12.7 

.10507.3-04 

231.6 

. 4017.V3—02 

17 

11.0 

.105*70-04 

233.7 

.400173-02 

11.0 

.10500-04 

233.7 

.40/300-02 

1* 

9.5 

.1 10* VO- 04 

235.7 

.4173173-02 

9.5 

.11000-04 

235.7 

.4017.73-02 

19 

6.2 

.10900-04 

238.0 

.400470-LT2 

8;2 

.10900-04 

238.0 

.40000-02 

20 

7.2 

. 1 0MV0-O4 

240.4 

.404770-02 

7.2 

.10500-04 

240,4 

.40000-02 

21 

6.2 

. 1 0000-04 

242.7 

.40000-02 

6.2 

.10C4XVO4 

242.7 

.40000-02 

22 

5.4 

.950470-05 

245.1 

.41VV.V0-02 

3.4 

.95000-05 

245.1 

.40/300-02 

23 

4 .7 

.91VVV0-O5 

247.5 

.41V.V/O-02 

4.7 

. 900730—05 

247.5 

• 40/3173—02 

24 

4.1 

.6LV.V_V.V05 

249.6 

.40047/1-02 

4.1 

.80000-05 

249,8 

.40000-02 

25 

3.6 

.70000-05 

252.3 

.417////0-O2 

3.6 

.70000-05 

252.3 

.400/30-02 

26 

3.2 

.69000-05 

254.6 

,4047.7/. -02 

3.2 

.69030-05 

254.6 

.40000-02 

27 

2.6 

.600.70-05 

256.9 

, 4 '/47/41-02 

2.8 

.60300-05 

256.9 

.40000—02 

2^ 

2.4 

.5600/0-05 

259.3 

. 4 04V. WO 2 

2.4 

. 58000—0 5 

259.3 

.40000-02 

29 

2.1 

. 54LXV.V05 

261 .7 

.400173-02 

2.1 

.54000-05 

261.7 

.40000-02 

30 

1 .9 

.54000-05 

264.1 

.4l7/4V/i-02 

1.9 

.54000-05 

264.1 

.4000*3-02 

31 

1.6 

.4517.73-05 

2G6.4 

.4047/4*1-02 

1.6 

.45000-05 

266.4 

.40000-02 

32 

1.5 

,40000-05 

206.6 

,4047.70-0.2 

1.3 

.40000-05 

268.6 

.40000-02 

33 

1.3 

.30047/1-05 

271 .2 

.4047/41-02 

1.3 

.38000-05 

271.2 

.40000-02 

34 

1 .1 

.3GOCO-C.5 

272.1 

.4047.70-02 

1.1 

.36000-05 

272.1 

.40000-02 

35 

1 .o 

.35000-05 

272.1 

.40000-02 

1.0 

*35000—05 

272.1 

.40000-02 



Table 7. Printer Output Using Card E (Fig. 5) 


PRESS 

sr-tciFic 

TEMP 

OZCNC 


humidity 

SCONCING 

►© 

CM/KCM 

K 

CM/KCM 

1013.2 

,22232 *02 

303.70 

.46700-07 

963. Z 

.21*49*02 

301 . 56 

.46700-07 

913.2 

.19663*02 

299,14 

.46700-07 

663.2 

.1631 7+02 

296.73 

. 46700-07 

©13.2 

.10061*02 

294.10 

.46700-07 

763,2 

.14639*02 

291 .41 

, 46700-07 

713.2 

,13236*02 

266.55 

. 46700-07 

663.2 

.11646*02 

263.27 

. 46700-07 

613.2 

.10010*02 

261 ,97 

. 46700-07 

563.2 

.64664*01 

276.43 

.46700-07 

513.2 

.73121*01 

274.06 

,46700-07 

*63.2 

.39364*01 

269.26 

.46700-07 

*13.2 

.43936+01 

263.75 

.46700-07 

363.2 

.31067*01 

257.54 

.46700-07 

313.2 

.19737*01 

250.15 

.*6700-07 

263.2 

.10364+01 

240.93 

.46700-07 

213.2 

.45164+00 

229.51 

.46700-07 

163.2 

.96410-01 

214.67 

.46700-07 

113.2 

.20772-01 

194.56 

.46700-07 

63.2 

,11597-01 

210.73 

.21775-05 

13.2 

.24221-02 

231.30 

.10500-04 


NU-CtR 1 , TRCPICAL SKRM AFTER RIEHLH954) 


PRT-5 

ANGLE .0 DEGREES FROM NADIR 
SURFACE TEMPERATURE 303.7 


PRESSURE 

FREC1P 

HEIGHT SPECIFIC 

TEMP 

CORRECT ICN 


WATER 

humidity 

SCLMDING 


►e 

CM 

FT CM/KGM 

K 

K 

963.2 

1 .0943 

1491.5 ,21449*02 

301 .56 

.23 

913.2 

2.1076 

3049.2 .19665+02 

299,14 

.77 

663.2 

3.0424 

4679.9 .18317*02 

296.73 

1.52 

613.2 

3.6934 

6391.6 .16681*02 

294.10 

2.40 

763.2 

4.6515 

8193.6 .14859+02 

291.41 

3.36 

713.2 

5.3266 

10097.6 ,13236+02 

288.55 

4.36 

003 . 2 

5- 9210 

12116.5 .116*6*02 

285.2/ 

5.4s: 

61 3 .2 

A .d^17 

r nviin*ri» 

Ml 07 

c 

563.2 

€.6647 

16569.5 .64864*01 

278.4^ 

7.56 

313.2 

7.2376 

19049.3 .73121*01 

274 .06 

8.66 

*63.2 

7.5406 

21735.8 .59384*01 

269.26 

9.76 

413.2 

7.7649 

24670.0 .43936+01 

263.75 

10.80 

363.2 

7.9234 

27907.8 .31067*01 

257.54 

11.77 

313.2 

8.0241 

31526.3 .19737+01 

250.15 

12.67 

263.2 

6.0780 

35634.2 .10564*01 

240.93 

13.49 

213.2 

8.1011 

40399.4 .45184+00 

229.51 

14.27 

163.2 

6.1061 

46107.9 .98410-01 

214.87 

15.03 

113.2 

8.1072 

53306.6 .20772-01 

194.56 

15.66 

63.2 

6.1078 

64660.6 ,11597-01 

210.73 

17.35 

13.2 

6.1079 

97933.5 .24221-02 

231 .30 

20.45 
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Dashed line encloses optional functions not needed for program execution. 

Fig. 7. Program Control Flown Chart 1. 




5i 


Label 

Execute 


Label 

Water 


Load 

Radiosonde 

Data 


■"^Pressure 1 
Step Size 


Expand Or 
Contract 
Data To 
Step Size 


Display 

Recomputed 

Radiosonde 

Data 













Program Listing 



52 


0 

ham: 

VERSION 

TYPE 

DATE 

TIME 

sea 

5IZE-FRE 

TEXT 4CYCLE W3RO PSRHXE 

LOCATION 


VAFRES 

SLANT 

SYM3XIC 

17 JUL 73 

11 *01 '55 

i 


5 

3 

0 

1 

1792 


TEMTAS 

slant 

SYMiXlC 

17 JUL 73 

11*01*42 

2 


16 

S 

Q 

1 

1797 


CCEF 

SLANT 

$YK>XlC 

17 JUL 73 

11*01*58 

3 


28 

5 

0 

1 

1813 


TRANS 

SLANT 

SYMBOLIC 

17 JUL 73 

11*01*47 

4 


26 

5 

0 

1 

1841 


KA IN 

SLANT 

SY|©CL 10 

17 JUL 73 

11*02*05 

5 


128 

5 

0 

1 

1667 


KAIN 

SLA NT 1 

SYFCCLIC 

17 JUL 73 

11*03*10 

6 


126 

3 

0 

1 

1995 


RACMCC 

SLA NT 1 

SYHXLlC 

17 JUL 73 

11*02*59 

7 


25 

3 

0 

1 

2121 


RA£>t=C 

SLANT 

SYPBX1C 

17 JUL 73 

It '01*51 

8 


27 

3 

0 

1 

2146 


TUX507 

SLANT 

ABSOLUTE 

17 jul n 

1102*54 

9 


563 




2173 


TLCCXIT 

SLANT 1 

ABSOLUTE 

17 JUL 73 

11 *03*59 

10 


546 




2736 


VAFRES 


relocatable 

17 JUL 73 

11*01*57 

11 

1 

4 




3282 


TEMTAB 


RELOCATABLE 

17 JUL 73 

11*01*44 

12 

1 

6 




3287 


COEF 


relocatable 

17 JUL 73 

11*01 ‘59 

13 

1 

20 




3296 


TRANS 


relocatable 

17 JUL 73 

11*01 *49 

14 

2 

20 




3317 


K0CEL1 


ELT SYKG 

13 SEP 72 

15*36*21 

IS 


7 

3 

0 

1 

3339 


M3CELZ 


ELT SYK© 

13 SEP 72 

15*36*22 

16 


7 

5 

D 

1 

3346 


M3 CEL3 


ELT SYFB 

15 SEP 72 

15*36*23 

17 


7 

3 

0 

1 

3353 


MODEL 4 


ELT SYKB 

15 SEP 72 

15*36*24 

18 


7 

3 

0 

1 

3360 


MXJEL5 


ELT $Y© 

23 APR 73 

09*57*13 

19 


T 

3 

0 

1 

3367 


MOCEL6 


ELT SYM3 

15 SEP 72 

15*36*25 

20 


7 

5 

D 

1 

3374 


MODEL 7 


ELT SYHB 

15 SEP 72 

15 ‘ 36*25 

21 


7 

3 

0 

1 

3381 


MODELS 


ELT SYI© 

17 JUL 73 

11*04*11 

22 


11 

3 

0 

1 

3388 


MCCEL9 


ELT 3Y« 

17 JU- 73 

It *04 * 12 

23 


7 

3 

0 

1 

3399 


MCOCIO 


ELT SYKB 

17 JUL 73 

11*04*12 

24 


6 

9 

0 

1 

3406 


MOOCH 


ELT STI© 

17 JUL 73 

11*04*12 

25 


6 

3 

D 

1 

3412 


MODE 12 


ELT SYI© 

17 JUL 73 

11 ‘04*12 

26 


11 

5 

a 

1 

3420 


KCEl 3 


ELT St^C 

IT JUL 73 

11*04*13 

27 


6 

3 

0 

1 

3431 


OZONE 


ELT SYK) 

07 APR 73 

lid* $9 * 10 

28 


24 

3 

0 

1 

3439 


FRT-5 


ELT SY*3 

09 JUN 72 

15*56*12 

29 


4 

3 

D 

1 

3463 


Rs-ia 

8-1 3 

ELT SY*« 

15 SEP 72 

15*36*39 

30 


4 

3 

0 

1 

3467 


default 

A 

ELT SThB 

31 JUL 73 

09*35“ 20 

31 


10 

3 

0 

1 

3471 


default 

B 

ELT SY*€> 

31 JUL 73 

G9* 35 *21 

32 


10 

3 

0 

1 

3481 


racm:c 

SLANT 

SYFDX1C 

31 JUL 73 

09*35*42 

33 


27 

5 

0 

1 

3491 


RArn:c 


RELOCATABLE 

3i JUL 73 

09 1 3 5 * 4 4 

34 

Z 

13 




3518 


TRANS 

SLANT 1 

SYKXLIC 

31 JUL 73 

09*35*49 

35 


26 

5 

0 

1 

3533 


TRANS 


RELOCATABLE 

31 JUL 73 

09*35*51 

36 

2 

20 




3559 


TRANS 

SLANT 

SYrtl-XtC 

31 JUL 73 

09*35*53 

37 


26 

S 

0 

1 

3581 


TRANS 


RELOCATABLE 

31 JUL 73 

09*35*55 

38 

2 

20 




3607 


MAIN 

SLANT 

SYfOXlC 

31 JUL 73 

09*35*58 

39 


133 

3 

□ 

1 

3629 


MAIN 


RELOCATABLE 

31 JUL 73 

09*36*04 

40 

4 

140 




3762 


RADMAF 

SLANT 

map srve 

31 JUL 73 

09*3605 

41 


2 

3 

0 

1 

3906 


TLtWOT 

SLANT 

ABSOLUTE 

31 JUL 73 

09*36*24 

42 


601 




3908 


TRANS 

SLANT 1 

SYFOXIC 

31 JUL 73 

09*36*27 

43 


26 

5 

0 

1 

4509 


TRANS 


relocatable 

31 JUL 73 

09*36*30 

44 

2 

20 




4535 


RACMOC 

SLANT 1 

SYtCCLlC 

31 JUL 73 

09*36*32 

43 


25 

3 

0 

1 

4557 


RADMIO 


RELOCATABLE 

31 JLL 73 

09*36*34 

46 

2 

9 




4582 


MAIN 

SLANT 1 

SYK5X1C 

31 JUL 73 

09*3639 

47 


132 

5 

0 

1 

4593 


MAIN 


RELOCATABLE 

31 JUL 73 

093645 

48 

3 

97 




4725 


TLC4X17 

SLANT 1 

ABSOLUTE 

31 JUL 73 

09 1 37 *04 

49 


556 

' 



4825 


RACMOO 

SLANT 

SYhCJCLlC 

20 AUG- 73 

18*50*07 

50 


27 

5 

0 

1 

5381 

* 

RACMCC 


RELOCATABLE 

20 AUG 73 

18*50*10 

51 

2 

13 




5408 


TRANS 

SLANT 

SYVCVXIC 

20 AUG 73 

18*50*11 

52 


26 

5 

0 

1 

5423 

♦ 

TRANS 


relocatable 

20 AUG 73 

18*50*13 

53 

2 

20 




5449 


VAi-Ktb 

SLANI 

STMAA.il. 

<0 AULr O 

id _>0 13 

>4 



3 

c 

i 

m i i 


VAFRES 


relocatable 

20 AUG 73 

18*50*16 

55 

1 

4 




54 76 


TEMTAS 

SLANT 

SYTOCLlC 

20 AUG 73 

18*5017 

56 


16 

5 

0 

1 

5481 


TEMTAB 


relocatable 

20 AUG 73 

18“ 50*19 

57 

1 

8 




5497 


COEF 

SLANT 

5YMDXIC 

20 AUG 73 

18*50*19 

58 


28 

5 

0 

1 

3506 


CCEF 


relocatable 

20 AUG 73 

18*50*21 

59 

1 

20 




5534 


MAIN 

SLANT 

SYICCLIC 

20 AUG 73 

td'SO^ 

60 


135 

5 

0 

1 

5555 

* 

MAIN 


RELOCATABLE 

20 AUG 73 

18* 50 *31 

61 

4 

142 




5690 


nsvoT 

SLANT 

ABSOLUTE 

20 AUG 73 

18*50*47 

62 


601 




5836 


TRANS 

SLANT1 

SYKXLIC 

20 AUG 73 

18*50*48 

63 


26 

3 

0 

1 

6437 


TRANS 


RELOCATABLE 

20 AUG 73 

10'5O'5i 

64 

2 

20 




6463 


RACMCC 

SLANT 1 

Srb&3LlC 

20 AUG 73 

18*50*52 

65 


23 

5 

0 

1 

6485 


ftACMOO 


RELOCATABLE 

20 AUG 73 

18*50*54 

66 

2 

9 




6510 


MAIN 

SLANTS 

SYKjCLIC 

20 AUG 73 

18*50*56 

67 


132 

3 

0 

1 

6521 


MAIN 


relocatable 

20 AUG 73 

18*51 ‘03 

68 

3 

97 




6653 


TLOOD7 

SLANT 1 

absolute 

20 AUG 73 

18*51 *21 

69 


554 




6753 


TESTl 


ELT 5YM3 

20 AUG 73 

18*51*22 

70 


3 

3 

0 

1 

7307 


*EXT AVAILABLE LOCATION- 
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FCR HAlN/SLANTl »MAlN,KAtN/SLANTl 

FC« 94L-Ci*/2G-ld'30 «0,0l 


HA IN FROORAM 

STCfiAtC OSCC’ CCCE'tl ) W.2S47 CATA*r.l KWM CLAf* COWDNiZ] CWtV.tfO 
ccw:n blocks* 

0003 21 001-164 
OCK 14 22 002570 

0005 23 030005 


EXTERNAL REFERENCES -ELCCK, NANt) 


0006 

VAFRES 

t»J7 

taflot 

CXilO 

ICFRMV 

OOI1 

6RIC1V 

0012 

FLTfC 

0013 

flterm 

0014 

tactem 

0015 

rackxj 

0016 

LIKCV 

0017 

NXV 

0020 

NYV 

0021 

NlNlni 

0022 

UvlAJi 

0023 

woe* 

0024 

NUCU$ 

0025 

NlOlS 

0026 

ALCC 

0027 

NS TCP 3 

0030 

>CKR3$ 

0031 

ASIN 

0032 

SIN 

0033 

SORT 

0034 

t-£RR4$ 


STORAGE ASSlOffrCNT -CLOCK. TYPE, RELATIVE LOCATION , NAMT) 


OOOI 

002044 

10L 

OOOI 

002223 

10526 

onno 

fW*K»23 

12F 

oooo 

0003 1 5 

16F 

oroo 

000354 

21 F 

0000 

0004 26 

2GF 

oooi 

0003 13 

3006 

oooo 

000514 

35F 

0001 

000737 

4720 

oooi 

001134 

5546 

oooi 

00074 4 

700L 

OOOI 

001032 

74 OL 

0003 R 

000000 

AD SCR 


oooi 

000576 

10L 

OOOI 

002277 

1070G 

nnnn 

000575 

i?r 

OOOO 

000320 

17 F 

oooo 

000363 

22F 

OOOI 

000267 

2640 

OOOI 

000132 

300L 

0001 

002272 

40L 

oooi 

002104 

5L 

0001 

000731 

COOL 

000 ft 

000771 

71 OL 

DGOl 

001742 

7736 

0003 ft 

000620 

ABSORC 


OOOI 

002161 

10L 

OOOO 

000551 

1 IF 

oooo 

Oty'&ip 

1 %F 

oooo 

000334 

1SF 

oooo 

0004 I 4 

23F 

oooo 

000427 

Z7F 

oooo 

000441 

31F 

DGOl 

000700 400L 

0001 

000722 

500L 

0001 

001236 6076 

OOOI 

001024 

72 OL 

0001 

001174 

800L 

0003 ft 

00 J 320 ACSCKO 


0001 

000013 

oooo 

000220 

noon 

000250 

oooo 

000344 

OOOO 

000416 

OJOG 

000431 

OOGO 

000465 

OGOl 

000640 

OOOI 

001037 

0001 

001267 

OOOI 

001514 

0001 

001366 

OOOO R 

000131 


1«JL 0001 

nr oooi 
\4F ryvyi 
19F DGOl 
Z4F 0000 
ZSF 0000 
32 F 0000 
4306 000 ft 
5246 COO I 
6260 0001 
7210 0001 
S10L 0001 
A MAX 


002206 10410 
001403 1100L 

*yvi 3 i 3 i?f 
000620 200L 
01-0422 25F 
000433 29F 
000506 33F 
G00643 435G 
001123 5450 
001417 662G 
001027 730L 
001307 615L 


0000 ft 000130 AMIN 
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DGOS 

R 

000001 

ANCLE 

0004 

R 

000764 

ATCENT 

0000 

R 000200 

BEGIN 

0004 

ft 001750 

CPRESS 

0004 

i 

002260 

CEFAUL 

0004 

ft 

001440 

CEL1RA 

00*73 

R 

000126 

CEST 

0000 

R 000153 

CEWAYG 

0000 

R 000U6 

CEVJT-NT 

0000 

R 

1X3014 5 

CCVTT 

OCK-14 

ft 

OOl 604 

EPRESS 

0004 

ft 

0021 1 ; 

ETEH“S 

GOOO 

R 000143 

EVHCLE 

0003 

R 001010 

ExPCI O 

0003 

ft 

000144 

EXPCN 

1 >003 

R 

000*32 

ExPctJC 

0003 

ft 

0011 54 

EXPCNO 

001X3 

R 000201 

FINIS 

0000 

R 000155 

FUNCTI 

GOOO 

R 

0001 37 

HEAC1 

0CAX1 

R 

0001 40 

HCAC2 

0000 

R 

IV.V.'i 141 

HCAC3 

OOtX 

R 000157 

HCIGHS 

0004 

R 000454 

HEIGHT 

0000 

R 

000160 

HIGH 

0004 

R 

0110620 

HUH 

o&:o 

R 

1X70154 

HUMIC 

0000 

ft 000174 

HUM i DP 

00110 

I 1X30135 

T 

0000 

I 

GOO 161 

IHIGH 

oooo 

t 

0001 76 

M 

OLXX* 


000664 

INJP* 

0000 

0170702 

INJP3 

0000 

X. 7.1 71 1 

I njt$ 

GOOO 


0X676 

1NJFI 

oooo 


000672 

iNJFf 

GOOO 

I 

COX72 

IN STR 

017.13 

R 000310 

INTER 

0000 

R 17.10000 

INTER V 

0000 

I 

000156 

J 

0000 

I 

0001 62 

ft 

ixvx 

l 

000136 

LABEL 

0000 

R O' 7X7.1 1 

LEADER 

0005 

I 000004 

LEVELS 

GOOO 

ft 

000004 

LIFE 

D 00* 

I 

cv.vxvx 

LIMIT 

OOOTi 

I 

00021 3 

LCCINT 

0000 

I 0X0215 

LCCSIZ 

0000 

I 0000 SC 

KtSSAG 

0000 

I 

000070 

NAME 

0017 

1 

oooiy.vji 

NXV 

0020 

I 

'yyyyy.'i 

NYV 

0004 

R 001130 

OTCENT 

17.7.14 

R 002424 

ozcue 

0000 

R 

0001 72 

C2CNEP 

OLV.VT 

R 

000175 

OZCKES 

gooo 

R 

'.V.Xil 50 

FART 

01X7.1 

R C 7X204 

PARTI 

0000 

R 000205 

PART? 

0000 

R 

000206 

PAuT3 

OOtV.l 

R 

000207 

PART4 

oooo 

R 

1X70 1 33 

Pn;rc 

0004 

R 001X7.70 

PRESS 

C77X 

R 000171 

FRESSP 

0000 

R 

0001 52 

PREStn 

GtXXl 

fl 

000214 

FRCWAN 

0053 

R 

0004 54 

ftEST-CN 

0005 

L 1X7X03 

SHELL 

0000 

R 000211 

SLCPEA 

0000 

R 

000212 

SLOPED 

0004 

R 

001663 

SPHERE 

GOOO 

R 

1X7.1202 

SPHERE 

GOOO 

R 17X132 

STEPSZ 

D0O4 

R 000310 

TEMPER 

0000 

ft 

0001 51 

tempo; 

0000 

R 

000173 

TEFCRE 

oooo 

R 

0001 63 

TEKT1 

avx 

ft 000164 

TEm"-2 

0001 

R 001720 

TERP 

0000 

R 

000210 

TERP 

OCH.VO 

R 

000203 

TEST 

0005 

R 

0X7X12 

THETA 

0000 

ft 000134 

trance 

0004 

R 001274 

TRANSM 

OOX 

R 

000144 

TWiXE 

oux 

R 

000147 

VAPEft 

0006 

R 

O.X7.X7J 

VAFRES 

0004 

R 000144 

WATER 

0000 

ft CA7.ll 7 7 

LEIGH 

O.VJ3 

R 

1X10644 

WEIGHT 

GOOD 

R 

0001 42 

vh:le 

o:/yj 

R 

000167 

X 

0X0 

R 000165 

XF 

0000 

R 000 1 70 

Y 

axo r 

000166 

TF 


OG1G1 

1* 

00103 

2 * 

anas 

3* 

00106 

4* 

otnoT 

5* 

oona 

6 + 

00111 

74 

05112 

6 * 

00113 

94 

00114 

104 

00114 

114 

00114 

124 

00114 

134 

00115 

144 

00115 

154 

00115 

16* 

00115 

174 

00115 

104 

00116 

194 

DD117 

204 

0012 a 

214 

00121 

224 

00122 

234 

0012,5 

244 

IAJ1 «iA 


GO* £5 

264 

00127 

274 

00131 

284 

00133 

294 

00135 

304 

00141 

314 

00141 

324 

00143 

334 

00145 

344 


COT' I LEKoDATA= SKXT ) 

EXTERNAL VAFfiES 
REAL INTER V 
REAL INTER 
RE/L LEADER 
REAL LUC 
LOGICAL SHELL 
RARAfETER SIZTkA- 100 
PARAfCTER SIZFTH=1G0 
CCH<W21/ 

* 

4 
* 

CCMOt/72/ 

♦ 

♦ 

♦ 

* 

ccwcw/n/ 

Integer detaul 

CIUENSICN LEADERS) 

DIHIN5ICN LlM>^2) 

DUtNSlOf hESSAG<ilCi) 

C1VENS1CN NA^E<23 ,1NSTk«Z0) 

Wihfc.Ni»JUN 

DATA SHElL/ ♦ 7ki£./ 

DATA AMIN/-2 .0/ 

Data amax/ij.O/ 

DATA STEFSZ/D.o/ 

DATA NAME/" J HAL0ACHV , CE5T/"MTT-GLD llOO m / ,FHXGl/ m 2146 m / 

DATA INSTR/'CNE HARD COPT - Cf>£ HlCRCFlLM - PLOT TO EOFILE " T-G2D 
*-/ 

DATA TRANGE/IOO,/ 

DATA« TEMPER^ I ) , 1=1 ,$IZrTH) /Sf ZFTH4900G./ 


AD SOT aS I ZTRA) , EXFCNaSl ZTRAJ , iNTERaSl ZTkA) , 
RESFCF*iS 1 ZTRA) , AESCRCal 0) , EXFCNCal 0) , 

V£ J GHTaS I Z7KA ) ,EXFOO*SlZ7RA) , EXPENDS l ZTRA) , 
ASSCROcSlZTRA) 

FRESSaSl ZFTH) , WATERaSlZPTH) , TEPF ERaS I ZP TH ) , 
HEICHTrtSlZPTH) ,HU+i5lZFTH) , ATCENTaSl ZPTH> , 

OTCENTasiZPTMJ , TRAKSPHSI ZPTH) .CELTRAaSlZPTH) , 
EFRESSaSlZrTHJ , CFftESSaS I Zr TH) ,ET£hCSaSl ZFTM) , 
CEFAULaSlZFTH) ,QZCNE*SI ZFTH) 

LIHIT « ANGLE, THETA, SHELL , LEVELS 
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001 4 7 35* DA7A<tOZCNC«I) ,1 = 1 , S I 2PTH) /Si 2PTH*4 4 67E~8/ 

00147 36* C 

00147 37* C COMPUTE RADIANS FROM ANCLES 

00147 3 8* C 

00151 39* RADaANG) -ANG/5? . ?9 57$ 

00 151 40* C 

00151 41* C CONFUTE SPECIFIC KUMIDITT FROM PRESSURE AtC VAPOR PRESSURE 

00151 42* C 

00152 43* RATMlX*V,P)=C22.*VAtF-ct.37d*V) ) 

001 52 44* C 

00152 45* C CCM=V7E THE ItsCfX CF REFRACT 1CN CF WATER- VAFOT AT THE TEMPERATURE 

00152 46* C AtC FKCSStflE CF A GIVEN LEVEL 

001 52 47* C 

00153 46* fiEFRAX«Ll)=l .♦«77.52G£-6)4FRES5‘xLn/TEFrER<iLl) 

00154 49* CAU- TAFLOTalHM) 

001 55 50* CALL I CF RMVaNAMC , OE ST » , INSTR ) 

00156 51* 100 CCKTINLE 

00156 52* C 

00156 53* C UT-in GCNTRCL CAROS 

OUt 56 54* C 

00157 55* REACTS, lt.EM^SOO) L IMIT .LABEL , FtACl , IHEAC2 ,HEAD3 

00157 56* C 

00157 57* C SET FLAGS TO AC JUST FREC IF I TATELE WATER CCPPUTATICN5 FOR 

00157 56* C REFRACTION AT HIGH ALTITUDES 

00157 59* C 

00166 60* lF<U_ACEL.EO.-SFHERC") SJ-tLL-.TR IE. 

00170 61* iF-LACeL.EO.-FLAT “ ) 5KTLL- -FALSE . 

00170 62* C 

00170 63* C ADJUST DISPLAY WlMXV CN TVE SC -4020 PLOTS 

out 70 64* c 

00172 65* IF«LABLL.EQ.“KAXCCT-) AHAV=L I H I T 

00174 66* IFaLACEL.EQ.-MlfJCO^-l AMIN=L I Ml T 

00174 67* C 

00174 68* C INPUT THE ANGLE OF OBSERVATION 

00174 69* C 

00176 70* IF*LABEL.EQ.'A1«LE “) ANGLE=L IHIT/IO. 

00176 71* C 

00176 72* C INFUT RESPONSE FlNCTlCNS 

00176 73* C 

00200 74* IF<tLAEEL.EQ. "RESPCN*) GO TO 200 

00200 75* C 

00200 76* C I7PUT RADIGSC^OE DATA 

00200 77* C 

00202 78* IF4LABEL.EQ.-WATEK *} CO TO 300 

00202 79* C 

00202 60* C IW=UT TEh~EkATURE RANGE CF 100 POINT TAELE 

00202 81* C 

00204 62* I F«L ABEL. EQ.' TEMPER") GO TO 400 

00204 63* C 

00204 C5* C 

002 05 66* XFtfLADEL.EQ. "►rSSAG-J GO TO 6GU 

00206 67* C 

00206 66* C CONFUTE A LIMB FUNCTION TO 50 DEGREES CENTERED AT A GIVEN NADIR 

DO 206 69 * C ANGLE 

00206 90* C 

00210 91* 


IF* LABEL .EQ. # LIK3 ■ ) GO TO 70D 



00210 

924 

C 


00210 

934 

c 

RECORD A PRESSURE INCREMENT TO ESTACLlSH THE STEP SITE FOR RADIOSONDE 

00210 

944 

c 

DATA FROM THE SURFACE TO CNE MILLIBAR 

00210 

954 

c 


00212 

964 


JFttLABCL .Etf .“STEP ■> STEFSZ=LlHlT 

00212 

974 

c 


00212 

964 

c 

O/ER-VJtl 7E MISSING RADIOSONDE DATA TO SATELLITE ALTITUDE WITH DEFAULT 

ft'*212 

994 

c 

DATA 

00212 

1004 

c 


00214 

1014 


lF«LAEEL*EQ.*CeFAUL*> CO TO 6CV 

00214 

1024 

c 


00214 

1034 

c 

IhTUT SPECIAL TRACE GAS TRANSMISSION FLNCTICN5 

0021 4 

1044 

c 


0021 6 

1054 


I FoLAEEL.EO. -TRACED GO TO 1100 

00220 

1064 


GO TO MM 

00221 

1074 

y.vj 

continue 

00221 

1064 

c 


00221 

1094 

c 

SET OBSERVATION ANGLE 

00221 

1104 

c 


00222 

1114 


THETA=RAD*ANGLE> 

00222 

1124 

c 

S 

00222 

1134 

c 

INITIALIZE PRECIPITATE WATER OFUTATlCN 

00222 

1144 

c 


00223 

1154 


LHXE-O. 

00223 

1164 

c 


00223 

1174 

c 

INITIALIZE EFFECTIVE PRESSURE G3K--UTATICN 

00223 

1164 

c 


00224 

1194 


EWHCLE=T« 

00224 

1204 

c 


00224 

1214 

c 

INITIALIZE EFFECTIVE TEMPERATURE CCPFUTAT ICN 

00224 

1224 

c 


00225 

1234 


TVM3_E-0. 

00223 

1244 

c 


00225 

1254 

c 

IhFLfT SURFACE RADlOSCfsCE DATA 

00225 

1264 

c 


00226 

1274 


REAC*5»12> PRESS*!) • DELFT * TEM^ER^l ) , HEIGHT* 1 ) , OZONE* 1 ) 

00235 

1264 


DEFAULT ) -0 

00236 

1294 


IFOHEAD2.M:. -KELVIN-) DEVFT-CEVFT+273 . 16 

00240 

1304 


lFa>CAC2.NE. -KELVIN-) TEMPER* 1 ) = TEMFEjR*1) +273. 16 

00242 

1314 


IFtiHEA03.NE*-OZC7£") 020^*1 ) -4 . G7E-3 

00244 

1324 


!F<iMLAC3.NE.-aZCNE “> FLCct2 , 1 » DEFALL*! ) ) -1 

00246 

1334 


DELFNT=CEW=T 

00246 

1344 

c 


00246 

1354 

c 

COMPUTE SURFACE VAPOR PRESSURE 

00246 

1364 

c 


00247 

1374 


VAP ER - VAPR ESc< CEVF T ) 

00247 

1364 

c 


00247 

1394 

c 

USE DEW-POINT GREATER THAN 77 DEGREES TO ASSUME 10 PERCENT 

00247 

1404 

c 

RELATIVE HUMIDITY 

00247 

1414 

c 


00250 

1424 


iFoiCEWFT.GE .350 . ) VAPCR= VAPkES*TEHPER* 1 ) ) 4 . 1 

00252 

1434 


IF*CELFT.GE.350.> FLC*3 , 1 jDEFAUL* I ) ) = 1 

00252 

1444 

c 


00252 

1454 

c 

COMPUTE THE SURFACE SPECIFIC HUMIDITY 

0O2S2 

1464 

c 


00254 

1474 


PART=RATMIX*VAFER »FRE5$*1 } ) 

00254 

1494 

c 
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00254 

1494 

c 

STCRE THE SPECIFIC HUMIDITY COMPUTED AT THE SURFACE 

002*4 

1 50* 

c 


00255 

1514 


HUFVtl) ^PART 

00255 

152* 

c 


002*5 

153* 

c 

THE AIUDS-CCNTIFCTCRS q f CAK&ZN DIOXIDE A*C OZCNC AT THE SURF ACE 

002*5 

154* 

c 


GOZ 56 

155* 


ATCENT<*1)=0* 

00257 

1 50* 


OTCENTal )=0* 

00257 

157* 

c 


00257 

15 a* 

c 

THE EFFECTIVE PRESSURE AT THE SURFACE IS EOOAU TO THE RACIC6CICE 

00257 

159* 

c 

SURFACE PRESSURE 

00257 

160* 

c 


GQ2GL) 

161* 


EFRESS'tl)=FRES5'tl> 

00201 

162* 


CFkESS*1)=FRESS*1) 

00261 

163* 

c 


00261 

164* 

c 

LAEEL THIS FACE CF AE&^BINC CAS CONCENT RAT I CNS. 

00261 

165* 

c 


00202 

166* 


LRITE«6,27) tiFCSSAG*!) ,I = 1,1Q| 

00270 

167* 


CEVJPT =DEVJFT-273 .16 

00271 

168* 


TEK^2F - TEFT CRd 1 ) - 2 73 . 1 6 

00272 

169* 


VR I TE<*6i 14) FRESSal) ,DEVFT,TEM-^ 

00277 

170* 


DO 1G 1=2 |LlMl T 

00277 

171* 

c 


00277 

172* 

c 

EXTRACT A LEVEL'S DATA 

00277 

173* 

c 


00502 

174* 


REACTS ,12> PRESSED , CEVR T , TEFFER* I ) .FEICHT*!! ,OZCF£*I> 

00311 

175* 


C€TAU_a| ) =0 

0031 1 

176* 

c 


00311 

177* 

c 

CONVERT FEIGHT TO FEET IF FtEDED 

0031 1 

17 a* 

c 


0Q5 12 

179* 


IFoM^ACI. EO. -UTTERS-) PC I GHT* I ) =*IE I GHT<x 1 )*3.2SG8 

00312 

160* 

c 


00312 

161* 

c 

COT-UTTE AVERAGE PRESSURE CF A LAYER 

00312 

162* 

c 


00314 

183* 


PRESUR=*FRESS<s I ) +FRESSa I -1 ) ) /2 . 

00314 

164* 

c 


00314 

165* 

c 

CONVERT FROM CENTIGRADE TO KELVIN IF NEEDED 

00314 

166* 

c 


00315 

187* 


XRiFCAC2.7£. “KELVIN") CEW=T^C£W-'T*273 . 16 

00317 

186* 


IFCXHEAC2.FE. "KELVIN") TEFPER*! ) =TEFrERctI ) +273.16 

OLt517 

189* 

c 


00317 

190* 

c 

1FPOSE A CONSTANT OZCTt FRCF1LE IF NEEDED 

0031? 

191* 

c 


00321 

192* 


iFdHEAC3.«:.*o2a£*> ozaTi^n =4.67E-d 

00323 

193* 


I F*HEAC3 . KE . "OZCrE ") FL£*2 , 1 .DETAUL*! > )=1 

00323 

194* 

c 


00323 

195* 

c 

Commute the average layer CEu point 

00323 

196* 

c 


00323 

197* 


DEVA VG DCVFNT * DELFT? /2. 

00325 

198* 

c 


OG325 

199* 

c 

CCMPUTE TEE VAPOR F RES SURE 

00325 

200* 

c 


DO 326 

201* 


V AP ER = VAFRE S«t C EV1A VG ) 

00326 

202* 

c 


00326 

203* 

c 

USE COW- POINT GREATER THAN 77 DEGREES TO ASSURE 10 PERCENT 

1X1326 

204* 

c 

relative humidity 

00326 

205** C 
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OOJ27 

2064 


1 F-CCWPT . CC . 350 . ) VAP CR = VAP RES'iTCHPCR'tl ) > * . 1 

00311 

207* 


1F«CEVA T.0C.350.) FUD<i3 f l ,CEFAUL<* 1 > ) M 

0033 1 

296* 

C 


00331 

209* 

c 

CCW-UTE THE Ml XI NO RATIO CF THIS LAYER 

00331 

210* 

c 


00333 

211* 


HUMIC-G22.*VAPER/«FkC$UR-VAFER) 

00333 

212* 

c 


00333 

213* 

c 

COMPUTE THE AVERAGE TEMPERATURE CF THIS LATER 

00333 

214* 

c 


00334 

215* 


TCMP0ft=*TE7eERaI ) ♦TEMPER ** I -1 ))/2* 

00334 

2164 

c 


00334 

217* 

c 

CCMFUTE TfC HEIGHT CF THE TOP CF THIS LAYER IF TCECEC 

CV*334 

216* 

c 


00335 

219* 


lFaHElGHTH) .LE.EO 

00335 

220* 


*te I GMT 'll > =HE I GHTti l -1) ♦ 29 . 3*TCPFCR*ALOGaFRESS'i I - 1 > /PRESS'* I ) 1 

00335 

221* 


*** 1 . ♦ .0.7061 *MUM 1 C ) *3 .2606 

00335 

222* 

c 


00335 

223* 

c 

STORE TVE ST-CCIFIC HUMIDITY AT EACH LAYER TOP 

00335 

224* 

c 


00337 

225* 


HUM* I > =RATMI XdVAFER , FRESU*) 

00337 

226* 

c 


00337 

227* 

c 

COFUTE TVe PkECI PI TABLE WATER AT THIS LEVEL 

00337 

226* 

c 


00340 

229* 


VM TERtf I ) TTRa I - 1 ) * <*> l > *SP EEREaS UT , 1 ) RESSa I - 1 ) -PR ESS^* I ) ) ) / 

00340 

230* 


* 9 80. 

00340 

231* 

c 


00340 

232* 

c 

OOFLfTE The ATMOS-CEXmeTERS CF CARDCN-ClCfclCG AtC QZCfC TO THIS 

00340 

233* 

c 

LAYER 

00340 

234* 

c 


□0341 

235* 


ATCENT<*1 ) -aFRESS^i 1 -PRESS* 1 M * . ?EO 

00342 

Z36* 


OTCENT«*I ) ^OTCENT^I -1 ) ♦nFRESS* I -1 } -PRESS* I ) ) **OZCK£r<i 1 1 ♦OZi^E* I -1 > ) * 

0034? 

237* 


* 237.968 

00342 

236* 

c 


00342 

239* 

c 

COFUTF TIC EFFECT I VE TEMPERATURES AND PRESSURES USED TO RELATE 

00342 

240* 

c 

hOCOMOUS PATH TRANSMISSION FUNCTIONS TO IMOlXErCCUS SLANT 

0034? 

241* 

c 

PATHS THROUGH THE ATMXFICRE 

0034? 

242* 

c 


00343 

243* 


EW<XX=EU*XE +«UA TER* I ) -WATER* I - 1 )) 4FRESLR 

□0344 

244* 


7VJ-CLE=TVi*XE+t*WATER« I ) -WATER* I - 1 ) ) *TEPFCR 

00345 

245* 


CFRESS^I ) =dFRESSal ) ♦FRE5S* I) ) /2. 

00346 

246* 


0=R£SS*I ) “EV4-CLE /WATER* 1 7 

00347 

247* 


ETT>eS*I ) =TVHCLE/WATER^ I ) 

00347 

246* 

c 


00347 

249* 

c 

COFl/TE CENTIGRADE VALLES FOR DISPLAY 

00347 

250* 

c 


00350 

251* 


CeWAVG=C€WPT-273 . 1 6 

00351 

252* 


m«=OR-TT>FERctI ) -273 . 1 6 

00351 

253* 

c 


003 3 i 

254* 

c 

PRINT lUTERf'EDI ATE RESULTS CF THIS PROCESS 

03352 

25C* 

c 

VR ! TE^C, 131 FRCSS^I) ,CEWAVG$ TO FOR , HE iCnTai ) , WATER <* i t ,ATCENTaI / , 

0035? 

257* 


* OTCENT»in ,EFRESS^n , CFRESS* II , ETEKF5*I > 

00366 

256* 


CCVPNT-DEVJPT 

00367 

259* 

10 

CONTINUE 

00371 

260* 

11 

FCRMAT<i 1 19 , A6 , 6X , 3A6I 

00372 

261* 

12 

FCRHAT«*1 X i F7. 1 t 2Xf F6 . 1 j 4X ,F7* 1 ,3x ,F1G. 1 ,ElCi.5> 

00373 

262* 

13 

FOkMAT<i 1 X,F7„ 1 |2X|FG.l ,4X r F7,l t 3X,FlD.l ,1X ,F7.4 , lX,F6. i ,F6.3 ,2X, 
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M3 7 3 

263* 


*F6.l, 

M37 3 

264* 


♦6X.F6.1 ( 6X,F5.n 

MS 7 4 

2654 

14 

FORMAT*' PRESSURE CEW POINT TEMPERATURE HEIGHT PREC1P ATWD6-CE 

00374 

266* 


*NT EFFECTIVE CONSTANT EFFECTIVE" , “ * ,F7 .1 t 2X,F6.1,4X,F7.i 

M3 74 

267* 


*,• FEET WA.ER CC*2 O ZONE FRESSUX E CONCENTRATION TEMP 

M3 74 

26*4 


4ERATURE-) 

M 375 

269* 

15 

FCRMATuF7.1 ,F7 *4) 

M370 

270* 

16 

FCRMAT*1X,F7.1 ,F7.4) 

00377 

271* 

17 

FCRHATaTN ANGLE ,F4.1 f 20H DEGREES FRCH NADIR ,/,21‘iH SURFACE TEMPER 

0037 7 

272* 


♦ATURCjFT . 1 ) 

00400 

273* 

16 

PCRMATalX,36HTE^ERATURE RANGE CF 100 POINT TABLE) 

00401 

274* 

19 

FORMAT^ 42H1WAVENUMCCR MIDPOINTS Arc PERCENT RESPONSE* 

00402 

275* 

21 

FORMAT* 3 7H 1 TEMPERATURE RANGE CF 100 POINT TAELE) 

00403 

276* 

22 

FORMAT* 5 1 HOP RE C I F HEIGHT SPECIFIC TEhP CORRECTION,/, 

00403 

277* 


*36H WATER HUMIDITY SCONCING,/, 

00403 

276* 


*42H CM FT CM/KCM K K) 

00404 

279* 

23 

FCRMAT*1H ,3A6) 

00405 

260* 

24 

FcuMAT*tF4 .1 1 1 AH DEG. TO NAG Ik) 

00406 

261* 

25 

FCRMAT«F5.1 ,13H SURFACE TEPF) 

00407 

2*2* 

26 

FORMAT* 10A6) 

0O4I0 

263* 

27 

FCRMATalHl.ir.AS,/) 

00411 

264* 

26 

FORMAT *4 El 0.5) 

00412 

265* 

29 

FORMATS X,F6.1 , |F6* 1 , 4*3X ,F7 . 4 ) ) 

0L1413 

266* 

31 

FCRMAT<tl H.) ,3AG i / , “ uAVENU-GER A**ATMDS-CENT*=X) **EF**C) • ,/, 

00413 

267* 


*- INTERVAL. ABC VEIGHT") 

aw 14 

266* 

32 

FORMAT* *0 INTERPOLATED Lire FUNCTION TO “,F7.Ct ( “ FEET *,FS*1, 

00414 

269* 


*- DEGREES FROM NAG JR" , Z<U , 1 X , 2l*F6 . 2) ) J 

00415 

290* 

33 

FCRMAT*1 X i F7 . 1 fZXpElO.5 »2X ,F6. 1 ,2X,EHi.5> 

01416 

291* 

35 

FORMAT*" 1 DEFAULTED STRATOSPHERE PROFILES" ,/, 

00416 

292* 


* " FRESSUTC OZOZ TE>? SPECIFIC-,/, 

00416 

293* 


* " SOLfCJNC HUMIDITY-, /, 

00416 

294* 


* ■ re GH/KGM K GK/KGM") 

00417 

295* 


CALL GRIDlVd.1 ,-50. ,50. , AMIN, AMAX, 1 0. , .5 ,0,0. 1 , 1 ,4,4) 

00420 

296* 


CO TO 100 

00420 

297* 

C 


0042(3 

296* 

c 

INPUT RESPONSE FUNCTIONS AS PERCENTAGES AT MIDPOINTS CF EACH 

00420 

299* 

c 

INTERVAL 

00420 

300* 

c 


00421 

301* 

zoo 

CONTINUE 

00422 

302* 


LEADERS) = HEAD 1 

00423 

303* 


LEADERS) =HEAD2 

00424 

304* 


LEADER rf-3 J "HEAC3 

00425 

305* 


*RITE*6,19) 

00427 

306* 


DO 30 1=1,65 

D04 32 

307* 

30 

RESPCNaI>=0* 

00434 

306* 


DO 40 I=J pLlMlT 

00437 

309* 


READ*5,15> INTER V.FUNCTI 

00443 

310* 


VwR I TEfltS ,16) INTER V , FUNCTI 

0044 7 

311* 


J=«t INTERV-1 2 . 5) /25 . 


*1^* 

4U 

RLSrOHJ) "FUNCT 1 

00452 

5 1 3* 


GO TO 11*0 

00452 

314* 

c 


00452 

315* 

c 

INPUT THE TEMPERATURE RANGE CF THE ©LACK-GCCY INTENSITY/ 

00452 

316* 

c 

TEMPERATURE SCALE IN CENTIGRADE 

00452 

317* 

c 


00453 

316* 

400 

REACTS, 15) TRANGE 

00456 

319* 


W7lTEa6,S0) 


\ 


NEW 

pew 

NEW 

NEW 
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004 60 

320* 


WUTE*6»16) TRANCE 

W463 

321* 


CO TO 100 

00464 

322* 

300 

coni iNue 

00465 

323* 


call pltnC'iG) 

004 66 

324* 


CALL FLTERM 

0046? 

325* 


aw 

00467 

326* 

C 


00467 

327* 

c 

IfFUT Cl SPLAY LABELS 

00467 

326* 

c 


00470 

329* 

600 

REACTS , 26) <t MESSAGED ,1 = 1,10) 

00476 

330* 


CO TO too 

004 76 

331* 

c 


004 76 

332* 

c 

FROCUCE A Llffi FLNCTlCN DISPLAY TO 50 DECREES CENTERED AT THE 

00476 

333* 

c 

GIVEN NADIR ANCLE 

00476 

334* 

c 


00477 

335* 

700 

CONTINUE 

00500 

336* 


ftlCHS=LlHlT 

00501 

337* 


lF*STEFSZ*GT.n.) CALL STEFEK 

00501 

338* 


lF«tSTEFSZ.GT .0.1 CALL CISTLY 

00501 

339* 

c 


00503 

340* 

c 

establish a tetter at ure/elack&xy intensity scale 

00503 

341* 

c 


00505 

342* 


CALL TAETE>taTE>rER*n .TRANCE) 

00506 

343* 


1=1 

00507 

344* 

710 

1 = 1*1 

00507 

345* 

c 


00507 

346* 

c 

DETER Hire IlsCEX to interpolation levels from the requested height 

00507 

347* 

c 


00510 

346* 


IFrtI .GT.10D) CO TO 730 

0051? 

349* 


IF*HET6HS,GE.»^:iGHTa!-t) .ATT^IGHS.Le.fOGHT^m GO TO 720 

00514 

35T«* 


tFoTtlGHTal) .L£.f<*> CO TO 730 

00516 

351* 


CO TO 710 

00517 

352* 

720 

HIGH=HEIGHS 

00520 

353* 


CO TO 740 

□□520 

354* 

c 


00520 

355* 

c 

DATA UNAVAILABLE FOR THE REQUESTED iNTERrCLAT ICN 

CO 5 20 

356* 

c 


00521 

357* 

T3P 

HICH=fE:iGHT<iI-l) 

00322 

358* 

740 

IHIGH=I-1 

00523 

359* 


DO 750 K=l,2l 

00523 

3€0* 

c 


00523 

361* 

c 

COTUTE THE LIKG ANCLE FkCH NADIR IN 5 DEGREE STEPS 

00523 

362* 

c 


00526 

363* 


LITC^K) =ANGLE-55 . *K*5 . 

00527 

364* 


THETA^RADciLlfOiK) ) 

00527 

365* 

c 


00527 

366* 

c 

CONFUTE THE LCVLk LAYER TEKTEKATURE 

00527 

367* 

c 


00530 

366* 


LIHIT=IHIGH 

ZZZZi 



V-f h 1 i -v? 

00532 

370* 


CALL RASf*CC«47ErT 1 J 

00533 

371* 


1F*LIMIT.EQ.1> TEfFl =0, 

00533 

372* 

c 


00533 

373* 

c 

COMPUTE THE UPPER LAYER TEMPERATURE 

00533 

374* 

c 


0053 5 

375* 


LINITMMIGHM 

00536 

376* 


CALL RADMCC^TCMF?) 


few 


h $ 6 § 5 
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00536 

377* 

C 


00536 

376* 

c 

INTERPOLATE A TEMrEfiATUlE TO THE REQUESTED HEIGHT 

00536 

379* 

c 


00337 

360* 

T30 

LltEWK+21) =TERFttHElGHT«LtMlT-l) ,MCl CHIul. llllT) , TEMPI ,TEMP2 ,HICH> 

00341 

361* 


WUTE<t6,32) HlGH,ANGLE,<iLtH><K> ,K=1 ,42) 

00551 

362* 


XP=-50. 

00532 

363* 


YP^O. 

00353 

384* 


DO 760 1-1,21 

00556 

365* 


X^LUf^l ) 

00557 

386* 


YsLlfOiI*?!) 

00560 

367* 


CALL LlSCVftNXVaXF) pNYVaYF) t NXVaX> ,NYV*Y)) 

oosei 

366* 


XP=X 

IX) 562 

369* 

760 

YP = Y 

00564 

390* 


GO TO 100 

00565 

391* 

600 

CONTlMJE 

00365 

392* 

C 


00565 

393* 

c 

READ DEFAULT PROFILES AND INTERFCLATC MISSING DATA 

00565 

394* 

c 


00566 

395* 


iFaSTEFSZ.GT.O) CALL 5TEFER 

00570 

396* 


WUTE*6p35> 

00572 

397* 


REACTS ,33 , EM>li>J) PRESS- .QZCNEP , TEfFRE , HUM I CP 

0L> 600 

396* 


W<ITE*6,33> PRESS- pOZOEF , TEKFRE , HUMICF 

00606 

399* 


OO 620 J=2, LIMIT 

00611 

400* 


R£AOi5,33 ( E7CrlOO> PRESUfi .QZCPCS , TENFCfi , HUMIC 

00617 

401* 


V«lTEaG,33) ESU< t OZa«, TENTER* HUN I C 

□0617 

402* 

c 


00617 

403* 

c 

SCAN REQUESTED LEVELS FOl DEFAULT PRESSURE WINDOW 

□0617 

404* 

c 


00625 

405* 


DO 810 1=2, LEVELS 

00630 

406* 


Irul-RESSiI) .LT.i-RESSP.ArL:. PRES^r I) .CE.PRESUM) GO TO 815 

00632 

407* 


GO TO 610 

00633 

406* 

815 

CONTINUE 

00633 

409* 

c 


00633 

41D* 

c 

INTERPOLATE MISSING TEJFEKATLRE 

00633 

411* 

c 


□0634 

412* 


IFcfTEfFERtl) .GT. 1000.) TEMPER* I) =TERFaFRESSP ,PRESU*,TEMPRE, TEfFCR, 

00634 

413* 


* FnESSal ) ) 

00634 

414* 

c 


□0634 

415* 

c 

INTERPOLATE MISSING SPECIFIC HUMIDITY 

00634 

416* 

c 


00636 

417* 


IF«iFLCa3 ,1 .DEFAULT!) ) ,EQ*1) KUMal ) =TERP*PRES5F ,PfiESlR ,HU^ICF ,HU*HlO 

00636 

416* 


* f FRESS*IJ) 

00636 

419* 

c 


□0636 

420* 

c 

INTERPOLATE MISSING QZCfE 

00636 

421* 

c 


00640 

422* 


IFaFLC«2 , 1 p CEFAULal ) ) *E9 . 1 ) OZCTC^I ) =TERF*FRESSP , FREStE ,OZCfEP , 

00640 

423* 


* CCOCS, PRESSED) 

00642 

424* 

810 

CONTINUE 

00642 

425* 

c 


00642 

426* 

c 

FETCH NEXT LEVEL DEFAULT PARAMETERS 

00642 

427* 

c 


00644 

426* 


PRESSP=FKESUR 

00645 

429* 


OZCWEP^OZONES 

00646 

430* 


tempre-temfcr 

00647 

431* 


humidp=humid 

00650 

432* 

820 

CONTINUE 

00652 

433* 


GO TO xoa 


NEW 
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00632 

434# 

C 

00652 

435# 

c 

ooesz 

436# 

c 

00652 

437# 

c 

G0652 

436# 

c 

00*52 

439* 

c 

W65Z 

440# 

c 

00632 

441* 

c 

00652 

442# 

c 

00652 

443# 

c 

00452 

444* 

c 

00652 

445# 

c 

00652 

446# 

c 

00652 

447# 

c 

OOG53 

446* 

1100 

00654 

449* 


00661 

450* 


00664 

451* 


00666 

432* 


00674 

433* 


00675 

434* 


00676 

455* 


00677 

456* 


00700 

437* 


00710 

436* 

uta 

00712 

459* 


00713 

460* 


00713 

461* 

c 

00713 

462* 

c 

0U713 

463* 

c 

00716 

464* 


00720 

465* 


00723 

466* 

to 

€0732 

467* 

11 

00732 

466* 


00732 

469* 


00733 

470* 

12 

00734 

471* 


00735 

472* 


00736 

473* 


00796 

474* 

C 

00736 

475* 

c 

00736 

476* 

c 

00736 

477* 

c 

00736 

476* 

c 

00741 

479* 


00742 

480* 


00743 

48J* 

c 

00742 

402* 

c 

00742 

463* 

c 

00744 

464* 


00744 

465* 

c 

00744 

466* 

c 

00744 

467* 

c 

00744 

466* 

c 

00745 

469* 


00745 

490* 

c 


INPUT TRACE CAS TRANSMISSION FUNCTIONS WHERE 

WEIGHT TRUNCATED TO AN INTEGER IS AN INDEX TO THE WAVENOCER 
INTERVAL FCft WHICH THIS ANALYTIC FUNCTION AFPLYS 

WEIGHT FRACTION IS THE GEOMETRIC KEAN WEIGHTING NIDIFIED BY THE 
PLANCK IAN FUNCTION FOR THE WAVE-NUCER INTERVAL FRACTION 
CF THIS ANALYTIC FUNCTION 

A6SCR.O = A ThCSe ARE PARAMETERS APPEARING IN ANALYTIC 

EXPCMD = B FUNCTIONS CF THE FORM A*^TMX-CENT**B) *«£P**0 

EXPOO s C FOR THE INTERVAL AhC WEIGHTING DESCRIBED AECWE 

CONTINUE 

WUTE«6,31> hCADl ■ HEAB2 , hEA03 
DO 1110 1=1, LIMIT 
IF#I *CT .100) GO TO 100 

REAO*3,2e,E7C=l001 VCIGHUI) ,AB$C«OH> pEXPCM*!) , EXF<>Oi! > 

IIsVCICHTaU 

WElC«=V£lGm*I)-II 

BEG IN= INTER41I) -1 2 . S 

FIN! S= INTER*! 1> +12.5 

W«ITE<t6,29) BE&INpFINlS.ABSCROin .EXPOO*!) .EXFCNDuO ,V£IGH 
CCNTINUE 
GO TO lOO 
SLGROUTUC CISFLY 

THIS SUBROUTIWT DISPLAYS T>C MX I FT ED DATA BASE 

WRlTEaG.il) 

DO IQ 1=1 (LEVELS 

WclTE*6,12) PRESS*!) pHUMal) , TEMPER*!) ,O20hC*I) 

FORMAT*- 1FRE5S SPECIFIC TE>r- OZ<>£“ ,/. 

♦ - HXIDITY SOLhClNG-,/, 

* • MB GM/KGM K GM/KGM - ) 

FORMAT*! X»F7 . 1 t lX>ElO«5»2X|F6*2, lXpElQ.S) 

LIMIT=LEVELS 

RETURN 

FUNCT ICN SPHERE**, L> 

THIS IS A CONDITIONAL FUNCTION WHICH KXIFIES THE PRECIPITABLE 
WATER AT SATELLITE ALTITUDES FOR TIC SPHERICAL SHAPE OF THE 
EARTH 

SPHCRE=t. 

IF*. NOT .SPELL) RETURN 

TEST IS THE SIMC CF TIC ANGLE CF OBSERVATION TO THE HORIZON 
TEST= 2089 6696. /*POfl 96 696. ♦HClGHTetU > 

NO COMPUTATIONS ARE MADE WHEN THE FIELD CF VIEW IS ABOVE THE 
HORIZON 


IF*THETA.CT.ASIN*TEST> > RETURN1 
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0074 5 

491* 

C 

mow is THE AVERAGE HEIGHT CF THE WATER VAPOR LATER 

0074 5 

492* 

c 


00747 

493* 


H 1 l CHTtiL ) ♦ HU 1 CHT it- 1> > /2 . 

00747 

494* 

c 


00747 

495* 

c 

ACC the average height c f the later to the radius cf the earth 

00747 

49G* 

c 


00750 

497* 


PARTI = HIGH ♦20698 696, 

00750 

496* 

c 


00750 

499* 

c 

MULTIPLY CY THE ►CAN INDEX CP REFRACT lCN AT THE AVERAGE HEIGHT 

00750 

500* 

c 

CF THIS LAYER 

170750 

501* 

c 


00751 

302* 


FAftT2=FARTl**RCFRAX*L) ♦REFRAX<4_-1 ) ) /£. 

00751 

503* 

c 


00751 

504* 

c 

MULTIPLY Or THE lhCEX CF REFRACT I CN AT THE 1CIGHT OF TFE SENSOR 

00751 

505* 

c 

AFC BY THE SIMI CF THE ANGLE CF OBSERVATION 

00751 

506* 

c 


00752 

307* 


FART 3=F Aft T 1 *REFRAX*L) *SlNl HCTA) 

00755 

506* 


FART4=5*2RT««FART2*FART2) -^FART3*FART3) ) 

00755 

509* 

c 


00753 

310* 

c 

COMMUTE THE MIAN FUNCTION ASSUMING W CE SERVING SENSOR IS JUST 

00753 

511* 

c 

ADlVt THE MIAN LAYER TO WHICH IT IS APPLIED 

00753 

512* 

c 


00754 

51 3* 


SPHERE =FART2/F Aft T4 

00753 

514* 


RETURN 

0G756 

515* 


FUNCTION TERP<xXl,X2 t TI,Y2 ( X) 

0CJ75G 

516* 

c 


00756 

517* 

c 

THIS FUNCTION PERFORMS LAGftANGlAN iNTERFCLAT ICN 

00756 

Sift* 

c 


00761 

519* 


SLCFEA=«0<-X2) /*X1-X2) 

00762 

5ZO* 


$LCFED=*X-X1 ) AtX2-Xl ) 

00763 

521* 


TtKP= SLOPE A*Y1 *SLCFEC*Y2 

00764 

522* 


RETURN 

00765 

523* 


SLEROLnim CPATHS 

00765 

524* 

c 


00765 

525* 

c 

THIS SUBROUTINE COMPUTES ANGLE DEFEMJENT OPTICAL PATHS 

00765 

526* 

c 


00770 

527* 


ELKLE^O, 

00771 

526* 


TV«CLE=0, 

00772 

529* 


DO lO 1=2, LIMIT 

00775 

530* 


OTCENT* I ) =OTCEHT«t I -1 > ♦aFRESS* I -1 ) -PRESS* I) ) **OZCNEaI > +OZCFO* I -1 ) ) * 

00775 

531* 


* 237.968 

00776 

532* 


WATER* I) =WATERa I - 1 ) +*HUMiI ) *SFH£RE*a 1 0 , 1 ) **FRESS*1 ~i ) -PRESS* I > ) ) / 

00776 

533* 


* 980. 

00777 

534* 


FRESUR=*PRESS*I ) +PRESS*I -1 ) ) /2 . 

OlOGO 

535* 


TE>POK=ci TEMPER* M ♦TEMPER*! -1 )) fZ . 

OlDOl 

536* 


ATCENT* I ) =«FR£5Sal ) 'PRESS* I > ) *. 260 

01002 

537* 


ELHa_E=EVHCLE+*WATER«n -WATER«I-1) ) *PRESUR 

01003 

532* 


TVJ^CLC^TW -ILC < tiWATER* I ) -WATER** t — IT/ *Tc MFCfi 

01004 

539* 


PPRETSSrtl) F/WATFRtftl 

01005 

540* 


ETEM=S* I) =TLHCLE/WA TER* 1 ) 

01006 

541* 


IF*I ,GT *1) 

01006 

542* 


*HEIGHT*l) =HEIGHT* 1-1 > *29. 3PTEK''OK*ALCG*FRESS*I'i> /PRESS*!)) 

01006 

543* 


***1 .♦.OQCi61*HUM*I>> *3.2808 

PIOIO 

344* 

10 

CONTINUE 

01012 

545* 


RETURN 

01013 

546* 


SUDnOUTI NE STEFER 

01013 

547* 

c 



-05 
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01013 

34*4 

C 

THIS SUBROUTINE EXPANDS CR CONTRACTS THE CATA TO THE NEEDED STEP 

OU‘il 3 

5494 

c 


ci me 

550* 


JM 

01017 

551* 


LCClNT-HlmrRES5al)/STCFSZ,SIZFTH) 

01017 

552* 

c 


01017 

353* 


TEST FOR END CF AVAILABLE CATA 

01017 

5544 

c 


01020 

555* 

5 

IF*J.6T,LCCINT> RETURN 

01020 

556* 

C 


0102*3 

557* 

C 

SET AVAILABLE DATA 

01020 

55 0* 

C 


01022 

559* 


FRCSSF-FRESSrtJ} 

01023 

5611* 


PRESLR“FRESSn J*1 ) 

01023 

561* 

C 


□ 1023 

562* 

c 

SET WANTED DATA 

01023 

563* 

c 


01024 

564* 


TEVT RE^TEFPER^J) 

01025 

565* 


TEM'OR = T t>F EK <i J ♦ 1 ) 

0102& 

566* 


OZCFEP^DZaCttJ) 

01027 

5G7* 


OZC*CS=OZCNE*J* 1 ) 

01030 

566* 


KJMlDF=lfUM4j) 

01031 

569* 


HUM I D- HLMn J ♦ 1 ) 

01CT32 

570* 


PREWAN^rRESS'U J -STEPSZ*J 

01032 

571* 

C 


01032 

572* 

C 

is wanted data within available data whccw 

01032 

573* 

c 


01033 

574* 


IFn(FREVi^N.LT.PRES5F.ANC,FREWAN.6E.FRESUx) W TO 10 

01035 

575* 


CO TO 40 

01036 

576* 

10 

CONTINUE 

01036 

577* 

c 


01036 

376* 

c 

HCW MANY LEVELS CF WANTED DATA ARE WITHIN W1MXW 

01036 

579* 

c 


01037 

580* 


LCC S 1 2=<tFRE SSF - FRE SLR J /STEF SZ 

01037 

561* 

c 


01037 

582* 

c 

HAKE SPACE IN AVAILABLE DATA FOR WANTED CATA 

01037 

583* 

c 


01040 

564* 


DO 20 K-LEVELS, J,-l 

01043 

565* 


FRESSrtK+LCCS I Z> -F RESS^tK ) 

01044 

566* 


TEM^ EK ax +LCC5 1 2 ) = TOR Ek*K) 

01045 

567* 


CZCNE*K*LCCSIZ) =OZCfsEaK) 

01046 

506* 

20 

HUM*K+LCCS 1 2) =HUtHK) 

01050 

569* 


level s= level s+ LCC S I Z 

01050 

590* 

c 


01050 

591* 

c 

LOAD CATA OVER WANTED INTERVAL 

01050 

592* 

c 


01051 

593* 


DO 30 K-l ,LOCSlZ 

01054 

594* 


j=j+i 

01053 

595* 


7EMF ER« i > -TEi\F»iF*”ES5F * FRESUv , TOFTX , TLKXR ,rRC\XAN) 




CaZs-4 ~ * 9 4 w 1 Wfc ’* ‘ 9 " *" • k - w ^ ■ 1 

viCJ * 

597* 


u , , 4 . ( r aTit E SiP p prrr yn $ * j £p # HLH 10 § FfiEVJAM) 

01060 

596* 


FRE5S-J J) =FREWAN 

01061 

599* 


PREWAN=FREV‘AN-5TEFSZ 

01062 

600* 


C£FAUL*J >=CEF AUL<tJ-l> 

01063 

601* 

30 

CONTINUE 

01063 

602* 

c 


01063 

603* 

c 

RESET BOTTOM LEVEL CF AVAILABLE DATA 

01063 

604* 

c 




010G5 
01066 
Ol 066 
010 GG 
01066 
01067 

01072 

01073 

01074 

01075 

01076 

01077 
OllOl 
Gl 102 
01103 


FREP 

RJRiHJR 


605* 


GO TO 5 

606* 

40 

CC«T I NUC 

607* 

c 


60S* 

c 

DESTROY TCF LEVEL CF AVAILABLE DATA 

600* 

c 


610* 


DO 50 6= J , LEVELS 

611* 


FRESS'iKM ) JT-KESS'tloZ) 

612* 


TEMTERoK* 1 > -TEMi“EE«tR*2) 

613* 


C2CMEaK+ 1 ) =<^':i.E«K+2> 

614* 


= Ht*VJ*+2> 

615* 


CEFAUL*aX*‘l) =CCFAU_aK+2) 

616* 

50 

CONTINUE 

617* 


LEVEL S-LEVFLS-1 

616* 


GO TO 5 

610* 


END 


DC CF CCM 3 lLAT ION ' NO DIAGNOSTICS. 


IFF*. 

0256 ^ 03 / 20 - 16' 51 
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FOfi.S MAIN/SLANT ,KAtN,MA|N/SLANT 
FCA S*L-Ci*/20-16’50 «0,0> 


HA IN PROGRAM 

STORAGE USED' CCCEul) 002555 CATA<iG) 000750 CLANK CCHMCt*l2) IW'V.V.1 
CC#^:w CLOCKS ' 


000.3 

Z1 

CO 1 464 

0004 

22 

□15550 

0003 

Z3 

axv;o5 


EXTERNAL REFERENCES «CLCCK, NAfC) 


0006 

VAFRCS 

0007 

taflot 

ooio 

ICFRMV 

OOll 

TABTEM 

O.i 12 

SETMIV 

0013 

GR1C1V 

0014 

FR INTV 

0015 

L1IFIV 

0016 

NXV 

0017 

NYV 

0020 

MvJCCCS 

0021 

APR NT V 

0022 

CHSIZV 

0023 

R1TE2V 

0024 

RADMZO 

0025 

fltnc 

0026 

PL TERM 

0027 

NINTH* 

0030 

NTDU3 

0031 

NICC$ 

0032 

NUCUS 

0033 

NlOl$ 

0034 

ALOG 

0035 

NS TCP $ 

0036 

AS IN 

0037 

SIN 

0040 

SORT 

0041 

PeRK4$ 

0042 

*CRR3$ 


STORAGE ASSlCNf-tENT tiOLOCK . TYPE* RELATIVE LOCATION, NA^F) 


DOOl 

00D6C.Q 

1C1L 

0001 

D02252 

10L 

0001 

002367 

10L 

0001 

000013 

idol 

0001 

002042 

1 00 4 G 

0OQ1 

002143 

10360 

0000 

000465 

1 IF 

0000 

000150 

1 IF 

0001 

001603 

1 1 00 L 

0001 

002414 

U05v 

ooot 

0*72431 

111 GO 

0001 

002505 

1 1 346 

0000 

00051 1 

12F 

oraa 

0001 53 

12F 

0000 

000162 

13F 

OOOQ 

000200 

14F 

0000 

OOG243 

15F 

0000 

000245 

16F 

GOOD 

000250 

17F 

0000 

ooo2r 4 

1<?F 

oooo 

0002 74 

19F 

DOOl 

001250 

200L 

0000 

000304 

21F 

0000 

000313 

Z2F 

0000 

000351 

23F 

OOOQ 

000353 

24F 

0000 

00*3357 

25F 

0000 

000363 

26F 

0001 

000271 

264 C 

0000 

000364 

27F 

OOOO 

000366 

zer 

0000 

0003 70 

29F 

0001 

000315 

300G 

0001 

DOOl 32 

3G0L 

0000 

000376 

31F 
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0000 


000422 

33P 

0000 


000430 

35F 

□Dili 


000603 

390L 

0001 


DO 2 500 

40L 

0001 


001330 

400L 

OOUl 


000631 

427C 

0001 


000722 

454G 

DtV.ll 


01X1767 

4646 

0001 


002312 

5L 

x.voi 


0111352 

500L 

OtXlf 


001241 

54 3G 

0001 


CO 12 70 

5610 

DlV.ll 


1X11273 

5660 

0001 


001301 

COOL 

0001 


001367 

6?3G 

oooi 


001436 

6520 

0001 


001467 

6710 

0001 


001 czo 

7256 

coot 


001374 

600L 

0001 


01*1566 

6 1 OL 

cira i 


001 507 

61 5L 

0003 

R 

000001 V 

ADSIR 

0003 

R 

01X1620 

AG5TC 

0003 

R 

00! 320 

ACSCRO 

0000 

R 

1X11X163 

amax 

CKV.V.T 

K 

000062 

AMIN 

0005 

R 

OOGGOl 

AIK.LZ 

1X7.14 

R 

004704 

ATCCNT 

0000 

R 

000126 

BEGIN 

0004 

R 

011610 

CFRESS 

0004 

I 

1*13560 

CEFAUL 

0004 

R 

007640 

CELTRA 

0000 

R 

000057 

CEST 

0000 

R 

0011104 

CEWAVG 

0000 

R 

000077 

CEWC'NT 

OOUO 

R 

0000 70 

DEW T 

0004 

R 

010624 

EFRESS 

IXX14 

R 

012574 

ETEMFS 

00011 

R 

1100074 

EUilLE 

0003 

R 

□ill 010 

EXm >0 

0003 

R 

000144 

tXPCN 

0003 

R 

000632 

EXrCfK 

01X13 

R 

001154 

EXFCfO 

1X7.10 

R 

000 1 2 7 

FINIS 

0000 

R 

000116 

PUNCH 

OtXXI 

R 

CXXX116 

HEADER 

OOGO 

R 

000070 

HEAC1 

DOOM 

R 

01X7.171 

HEAC2 

GOOD 

R 

Ci'XOO?2 

HEAD3 

17004 

R 

002734 

HEIGHT 

coon 

R 

1X10132 

HIGH 

LXV04 

R 

003720 

HLM 

IX.XV.1 

R 

1XXU05 

HUMIC 

17X7.1 

R 

1X70122 

HUMICF 

o.voo 

I 

00LV.1GG 

I 

tXV.V.l 

I 

000107 

IRJNCT 

LXXX1 

I 

1X7.1 1 24 

II 

001X1 


CXX1GS6 

INJFJ 

0000 


0X0652 

iNjrs 

0000 


000675 

INJF3 

c xxv,) 


000662 

INJFJ 

0000 


iXV.i 066 

IN JPi 

OtXXl 

I 

000023 

INSTR 

0003 

K 

000310 

INTER 

0000 

R 

0000170 

INTERV 

0000 

I 

1X7.1 1 06 

J 

0000 

I 

1X701 10 

JFUNCT 

0000 

I 

17X1145 

K 

D 000 

I 

CXXV.167 

LABEL 

0000 

R 

tXXXXU 

leader 

CMOS 

1 

000004 

LEVELS 

0000 

I 

OOOH5 

LIM 

0005 

1 

000000 

LIMIT 

0X7.1 

I 

□00 1 42 

LCD l NT 

otvoo 

I 

1X70144 

LCCSIZ 

OCXVJ 

1 

GOOOG4 

MESSAO 

0000 

I 

01X7021 

NAre 

0017.1 

I 

OHO 1 1 4 

M.AST 

0016 

1 

X.7XXX0 

NXV 

0017 

I 

000000 

NYV 

<7.7.14 

R 

0056 70 

OTCENT 

txvo 4 

R 

014544 

c*:tc 

0000 

R 

ooo: 20 

OZCNEP 

oxv.i 

R 

0170123 

OZCLES 

OtXVO 

R 

000 I 0 1 

FART 

0000 

R 

000133 

PARTI 

oooo 

R 

17.70134 

FART2 

OOOO 

R 

17.701 35 

FART3 

GOOD 

R 

CAXH 36 

FART4 

0000 

R 

1X7.7064 

fh:ne 

0004 

R 

000000 

PRESS 

0000 

R 

000117 

FRESSP 

txxx.1 

R 

17.701 03 

FRE$U< 

0X7.1 

R 

1X7J143 

F* REWAN 

0003 

R 

1X104 54 

KCSPCN 

0005 

L 

0CXV.X13 

SHELL 

0 goo 

R 

1700140 

SLCPEA 

exxx* 

R 

000141 

sloped 

0X01 

R 

002022 

SPHERE 

0000 

R 

0001 30 

SPHERE 

OQUG 

R 

000061 

STEPS? 

C7X.T4 

R 

001750 

TEMPER 

OtXXl 

R 

1X7.1102 

TEfTOT 

017X1 

R 

0170121 

TEHT RE 

017.X."* 

R 

000137 

TEKr 

OGGI 

R 

0021 21 

TERP 

0000 

R 

000131 

TEST 

coos 

R 

1XXX7.12 

THETA 

0X70 

R 

'X.XV065 

T RANGE 

0004 

R 

□06654 

TRANSM 

exxo 

R 

000073 

TWkXE 

0000 

R 

0001 LX# 

VAFER 

LV.X16 

R 

017.7X70 

VAFRES 

0004 

R 

1X70764 

WATER 

LXX7.J 

R 

□00125 

WEIGH 

0003 

R 

COO 644 

HEIGHT 

0000 

R 

000073 

VH1£ 

0000 

R 

1X70112 

X 

O.VXJ 

R 

0.70113 

XP 

0000 

R 

G001 1 1 

r 


OOlOl 1* CCM’lLrRuDATA-SH^T) 

GOltTS 2* EXTERNAL VAPRES 

QUIDS 3* REAL INTEKV 

00106 4* REAL INTER 

COIOT 5 * REAL LEADER 

OOliO 6* LOGICAL SHELL 

Will 7* PARA St TER SlZTRA-lOO 

00112 fl* PARAMETER SIZFTHt5LXi 

00113 9* GO^CN/Zl/ AESCRdSlZTRA) , EXPENDS 1 7TR A) , l NTER <tS I ZTR A) , 

□0113 ID* * RESFCN^SIZTRA) ,A£SCRCalO) , EXFCNCalO), 

00113 11C * If IGHToSIZTkA) .EXFOOiSIZTRA) , EXPCW^SIZTRA) , 

00113 12* * ABSORBS I ZTKA) 

□0114 13* CCJMCH/Z2/ PRESSES ! ZPTH) , WATER*SIZFTH> , TE>E EF*S 1 ZF TH J , 

00114 14* * IE I GHT<i S I ZF TH ) ,HVjr>fci 5 1 ZFTH) , ATCENT*S I ZFTH) , 

00114 15* * OT CENTOS 1 7FTH) f TRANSMfcSlZPTH) .CELTRArfSlZrTH) , 

00114 16 * * EPRESS*SI7FTH) ,CFRESS*SlZFTH) ,£TEMFS*SIZFTH> , 

00114 17* * DEF AUL *i S I ZP TH > pOZafrt SI ZFTH) 

00115 If* CCMMZN/Z3/ LIMIT r ANGLE , THETA, SHELL, LEVELS 

00116 19* 01 HENS ION LEADER o3J 

OOltT PIf*>J?JCN HES^rwfiri) 

00120 21* DIMENSION HEADERS) 

00121 22* ClfCNSlCN NA*-C<t2) ,lNSTR*2tf ) 

00122 23* 0 I MENSICN CEST«2) 

00123 24* DATA SHELL/, TRUE./ 

00123 25* DATA STEFSZ/D.D/ 

00127 26* INTEGER DEPALE 

00130 27* CATA AMIN/-2.0/ 

00132 28* DATA AMAX/11.0/ 

00134 29* DATA NAME/'R CLARK V,CE ST/ "HIP T-D2I3V ,PHCNE/*’Z1 46V 



00140 

30* 


PAT A INSTR/'CNE HARD COPY - CNE MICROFILM - PLOT TO ECFlLE ' T-OZG 

00140 

31* 


*"/ 


CV.U42 

32* 


DATA TRANCE/ 100./ 


00144 

33* 


DAT Adit MrERul ) ,1 = 1 , S I ZP TH > / S 1 7f THv^LVTO . / 


00146 

34* 


CAlArtOZCfCal) f 1 = 1 ,SlZPTH)/SlZrTH*4.67E-6/ 


00146 

35* 

c 



OOl 40 

36* 

c 

COMPOTE RADIANS from ancles 


00 146 

37* 

c 



OOl 50 

38* 


RAOlANC) =ANG/5T, 29578 


OOl 50 

39* 

c 



00150 

40* 

c 

COMPUTE SPECIFIC HlWlDITY FROM PRESS IXvE AND VAPOR 

PRESSURE 

00150 

41* 

c 



cm si 

42* 


RATMlXdV,P>=622.*V/aP-a,378*V) ) 


00151 

430 

c 



00151 

44* 

c 

COLETTE THE I SC EX OF REFRACTION CF WATER-VAPOR AT 

THE TEMPERATURE 

00151 

45* 

c 

AhC FkESSUTC CF A GIVEN LEVEL 


043151 

46* 

c 



00152 

47* 


REFR AXdE 1) - 1 . +d77 . 526E-6) *PRESS^iL 1 ) /T EFFERdL U 


CGI 53 

46* 


CALL TAFLOTdlHM) 


00154 

49* 


CALL ICFRMVaN/.rt: ,CEST , FHZKi , 1 NSTR) 


00135 

50* 

ICO 

CC7H iNUE 


00155 

51* 

c 



COISS 

52* 

c 

INPUT CONTROL CARDS 


00155 

53* 

c 



00156 

54* 


REACdS'll ,EKD=500) LI MI T, LABEL , HEAD1 ,1GAC2,K€AD3 


00156 

55* 

c 



OOl 56 

56* 

c 

50 FLAGS TO AC JUST FRECIPI TATE-LE WATER COMMUTATIONS FOR 

00156 

57* 

c 

REFRACTION AT HIGH ALTITUDES 


OCU56 

56* 

c 



00165 

59* 


1F*LA£EL*ES. ’SPHERE-) SrFLL- * TRUE . 


00167 

GO* 


IFdLABEL.E«3.-FLAT 5PCLL- .FALSE * 


00167 

61* 

c 



001G7 

62* 

c 

ADJUST DISPLAY WINDOW CN T>d SC-4020 PLOTS 


00167 

63* 

c 



00171 

64* 


IF*LAEEL , EO . “ MAXCCR " ) AMAX=L IHI T 


00173 

65* 


IFdLADEL.EO.-MlNCCR*) AMIN=LIMIT 


00173 

66* 

c 



00173 

67* 

c 

input the angle cf observation 


00173 

66* 

c 



00175 

69* 


IF*LAEEL *EQ. “ANGLE “) ANGLE=L IHI T/10, 


00175 

713* 

c 



00175 

71* 

c 

INPUT RESPONSE FUNCTIONS 


00175 

72* 

c 



OOl 77 

73* 


I F*LA BEL. EQ." RE SFCN") GO TO 200 


00177 

74* 

c 



00177 

75* 

c 

INPUT RADIOSONDE DATA 


00177 

76* 

c 



00301 

77* 


IF^rEL . EO . p WATER *! GO TO 300 


00201 

76* 

c 



00201 

79* 

c 

IH=UT TEMPERATURE RANGE CF 100 POINT TABLE 


00201 

ea* 

c 



00203 

81* 


I F*LABCL.EQ.“ TEMPER*) GO TO 400 


00203 

82* 

c 



00203 

63* 

c 

EXECUTE WITH CATA FROM PREVIOUS SOUNDING 


00203 

84* 

c 



00205 

85* 


I RELABEL. EQ.*EXECUT*) GO TO 390 


00205 

66* 

c 
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OOZtiS 

67* 

C 

Recede a fressire Increment to establish t>c sicf size fcr radiosonde 

00205 

66 * 

c 

DATA FROM THE SURFACE TO ONE MlLlICAR 

00203 

69* 

c 


CIOZOT 

90* 


I F«L ADEL. Ed. “STEP -> StEFSZ=LlMlT 

cum? 

91* 

c 


00207 

92* 

c 

INPUT KISSACE LlhCS DESCRIBING OUTPUT 

00 WT 

93* 

c 


ooai 1 

94* 


lF«tLACO_ .EO.-WZSSAQ-) CO TO 600 

cozi i 

95* 

c 


0021 1 

96* 

c 

CNER-LRITE MISSIS RACK^IfCE DATA TO SATELLITE ALTITUDE WITH DEFAULT 

002 i i 

97* 

c 

DATA 

ooet i 

96* 

c 


EU213 

99* 


iFaLADEL.Ed.-CCFAUL*) CO TO 800 

O02I3 

100* 

c 


00213 

IOI* 

c 

IrCUT SPECIAL TRACE CAS TRAhiSMl SSlCN FUNCTIONS 

00213 

102* 

c 


0021 5 

103* 


IFoLACEL.EQ.-TRACE-) GO TO 1 100 

00217 

104* 


CO TO 100 

00220 

105* 

300 

CCNT1MJE 

00220 

106* 

c 


00220 

107* 

e 

RECORD M>CCR CF LEVELS FCu IhTUT RADICOCNCE 

00220 

1 G 6 * 

c 


00221 

109* 


LEVEL S^LlHl T 

00221 

HO* 

c 


002 21 

111* 

c 

SET OBSERVATION ANGLE 

00221 

112* 

c 


00222 

113* 


THTTA=RADfxANGLE) 

00222 

114* 

c 


00222 

115* 

c 

INITIALIZE FRECIFI TABLE WATER COMMUTATION 

00222 

13 6 * 

c 


00223 

117* 


WHOLE -0. 

00223 

116* 

c 


00223 

119* 

c 

initialize effective PRESSURE CCPFUTAT ICN 

03223 

120 * 

c 


00224 

121* 


EVJHGLE-0 

00224 

122 * 

c 


00224 

123* 

c 

initialize effective temperature opi/tation 

00224 

124* 

c 


00225 

125* 


TVMJLE=0. 

00225 

126* 

c 


00225 

127* 

c 

IN~-LTT SURFACE RACIOSC^CE DATA 

00225 

126* 

c 


00226 

129* 


REACTS ,12J FRESSrtl) , CEWF T , TEMPER* 1 ) , HE I GHTrtl ) .OtZCMInl ) 

00235 

130* 


CEFAU_<tl)=0 

00236 

131* 


I F*HEAC 2 .NE. -KELVIN-) DEWFT=DEVTT*273 . 1 6 

OG24Q 

132* 


I F*HEAC2 . NE „ -KELVIN" ) TEM“ERal > -TEKPEfial ) +273 . 1 6 

0024"^ 

133* 


IF w* d ADD . RE » "CCVT FI " ) FLC" £ , 1 ; A'.Lal ) ) " 1 

00244 

134* 


1 FaHF a03 . NF . •n7nrjr-> r»7i-toc* n-/< e?tr_4 

00246 

135* 


DEWFNT=CEW=T 

00246 

136* 

c 


0024 6 

137* 

c 

COTUTC SURFACE VAPOR PRESSURE 

0024 6 

138* 

c 


00247 

139* 


VAFER - VA F ft E S*C EWF T ) 

0024 7 

140* 

c 


0024 7 

141* 

c 

USE DEW-POINT GREATER THAN 77 DECREES TO ASSUME 10 PERCENT 

00247 

142* 

c 

RELATIVE HUMIC ITT 

0024 7 

143* 

c 




00250 

144* 


iFaCEVJH . GE .350.1 VAFER^VAPRES'tTEMFERo.l ) > *. 1 

W252 

145* 


|FtfDtV*T.GC.350.> FLO.3 , 1 , CEF AUL«* 1 1 > = l 

00252 

146* 

C 


00252 

147* 

c 

COMPUTE THE SURFACE SPECIFIC HUMIDITY 

00252 

140* 

c 


Cy.'i 254 

149* 


FART=RATMl XitVAPCR ■ PRESS'* 1 1 ) 

00254 

150* 

c 


002 54 

151* 

c 

STORE THE SPECIFIC HUMIDITY CONFUTED AT THE SURFACE 

OJ254 

15?* 

c 


00255 

153* 


M>*il)=FART 

00255 

154* 

c 


00255 

155* 

c 

T>C ATMTS-CENTlFtlERS CF CARDCN DIOXIDE ANC OZCt£ AT THE SURFACE 

00255 

15G* 

c 


00256 

157* 


ATCENT«1>=0. 

00257 

158* 


OTCENTal)=0. 

00257 

159* 

c 


00257 

160* 

c 

THE EFFECTIVE FRESSURE AT T Ft SURFACE IS EOUAL TO UC RADIC6CNCE 

OLJZ57 

161* 

c 

SURFACE FRESSURE 

00257 

162* 

c 


002 GO 

163* 


EFRESS^l ) ^FRESSal) 

00261 

164* 


CFRESSrtl ) -FkESSh1 ) 

00261 

165* 

c 


C732G1 

186* 

c 

LABEL THIS FACE CF ACSZRDINO CAS Gl>KENTRAT ICNS. 

00261 

167* 

c 


00262 

168* 


VR !TE«a6 , 27) «OtSSAC<xI ) * I "1 , lO) 

002 70 

169* 


CEVF T -CEW'T-27 3 . 1 6 

00271 

mi* 


TOFCR-TEfFERril 1 Z73.16 

00272 

171* 


UFITE*6,14) FRESS^il) pCEVFT.TEMKR 

00277 

172* 


DO 10 I^ZtLlMlT 

00277 

173* 

c 


00277 

174* 

c 

EXTRACT A LEVEL -S DATA 

00277 

175* 

c 


00302 

176* 


REAC*i 5 ,12) PRESSED , DCVFT , TE>FER«I ) ♦ FE I OKT* I ) , OZONE* I) 

00311 

177* 


DEFAULT I ) ”0 

□0311 

178* 

c 


□D311 

179* 

c 

CONVERT 1CICHT TO FEET IF NEEDED 

0031 1 

160* 

c 


00312 

161* 


iF-FeAOl . ES . “ KETERS" 1 hC I GHT*1 ) =HCl GHTd I > *3 . 2800 

00312 

1624 

c 


00312 

163* 

c 

CCFFUTE AVERAGE FRESSURE CF A LAYER 

0031 2 

184* 

c 


00314 

165* 


FRESUR^aFRESSdl) >rRESSaI'l) ) /2. 

00314 

166* 

c 


00314 

167* 

c 

CONVERT FRjCH CENTIGRADE TO KELVIN IF NEEDED 

00314 

186* 

c 


00315 

169* 


IFaHEAC2 .NE * "KELVIN" ) CELFT =CCVFT+273 . 1 6 

00317 

190* 


1F«*>CAC2 .NE . “KELVIN") TEfF ER* 1 ) =TEH~-EfttfI > +273 . 16 

□0317 

191* 

c 


00317 

1 92* 

c 

It* \x*_ A UUM 5 IAN 1 CVCNt FRCMLL IK NEEDED 

00317 

193* 

c 


00321 

194* 


IF«HEAD3.W:. -OZONE "> OZONE*! ) =4 .67E-S 

00323 

195* 


IF«*HEAD3.KE. -O20NC ") FLCh2, 1 , DEFAULT I) 1 = 1 

00323 

1 90* 

c 


00323 

197 * 

c 

COMPUTE THE AVERAGE LAYER DEW POINT 

00323 

196 * 

c 


00325 

199 * 


CFWA VG-dCEV-RNT ♦ CEVJF T > /2 . 

00325 

200* 

c 




1X1325 

201 * 

C 

COMPUTE THE YAPCft PRESSURE 

00325 

202* 

c 


00326 

203* 


VAPEfUVAFRESuCEWAVG) 

00326 

204* 

c 


00 32G 

205* 

c 

USE CEW-rolNT GREATER THAN 77 DEGREES TO ASSUME 10 FERCENT 

00326 

206* 

c 

RELATIVE HUMICITY 

00326 

207* 

c 


00327 

206* 


lFttCEW-T.CE.3S0.) VAFER-VAPRES*TEMFER«il > ) *. 1 

00331 

209* 


IF*CCWT T »GE*3S0.) FL D*3, 1 ^CEFAULal ) ) =1 

00331 

2 10* 

c 


00331 

211* 

c 

COMMUTE THE MIXING RAUOCF THIS LAYER 

00331 

212* 

c 


00333 

213* 


H-MI D-622 . *V AF E R /ctp RE S UR - V AF ER ) 

00333 

214* 

c 


00333 

215* 

c 

COTUTE T>-e AVERAGE TEMPERATURE CF THIS LAYER 

00333 

216* 

c 


00334 

217* 


TEWKR^* TEMPER* I ) ♦TEJ^FER*! -1 )) /2 . 

00334 

218* 

c 


0-1334 

219* 

c 

COrUTE T PEIGHT CF THE TCP CF THIS LAYER IF fCEDEO 

□0334 

220* 

c 


00335 

221* 


lF*>fE!GHT*I) .LE.O.I 

0.1335 

222* 


WEIGHT** I > =HElGHTciI-l) +29.3*TEM"-OR *ALCC«FRES&xI'l ) /PRESSED > 

01335 

223* 


*** 1 . ♦ .00061 a 1LM I D) *3 . 2808 

00335 

224* 

c 


00335 

225* 

c 

STCFE TVC SPECIFIC HUMICITY AT EACH LAYER TCP 

00335 

226* 

c 


00337 

227* 


KM4l)=RATHlXaVAFER,FRESU<) 

00337 

228* 

c 


00337 

229* 

c 

COf^LTTE TVe FRECIFITAClE water at this level 

00337 

230* 

c 


00340 

231* 


WATER* I ) =WATFR*I-1 ) ♦*HUM*I > *SFHEREa|10 , 1 J **PRES5* I '1 ) -PRESS*] )))/ 

00340 

232* 


* 980, 

00340 

233* 

c 


00340 

234* 

c 

CCFFLTTE TTC ATMDS-CENT I UTTERS CF CARECN-ClCfcl DC AlsC GZChE TO THIS 

0.1340 

235* 

c 

LAYER 

00340 

236* 

c 


00341 

237* 


AT CENT a I ) =aFRES5* 1 ) -PRESS* I ) ) * . 260 

00342 

238* 


OTCENTal > =OTCEWT* I -1 ) ♦ap RESS* I - 1 > -PRESS* 1 ) ) *dOZCtsE* 1 ) ♦CtfCNE* I - 1 ) ) * 

00342 

239* 


* 237.968 

00342 

240* 

c 


00342 

241* 

c 

COMPUTE THE EFFECTIVE TEM' ERA TUFFS ARC PRESSURES USED TO RELATE 

00342 

242* 

c 

HOOGEhCOUS PATH TRAN SMI SSI CM FUNCTIONS TO Hi tTHXFHFOUS SLANT 

00342 

243* 

c 

PATHS THROUGH THE ATMOSPHERE 

DD34 2 

244* 

c 


00343 

245* 


ewhole~e\j^:le+*water*i ) -water* i-i > )*presur 

00344 

24 6* 


TVM^E-TUiXE+*WATER* I ) -W-\TER*I -l ) 3 *TEMKR 

00345 

247* 


CFRESS*! > =otFKE5S*l > ♦PRESS* 1 ) ) /2. 

IXJ34 6 

246* 


tf-Kt *s*i ) =Eva^:lefwater» i > 


n *o-*- 


r^C^T \ r/WATTS^n 

0034 7 

2Vl* 

c 


00347 

251 * 

c 

COMPUTE CENTIGRADE VALUES FOR DISPLAY 

00347 

252* 

c 


00350 

253* 


CEWAVC^DEVF T-273 . 1 6 

00351 

254* 


TEKFORr TEMPER* 1 ) -273 . 3 6 

00351 

255* 

c 


00351 

256* 

c 

PRINT INTERMEDIATE RESULTS CF THIS PROCESS 

00351 

257* 

c 




00352 

2364 


U UTE*6,13> PRESS*!) , CCWAVG , TEMPCR , HE I GHT* 1 1 .WATER*!) , ATCENT*! > , 

00352 

2534 


4 OTCENT* 1 ) iEFRCSS*!) ,CFRESS*I) .ETEMPSa ( ) 

60306 

2604 


CEWPNT=DCWPT 

00307 

2614 

10 

CON! I NUE 

00371 

2624 


OO TO 100 

00372 

2634 

11 

FORMAT* I 10 , A6 , 6X , 3 AG) 

00373 

264 4 

12 

FCRMAT*1 X,F7 *1 ,2X,F6.1 ,4X,F7. 1 .3X.F10.1 ,ElO. 5) 

00374 

2654 

13 

FCRMAT*IX,F7.t ,2x,P6.1 ,4X,F7«1 ,3X ,F10.1 , 1 X ,F 7 . 4 , 1 X ,F6.1 ,F6.3 ,2X, 

00374 

2664 


4F6.lt 

00374 

2674 


46X.F6.1 ,6X,F5.1) 

00373 

2664 

14 

FORMAT*" PRESSURE DEW POINT TEMPERATURE HEIGHT FREClP ATMD6-CE 

00373 

2694 


4NT EFFECTIVE CONSTANT EFFECT I VE" ,/,• • ,F7 . 1 , 2X ,F6. 1 *4X .F7 . 1 

00375 

2704 


FEET WATER 00*2 OZCNE PRESSURE CONCENTRATION TEMP- 

00375 

2714 


ER A TORE- ) 

00376 

2724 

15 

FORMA7*F7.1 ,F7.4) 

00377 

2734 

16 

FCRMAT«1X,F7.1 ,F7.4> 

00400 

2744 

17 

FCRMAT*7H ANGLE ,F4.1 ,20H DEGREES FROM NADIR ,/,20H SURFACE TESTER 

00400 

2754 


4ATURE,F7.1) 

00401 

2764 

18 

FORMAT*! X,36HTEW-ERATURE RANGE OF 100 POINT TAELE) 

00402 

2774 

19 

FORMAT* 42H1WAVENUCDR MIDPOINTS APC PEkCOH RESPONSE) 

00403 

2764 

21 

FORMAT* 37H1 TEWERATURE RANGE OF UV.l POINT TAELE) 

00404 

2794 

22 

FORMAT*' PRESSURE PR EC IF HEIGHT SPECIFIC TEW' CORRECT I 

00404 

2604 


4CN-,/, " WATER HUMIDITY SILNClNG" 

00404 

2614 


4 -1© CM FT GM/KGM K K") 

00405 

2624 

23 

FCRMAT*1H ,3A6) 

03406 

2634 

24 

FCRMAT*F4.1 ,14H DEG. TO NADIR) 

00407 

2644 

25 

FCRMAT«F5.i ,13H SURFACE TE>F) 

00410 

2654 

26 

FORMAT* 10A6) 

00411 

2664 

27 

FORMAT* 1 HI ,lOAG,/> 

00412 

2674 

26 

FCRMAT*4E10.5> 

00413 

£664 

29 

FCRMAT*1 X , F6 . 1 , - , K6 . 1 , 4*3 X . F7 . 4 ) ) 

00414 

2894 

31 

FCRMATnl HL1 ,3A6, /| * WAVENUCER A4aATMJ6-CENT4=iE) 4a£F^*I) * ./, 

00414 

2904 


4- INTERVAL ABC VEIGHT*) 

00415 

2914 

33 

FORMAT*! X,F7. 1 ,2x,ElG.5 ,2X.F6.1 ,2X,EI0.5) 

03416 

2924 

35 

FORMAT* * 1 DEFAULTED STRATOSPHERE PROFILES- 

03416 

2934 


4 “ PRESSURE OZCNE TEW SPECIFIC-,/* 

□0416 

2944 


4 • SOURCING HUMIDITY*,/, 

00416 

2954 


4 "1© GM/KCH K GM/KGM- ) 

00417 

2964 

390 

CONTINUE 

00420 

2974 


IFaSTEFSZ.GT.G.) CALL STEFER 

00422 

2984 


IF*STEPSZ.GT.O.) CALL CISPLY 

00422 

2994 

C 


00422 

3004 

c 

ESTABLISH A TEWERATURE DlACKDUCY INTENSITY SCALE 

00422 

3014 

c 


00424 

3024 


CALL TABTEM*TEWER*1) .TRANCE) 

00424 

3034 

c 


00424 

3044 

c 

PRODUCE A DISPLAY TO SPCW TEMPERATURE CORRECTION AS A FUNCTION CF 

00424 

30^4 

c 

L!fC CF SICHT ANCLE,. HE KMT TM fETEFS, AfC FT;EC I P I TASLE WATER. 

00424 

3064 

c 


00425 

3074 


ViR I TE *6 1 27) oMESSAG*I) ,1 = 1 ,10) 

00433 

3084 


IF«SHELL) LEAQER«3) =- SPHERE" 

00435 

3094 


WRITE*6,23) LEADER*1> ,LEADER*2> ,LEACER*3> 

00442 

3104 


WRITE*6,17) ANGLE, TEMP ER*1) 

00446 

3114 


Vf? I TE*6 , 22 > 

00450 

3124 


CALL SETMlVaJCO, 100,100,50) 

00451 

3134 


CALL CR 101 V*1 , AMIN, AMAX , 0. ,60000. , l . ,70000. ,0,0,1,1,3,31 

00452 

3144 


CALL PR INTV* 16 , LEADER , 703 ,46) 


ICW 


icv 

*£W 

rew 

NEW 
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FO£,S RAC ►KC/SL ANTI , RACKC , RACMIC/Sl ANT 1 

FOR 94L-n«/2L'r-l«'Stt 


SLBRCUT I hC RACMIC ENTRY POINT LXV.TZltt 


STCftAGL USCC ' CCCC<tl ) 000222 CATA^O) »V.XILT<6 ELAfC; COMt-Cfto?) IXVOBOO 
COtKN BLOCKS 1 


0003 

21 

CM 14 64 

W04 

22 

002570 

0005 

23 

(XXX7.15 


EXTERNAL REFERENCES -BLOCK, NA^C) 

0UG6 TRANS 
0QO7 TEMTAB 
DO ID COG 

0011 EXP 

0012 NERR31 


STCRAOC 

ASSIGN* -ENT -BLOCK , TYPE, RELATIVE LOCATION, NAME) 










0001 


000005 

1216 

POOL 


0X7.7020 

1276 

root 


000126 

1436 

0001 


000157 

30L 

0003 


000000 

ABSCft 

0003 


000620 

AE5CRC 

0003 


00 J 520 

absiro 

0005 


1X7.7.7.11 

ancle 

0004 

ft 

000764 

AT CENT 

0000 

ft 

01X701 3 

atmds 

0004 


601 7aO 

CftrtLSS 

■VV'uJ 



CfcT ALL 

l>.Xi4 

* 

001440 CELTRA 

ocoo 

ft 

0000 1 4 

OETECT 

GOTO 

ft 

000015 

CftTEMr 

COC4 


001604 

STRESS 

COE*4 


002114 

ETE>FS 

1X703 


00 10 10 

EXFO-O 

D003 


rv~j 

E XftCN 

0003 


OLA7G32 

EXPCNC 

0003 


C01 1 54 

EXPOND 

CX7.T4 


000454 

HEIGHT 

0004 


000620 

HUM 

0000 

1 

0001X14 

I 

0000 


000024 

!njp$ 

0000 

R 

COLVOOl 

INTEN 

CXX13 

R 

000310 

INTER 

0000 

R 

LXV.vvvi 

INTFRV 

0000 

I 

CM0D02 

J 

0005 


17.7.7.7.14 

LEVELS 

0005 

1 

cxvjtxvr 

LIMIT 

0004 

R 

001130 

OTCENT 

0004 


002424 

0004 


0LX7JX7J 

PRESS 

0003 

ft 

0004 54 

RESTVN 

OOOO 

R 

EXVTOl 1 

sFcce 

0005 

L 

1X70003 

SHELL 

0004 

R 

1X703 10 

TEMPER 

0000 

R 

OLXH312 

TEhFOl 

gods 

ft 

0X7.7X12 

THETA 

0006 

R 

DtXV.XVO 

TRAMS 

0004 

R 

001274 

TRANSM 

tXX '14 

R 

IXX1144 

WATER 

0000 

R 

0004X16 

WATERC 

0000 

ft 

0t7.X7.i7 

WATEkO 

CXX’jO 

R 

000005 

waters 

GOTO 

ft 

CXXXXB 

UAVEKO 

1X703 


000644 

VE1GHT 

0000 

R 

OOOOUl 

VH3_E 






DO101 

1* 


SUBROUTINE 

RACMXJ— DLTE>F ) 


ooioi 

2* 

C 




□aim 

34 

c 

THIS SLBRCUTIMI INTEGRATES THE RADIATIVE TRANSFER FUNCTION USING 

aoi oi 

4* 

c 

THE TRAFEZOIC RULE. 


□OlOl 

5* 

c 




D0103 

6* 


LOGICAL shell 


00104 

7* 


F^TFRNAJ TRANS 


00105 

8* 


REAL Iff! Eft 



00106 

9* 


REAL I NT Eft V 

, INTEN 


□0107 

to* 


PARAMETER Sl2TRA=lGO 


oano 

11* 


PARAMETER SI ZFTH-10G 


00111 

12* 


C0MMDN/Z1/ 

AESOR-SI ZTRA) , E7FGN-SIZ7RA) , 

INTER-SI Z7RA) , 

oain 

13* 


* 

RESPONDS I ZTRA) ,AB&}“C«lO) , 

EXFCNC-10) , 

OOlll 

14* 


* 

WTlGHTaSlZTRA) .EXFCHDttSl ZTRA) 

, EXPONAS I ZTRA J , 

OOJ 1 1 

15* 


* 

ABSORCHSIZTRA) 


□0112 

16* 


COHHCN/Z2/ 

PRESS— SI ZPTH) , WATEft-SIZPTH) r 

TEMFER-S JZrTH) , 



00112 

174 


4 HClCMT*SlZrTH> ,HUM'i5l ZPTH) , ATCENTaSl ZFTHJ . 

cm 1 ? 

164 


4 OTCENTaSlZPTH) , TRANSMnStZFlM) .CELTRAaSl ZPTH) , 

00112 

194 


* EFRESSwSlZRTH) .CPRESIVtSlZFTm ,ETEMPS«SlZPTH) , 

00112 

204 


4 CEFAULrtSlZrlHl t OZCNEa$lZfTH) 

00113 

214 


COU4^/Z3/ LIMIT, ANGLE, THETA, 5HCLL, LEVELS 

00113 

224 

C 


om 1 3 

234 

c 

ARITHKCTIC FUNCTIONS 

001 1 3 

244 

c 


0.1114 

254 


B&XYaV.T) =<t6.9349E-H*V ^) /aCXF«ici| ,43654V) /T)-1.0) 

00115 

264 


TEMPED, V,W) =«l ,43654V) /ALCCa«'aW^ .93460 1 34V443) /C) +1 .) 

00115 

274 

c 


cm 1 5 

264 

c 

INTEN IS THE SUMATlCN C F PRODUCTS clRACI ANT INTENSITY lifts 

cm is 

294 

c 

WAVE NUCER INTERVAL) 

00115 

304 

c 


00116 

314 


INTEND. 

00116 

324 

c 


00116 

334 

c 

InTERV is the wave nuhcer interval 

0.1116 

344 

c 


00117 

354 


INTERV-25, 

00120 

364 


DO 30 J=1 ,65 

00120 

374 

c 


00120 

364 

c 

SMF INTERVALS HAVING A ZERO RESPONSE FLNCTICN. 

00120 

394 

c 


00123 

404 


1 F^RESFCNa J ) . LE ,0. ) CO TO 30 

00123 

414 

c 


00123 

424 

c 

WAVENO IS THE KEAN VAVENU-CCR CF THE INTERVAL BEING EVALUATED 

00123 

434 

c 


00125 

444 


WAVE>0=lNTEK«J) 

01126 

454 


CO 10 1=1 ►LIMIT 

00126 

464 

c 


00126 

474 

c 

COMPUTE Trt CFTlCAL PATHS CF ABSORBING CASES 

C£U 26 

464 

c 


00131 

494 


wATEKS=<iWATER4LlMlT> -WATER^I ) ) /COG* THETA) 

00131 

504 

c 


00131 

514 

c 

FOR SATELLITE ALTITUDE COFUTAT 1CU5 THE OPTICAL PATH CF WATER 

00131 

524 

c 

VAPOR MJST INCLUDE THE INDEX CF REFRACTION CHANGES 

00131 

534 

c 


00132 

544 


IFaSHELL) WATERS - WaTERoL IMI T) -WATER* I ) 

00134 

554 


WATERC=*ATCENT<*LlMtT) -ATCENT * I) } /COSa THETA) 

00135 

564 


WATERO=aOTCENT<iL IMIT) -OTCENTcxI ) ) /COS-iTHETA) 

00135 

574 

c 


00135 

564 

c 

COFUTE THE TRANSMISSIVITY FCR EACH LEVEL 

□0135 

594 

c 


00136 

604 

10 

TRANSHt I ) = TRANS*WAVEND , WATERS , WATERC , WATERO, I ) 

00136 

614 

c 


00136 

624 

c 

initialize the upuellinc radiation 

00136 

634 

c 


00140 

6it 


V.HtT-E^C' . 

00140 

654 

c 


00140 

664 

c 

FOR this INTERVAL COMPUTE THE EXTINCTION OF A SURFACE IMAGE. 

00140 

674 

c 


00141 

664 


SFC B B - TR AN SKrf 1 ) ftCBOCYrtWAVEND, TEKFERr*! > ) 

00141 

694 

c 


00141 

704 

c 

USING WATER VAPOR, CARBON DIOXIDE AND OZONE AS THE IMPENDENT 

00141 

714 

c 

VARIABLES INTEGRATE THE CCNTRICUT ICNS CF EACH ATMOSPHERIC LAYER 

00141 

724 

c 

TO THE IMAGE IN THIS INTERVAL 

00141 

734 

c 
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VOX 42 

74* 


03 20 1 = 2 T 

ooi 4? 

75* 

c 


IVjHZ 

76* 

c 

FIND THE CHANCE IN TRANSMISSION BETVCEN LAYERS 

Wl 42 

77* 

c 


001 4 s 

76* 


CELTRA«x|> -TRANSMal 1 ~TRANSMil~l ) 

00145 

79 V 

c 


0014 5 

60* 

c 

ribC AVERAGE TEMPERATURE CF A LAYER 

00145 

61* 

c 


O'-*! 46 

62* 


TEHFOL^a TEMPER* 1 ) 4 TEMPER*! * 1)) /Z . 

1X1146 

63* 

c 


00146 

64* 

c 

CCJKFUTL THE EMISSION FRCH THAT LAYER 

00146 

65* 

c 


00147 

66* 


ATHDS^CGCCTuWAVEND, TEMrCRl 

0014 7 

87* 

c 


00147 

66* 

c 

VHXE IS THE SUM CF fROCUCTS *ATMj6FHER IC INTENSITY TIMES CHANGE 

ail 4 7 

69* 

c 

IN TRANSMISSION) 

an 47 

90* 

c 


GOl so 

91* 


vHlLE=VJ^lLE+aATHC^CO_TRA<i I ) ) 

an 51 

92* 

20 

COMTINU C 

ansi 

93* 

c 


aii si 

94* 

c 

INTEGRATE LFVELLING RADIATION FCfi THIS INSTRUCNT 

00151 

95* 

c 


an 53 

96* 


l KTtN= I NT EW+ **W CL E * SFC CC ) <RC SFCNi J ) 

an 54 

97* 

30 

CONTINUE 

an 56 

93* 


DETECT = INTET** 1 NTERV 

00136 

99* 

c 


00150 

ICO* 

c 

EXTRACT A HATCHING DETECTOR TEMPERATURE FROM THE CALIBRATION TABLE 

01156 

1Q1* 

c 


00157 

102* 


CALL TEHTACctCCTECT t CRTE>F) 

00160 

1CI3* 


DLTEM^TEM LRal > -ORTEHf 5 

an go 

104* 

c 


GO 160 

105* 

c 

DISPLAY TEKT-ERATURE CORRECTION AS a FUNCTION CF FRECl Pi TABLE WATER 

an ten 

106* 

c 

HEIGHT CF THE OBSERVATION PLATFORM. 

ooieo 

107* 

c 


DO 161 

106* 


RETLKN 

00162 

109* 


EFsC 


E>C CF 

COMPILATION* NO DIAGNOSTICS. 
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roR ,s RACM1C/SLANT ,RACh:C.RACM:C/$LANT 
rOR 9«L- l ."*tf/2lt-10 > S r .j 


SUGROUT 1 NC RADMIC ENTRY FOlHT W0276 


STORAGE USCC * CCCE*1 > '-XX13H6 CATAa'.l) IX'OGGS BLAW CCNM;N«2> OT/ti 

cchm:n dlccks* 

eve'll 21 001464 

0004 22 015530 

CAX15 23 DLXXXi5 


EXTERNAL REFERENCES «CLCCK , NAKI> 


0036 

TRANS 

0007 

TEMTAB 

001 Q 

L 1NEV 

coil 

NXV 

0012 

NYV 

0013 

CC6 

aoi4 

EXP 

OOj 5 

NWCUS 

0016 

Nice$ 

0017 

NERRit 


rtELOCX, TYPE i RELATIVE LOCATION, NA^t) 


STORAGE ASSIGNMENT 
OGOO 000023 nr 

twos ocjoooci aes:r 
OGOO r noa.us AT Mis 
DGLVD R 000020 DLTEYF 

0003 000144 CXFCN 
GCXXl I LXVJGG6 I 
COCO 1 01X7004 J 

0004 R 005670 OTCENT 
OLVJ5 L CtXXTOi SJCLL 
DOTH R GG6G54 TRANSM 
OOGO R 000005 WAVEfO 
0000 R GOLTOZZ Y 


0001 GG0021 1250 

0003 000620 AESIkC 

0004 011610 CFRESS 

0000 R IXXX317 DRTEMP 

0003 000632 EXrCNC 

00*0 000042 1NJFS 

0X1 5 000004 LEVEL 5 

0004 014 544 OZCt-E. 

0004 R 001750 TEh^CR 
DlXU R 000764 WATER 
00113 000644 WEIGHT 

0000 R QOOJG3 Yr 


00.11 000034 133G 

0003 0O132 , 0 AfiSCwO 

0004 0 135a! ce F A UL 

0004 010624 EFR ESS 

0003 001154 EXFttO 

0000 R Gi.AV.101 1N7EN 

0005 I LX7XXX1 LIMIT 

0004 R 01VXV.YJ PRESS 

oooa r rooou temfor 

0000 R 1V.VJ010 WATERC 
0000 R 000012 VHLE 


axil DOLT 142 147G 
OTA! 5 OXXXll Angle 
0004 R 0.17640 CELTRA 
0004 012574 ETEW-S 

0004 R 002734 HEIGHT 
0X33 R 000310 INTER 
0011 I LXXXXX1 NXV 
0003 R 0X1454 RESFCN 
LXX15 R LXXXX12 THETA 
CXV.X3 R 00001 1 WATERO 
0300 R OXX321 X 


0301 CXXj 1 73 30L 
0004 R 0.14704 ATCENT 
0000 R LXXX.16 DETECT 
0003 001010 EXrCtO 

0304 R 003720 HUM 
0000 R 0LX1000 INTEkV 
00*12 { LXXXXX1 NYV 
00013 R iXV.Xi 1 3 Sr CEB 
0006 R LV.XXXXi TRANS 
0X10 R 000X17 WATERS 
0000 R 0X1002 XT 


nrnoi 

i *■ 


SlRPrirrTKT RAOMrT 

CGI01 

2* 

c 

THIS SlCRCUTlNE INTEGRATES THE RADIATIVE TRANSFER FLJJCTICN USING 

oaioi 

34 

c 

oaioi 

44 

c 

THE TRAFE20ID RULE. 

00101 

5* 

c 


00103 

6* 


LOGICAL SHELL 

00104 

7* 


EXTERNAL TRANS 

00105 

e* 


REAL INTER 

00106 

do 


REAL lNTERV , INTEN 
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Gilt 07 tli* 

Willi n* 

Will 12* 

W III 13* 

Will 14C 

COl 11 15* 

00112 16* 

0,1112 17C 

OO 112 16* 

00112 19* 

00112 20* 

00113 21* 

00113 22* 

00113 23* 

00113 24* 

00114 25* 

001 1 5 26* 

00116 27* 

00120 26* 

00120 29* 

00120 30* 

00120 31* 

GO 1 20 32* 

00122 33* 

00122 34* 

□0122 35* 

00122 36* 

00123 37* 

00124 36* 

00124 39* 

OOl 24 40* 

00124 41* 

00127 42* 

00127 43* 

00127 44* 

00127 45* 

00131 46* 

00132 47* 

00132 46* 

00132 49* 

00132 50* 

00135 51* 

00135 52* 

00135 53* 

00135 54* 

00135 55* 

00136 56* 

CC243 37* 

00141 56* 

00141 . 59* 

OOl 4 1 60* 

00141 61* 

00142 62* 

00142 63* 

□0142 64* 

00142 65* 

00144 66* 


PARAMETER SUFTHsVAO 
PARAMETER SJ ZTRAr 1 i.Vj 

COMHtN/ZI/ ABSCKuSlZTRA) * EXFCNuSlZTRA) , INTERNS l 2TRA) , 

* R£SP<N4SIZTkA) .ABSCROtlO) , exFC«C*Ui) , 

* VC 1 CHI u 5 1 Z Tri A } .ExPCHTuSIZlriA) , C XrCMDu S I Z TR A) , 

* AESO~CwSlZ7RA> 

CCHHX1/2Z/ FRCS*i$lZnHl , WATEfWS I ZPTM) , !£!'," CRaS 1 7PTH> , 

* HElCHlciSlTFIM) pHUViSlZPTH) p ATCENT*SIZPTM) , 

* OTCENTaSlZFlH) , TRAN$Ma$I Zr THj ,CCLTRA«M ZPTH) , 

* EFRESS4S1 ZrTH) pCPRESS'f Si ZP TH) ,ETEMP$*iSI ZPTH) , 

* CETAUL*SI ZP1H) ,OZa,E4 SI ZPTH) 

ccwm:n/23/ limit, angle » theta, shell, levels 

c 

C ARlTHhCTlC FUNCTIONS 

c 

EOaCYaV.T)-ae.9349e-13*V**3)/*EXr*al . 4365*V) /T) -1 .0) 

TEVTa-C » V , W) — u 1 .436S*V)/ALCC4<triW^.S34GE-13*V**3)/C)*l.> 
|FctLlMlT.E<3.2> XF=0.a 
lF4tLlMlT.E-2.2J YP=0.0 
C 

c INI EH IS THE SUT-MATICN CF PRODUCTS ctRACIANT INTENSITY TlKCS 
C WAVE NLFCER INTERVAL) 

C 

INTEND. 

C 

c INTERV IS TIE WAVE NUMBER INTERVAL 

C 

INTER V=25. 

DO 30 J*I ,65 
C 

c SKIP INI LwVALS HAVING A ZERO RESPONSE FUNCTION. 

C 

lF<*R£SFCNtf J) .LE.O.) GO TO 3G 
C 

C WAVEND IS THE KEAN WAVENUMBER CF THE INTERVAL BEING EVALUATED 

C 

WAVENO= IhJTERa J ) 

DO 10 1-1 , LIMIT 
C 

C COMPUTE THE CTTlCAL PATHS CF ABSORBING GASES 

C 

WA TER S^aWA TER *L I MIT) -WATERED) /COS* THETA) 

C 

C FOR SATELLITE ALTITUDE CCHTUTAT IONS THE OPTICAL PATH CF WATER 

C VAPOR MUST INCLLCE THE INDEX CF REFRACT I CM CHANGES 

C 

IFaSHELL) WATERS = WATERS I MIT) -WATER*!) 

* ia yEHC — A T C E N T L I H I T) **ATCE.* > «T'i ! ) ) /COS* THETA) 

WATERCcdOT CENTAL I MI T) -OTCENTa I ) ) /CO&iTHETA) 

C 

c CONFUTE THE TRANSMISSIVITY FOR EACH LEVEL 

C 

10 TRANSKWI ) = TRAN5<tWAVEN0, WATERS pWATERC jWATERO, I > 

C 

c INITIALIZE THE UFWCLLING RACJATlCN 

c 


vhcle=o 



0014 4 

67* 

C 


Ctft 4 4 

66* 

C 

FOR THIS INTERVAL CCNFUTE THE EXT 1 NCI KM CF A SURFACE IMAGE. 

0014 4 

69* 

c 


DM 4 5 

70* 


SFCErC=TRANSMal ) *CCCCY«* WAVENO, TEMPERal ) > 

IV.T 14 5 

71* 

c 


0014 5 

72* 

c 

USING WATER VAPOR , CARDCN ClC'XlCE AND O^ZNE AS THE iNCCf ENOENT 

CKU4 5 

73* 

c 

VARIACLES INTEGRATE THE CCNTRICUTICNS CF EACH ATMOSPHERIC layer 

0014 5 

74* 

c 

70 THE IMAGE IN THIS INTERVAL 

00145 

75* 

c 


Ctfi46 

76* 


DO 20 I=2,LIMIT 

0*0146 

77* 

c 


00146 

76* 

c 

FI M3 T!£ CHANGE IN TRANSMISSION BETVCEM LAYERS 

0014 6 

79* 

c 


00151 

60* 


CEL TRAci I ) ^TRANSMxl >-7 kANSMuI -1 ) 

rtf! 51 

81* 

c 


00151 

82* 

c 

FIND THI AVERAGE TEKrERATUfiC CF A LAVER 

00151 

63* 

c 


00152 

84* 


TE W-CR=*TE>^Ek<i I > ♦ TE^ER*I -1 > > /2 . 

00152 

65* 

c 


00152 

66* 

c 

COMPUTE TX EMISSION FROM THAT LAYER 

00152 

67* 

c 


00153 

88* 


A TMOS= DDX Y*ANA VE N 3 , TEfPCR ) 

an 53 

69* 

c 


OQ153 

90* 

c 

WirLE IS TVt SUM CF FRCCUCTS aATK26PHERlC INTENSITY TUCS CHANGE 

00153 

91* 

c 

IN TRANSMISSION) 

00153 

92* 

c 


roi 54 

93* 


VHa.ErWHXE+«AATM:C^CCLTRA/tI) ) 

00155 

94-* 

20 

CCNTIMJE 

00155 

95* 

c 


OU 155 

96* 

c 

jnIEGRaTE LrW-LLiNG RA2IATICN FOR THIS INSTRUCT IT 

00155 

97* 

c 


03157 

98* 


inten-inttthc^vj cLE+sFcee)*REsrcmj> 

□0160 

99* 

11 

FCRMATctl x j F7 . 1 p ZX , F7 . 4 , 1 X»F6 * 1 »E10. 5 i 3X»F7 .2, 4X p F7, 2) 

00161 

103* 

30 

CCNTINLF 

00163 

lOl* 


CETECT= l NTEN* t NTER V 

00163 

102* 

c 


00163 

103* 

c 

EXTRACT A HATCH ING DETECTOR TEMPERATURE FROM THE CALIBRATION TABLE 

00163 

104* 

c 


00164 

105* 


CALL TEHTABdCETECT.DRTEM--) 

00165 

106* 


d_TErt==T£>FER*l) -DR TEX 

00165 

107* 

c 


00165 

108* 

c 

DISPLAY TEMPERATURE CCRRECTlCN AS A FUStCTlCN CF PRECIPITABLE WATER 

00165 

109* 

c 

HEIGHT CF TX CCSEKVATICN PLATFORM. 

00165 

1 io* 

c 


00166 

111* 


tJUTE«6,ll>FKESS*LlMlT) , WATERrtLlMl T) .XlGHTaUMlTl .HLMnLlMlT) , 

00166 

112* 


* TEHFER*LlMlT> ,DtTE>P 

00176 

113* 


X=DLTEKf= 

mi 77 

114* 


Y=HE1GHT**LIMIT) 

00200 

115* 


CALL LIXV*NXV*XF> , NYVdYF) ,NXVaX) r NYVciY)> 

00201 

116* 


XF=X 

00202 

117* 


YF=Y 

CXI 20 3 

116* 


RETURN 

00204 

119* 


ENO 


END C F CCWP ILATICN' 


NO DIAGNOSTICS 
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FOfc,S 1RANS/SLANU , TRANS , 1 RANS/SLAN11 
FOR 9 <4.-06/20-1 6’ SO <* Vi 1 1*1) 


FUNCTION TRANS ENTRY POINT IVX .,66 


STORAGE USED ' CCCC*tl> 0X510 CATAt/O) 0 , Xl03 ELANS. CCt-VOfriZ) 0l>_ rl i 


COHMX BLOCKS' 


0003 21 OMASA 

t»T4 72 002570 


EXTERNAL REFERENCES oEUOCK, NAPE) 


COOS 

PCXP 61 

0 X 6 

ALOG 10 

0007 

EXP 

0010 

I^RRSS 


«CLOCK * TYPE, RELATIVE LOCATION, KMC) 


STORAGE ASSIGNMENT 

OOOl 000447 UX. 

rooi 000146 5 SOL 

0003 R 001320 AESCRO 

0004 R 001750 O-kFSS 
0X4 ft 002114 EiLfi-5 
0004 ft G004 54 hCIGJfT 
OOOO 000056 lNJr$ 
0004 R 1X2424 G2CFC 
0004 R 000310 TEPEE R 
OOCVJ R 00.-015 THANSO 
OOO O R 000021 LEIGH 


OOOl 


000045 

1 OOL 

0001 


P00233 

fcV.V.X 

0004 

R 000764 

A I CENT 

0X0 

R 

000012 

D 

0003 

R 

001010 

EXFCI-O 

CUG4 

R 

cvxl20 

HU-1 

0003 

R 000310 

INTER 

oooo 

R 

0170005 

FRATIC 

0000 

fi 

OtXV.VXi 

TRANS 

COCO 

R 

000014 

TRAN5R 

0003 

R 

0L7JS44 

UEIGHT 


0 X 1 


0 X 266 

1217 .X. 

0 two 

R 

0 , 70 X 7 

A 

OL 7 .VJ 

R 

0 . 7 X 10 B 

0004 

R 

00220 * 

CEF ALL 

0 X 3 

R 

000 I 44 

£ 

£ 

Li 

row 

I 

O'OVTl 

I 

0000 

I 

00.7017 

J 

0000 

R 

OOOO.T 4 

■n 

X) 

> 

0 

0000 

R 

0 X 006 

TRANS C 

oooo 

R 

000016 

TRASLM 


0001 000400 176G 

0003 R 00L7.VJ0 ABSOR 
OOOG R 000013 BETA 

0004 r an 440 celtra 

0003 K 000C30 EXFCNC 
ro/io i 000020 n/:ex 
0000 R 000002 CNE 
XG4 r xxkx press 

0004 R 001 274 TRANSM 
0000 R 000003 TRATIO 


0001 


0*00353 

21 SOL 

0003 

R 

000620 

ABSTRC 

0000 

R 

0X00 1 1 

C 

0004 

R 

0X1CO4 

EPRESS 

0003 

R 

001154 

Exr;?^ 

COG€i 


0.7X52 

INJF* 

0X4 

R 

001130 

OT CENT 

0003 

R 

0004 54 

RESfON 

0001 

R 

0X471 

TRANSO 

0004 

R 

000144 

WATER 


OOlOl 

14 


00101 

24 

C 

oatoi 

3* 

C 

DD1Q1 

44 

c 

ocuai 

5* 

c 

DU101 

64 

c 

OOlOl 

7* 

c 

nni r t * 

ft* 


00104 

9* 


00105 

in* 


00106 

li* 


00106 

124 


00106 

134 


00106 

144 


00107 

154 


0D107 

164 


00107 

17* 



FUNCTION TRANSM ,W,w: p WD r LEVEL) 

THIS FUNCT ICN COMMUTES DC PERCENT TRANSMISSION THROUGH WATEit VAPOR 
IN TVE MANNEn SUGGESTED BY DAVIS ArC VIEZEE FCR DIE WAVE NU’CER RANGE 
25 y 0 2150 i/CH. IT IS SUPPLEMENTED BY MOSKALENKO TRANSMISSION 
FUNCTIONS FOR OZONE 

RFAf INTER 
PARAMETER 5IZFTH=100 
PARAMETER SIZTkA-ILO 

COMHCN/Z1/ ABS^rtSlZTRA) , EXPOWSIZTRA) , IMTER<*S121 RA) « 

* RESPON^SIZTRA) ,ADSD<C«10) , EXFCNCrID) , 

* LElGHTaSlZTRA) .EXFCWiftl ZTRA) ,EXFCN:«SIZTRA> , 

* ABSORCnSlZTRA) 

CCKMDN/Z2/ FkESS^SIZFTH) , WATERdSl ZPTH) , TEMFERctSlZrTH) , 
t HEIGHT^SI7FTH) ,HUt'V^5lZPTH) » ATCENT<*$I ZPTH) , 

4 OTCENTaSlZrTH) , TRANSAS I ZPTH) .CELTRActSl ZPTH) , 
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0010? 18* * EPRESS*SIZPTH> .ePRCSS'rStZPTH) .ETEMPS'iSlZFTH) , 

tv.mi? 19* * CeFAUL«5T2FTH) ,OZ‘:cC«iSlZr-1M) 

wuin 2C»* irans=i. 

Will 21* IFUW.LE.0.1 CO TO 215'.i 

OOlll 22* C 

will 21* C USE WAVEnumCer to compute an index 

com 24* c 

CO 115 25* I=*F-12.S)/25. 

Gtni4 26* CNE = 1. 

00114 27* C 

00114 28* C COMPUTE PRESSURE AfC TEMPERATURE RATIOS AT THESE ST AMI ARC VALLES 

00114 29* C 

00115 50* TRATICC1* 

00116 51* PRAT 10=1. 

00117 32* PRATIC=1. 

CO 1 2'J 33* iFctLEVQ- .EO.l) CO TO 100 

00122 54* TRATlC*=‘iETEFr^EvEL>+ETEKrSrxLEVLL-l>l/546.52 

00123 35* f RATIC^EPRCSS^LEVEL) + EFRESS' , 4XVEL-1 > ) /2V326. 4 

00124 36* PRAT lC=aCFRESS'iLEVEL>*CrRESS*LEVELM> 5/2026.4 


00125 

37* 

ICO 

corn 1 HUE 

00126 

3a* 


IF* I .GE.32) CO TO &JO 

00130 

39* 


TRANS C=l. 

00131 

40* 


lF<x 1 4 LT .22) CO TO 550 

00131 

41* 

C 


DJJ31 

42* 

c 

COPE LTTE TRANSMISSION DUE TO CARDIN DIOXIDE 

0L1131 

43* 

c 


00133 

44* 


A- . 4 -<x . 1 5 ) *ALOC 1 'VtP RA T I C ) - 1 .6 

00134 

45* 


ADSCRCn I -2 1 ) . 6 7* aPRAT IC**a . 8 J ) ) *PRAT I C**A 

D.U35 

46* 


C=TRAT lO^JOtFCNC* J -21 ) 

00136 

47* 


D-a-4) ^tPRAT . 8 ) ) 

DU137 

46* 


TRANSC-EXFa-ciC*^iriC*C-H>C)**.53 -aE) ) ) 

00137 

49* 

c 


00137 

50* 

c 

SEARCH AM? INCLUDE OZCr£ TRANSMISSION FOR THIS INTERVAL 

00137 

51* 

c 


CX-il 40 

52* 


TRANSC=TlRANSC*TRAr^Oi 1) 

OO 14Q 

53* 

c 


00140 

54* 

c 

COT LTTE TRANSMISSION CUE TO WATER VAPOR 

00140 

55* 

c 


00141 

56*. 

550 

continue 

DJI 42 

57* 


BETA= 4 76+<i. 5a+« # 48*aPRATlO**2) ) )**. 5 

00143 

sa* 


A^ASSOkciI) ^Wi^TRATIOl^EXFONal ) ) 

00144 

59* 


B-l .♦^.17 *CETA««FkATIO0‘4«-1»1))*A> 

00145 

60* 


C=BETA4«FRAT lO“«- . 1 ) ) *A 

DO 146 

61* 


TRANSR-EXFa-C**C**<t-*S) ) ) 

□0140 

62* 

c 


00146 

63* 

c 

MULTIPLY THE CAREON-D IQX ICE AMI WATER VAPOR TRANSMISSIONS 

00146 

64* 

c 


DJI 4 7 

65* 


TRANS^TRANSR*TRANSC 

COi so 

€6* 


CO TO 2150 

J-* 

w“ 

c 


DJI 30 

cs* 

c 

COMPUTE TRAN SMI 33 1 Of 4 IN T* IE WATER VAPOR *W3J«w»Ca«»* . 

OOI 50 

69* 

c 


OOl 51 

711* 

800 

CONTINUE 

OO! 52 

71* 


iFal ,GE,48> CO TO 120C1 

DJI 54 

72* 


TRANS=EXP*-«AGSCRrt I ) *W*FRAT I O) **CXPCN<t I ) > 

00154 

73* 

c 


00154 

74* 

c 

SEARCH AND INCLUDE OZONE TRANSMISSION FOR THIS INTERVAL 
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00453 

315* 


03 31 1 J=1 ,64 

0U4 5G 

316 * 


1FUNC T = *RC SrCMa J ) *1 00 *1+46 

00457 

31 74 


JFUNCT = uRCSFCN'i J -♦ 1 ) *1 00 . 1 ♦ < B 

00460 

316* 

311 

CALL L I NE Vti920* J * 1 FUNC T ,921 + J , JFI^CT > 

004 62 

319* 


Y=0. 

0*0463 

320* 


03 312 J=1 »5 

004 66 

321* 


V=YMOOOO. 

00467 

322* 


X=ll . 

004 TO 

323* 


XF=-2. 

00471 

324* 

312 

CALL LlHEVuNXVaXP) ,NYV«V> .NXVdX) .NYVuTl) 

1X1473 

325* 


HEADERS) =- • 

CX14 74 

326* 


ENCODE '*2 4, HEADER) ANCLE 

00477 

327* 


CALL PRINT Veil 6 .HEADER ,703,64) 

00500 

328* 


HEACLR-t3) = - “ 

00501 

32** 


ENCC0Cti25 , HEADER) TElTERtti) 

00504 

330* 


CALL F R 1 HI V«U 8 , HE ACER , 703 , 80 > 

00505 

33t* 


HCADERalJr^TE^ C“ 

00306 

332* 


HEADER^ 2) -"‘ORkECT - 

00507 

333* 


WIACER<0> = -ICN C“ 

aisto 

334* 


CALL FR I NT Va IB , HEADER , 500 , 60) 

00511 

335* 


HEADERS 1 > ** ALT I TU“ 

00512 

336* 


HEADERS) = *CE* 

0031 3 

337* 


FCACEK^3) = “fEET** 

0051 4 

336* 


CALL APRNJ VaO ,-14,18, HEADER , C 0 ,9'Xi) 

0051 5 

339* 


fcadgwi)=*suvFa<:- 

00516 

340* 


►CADEk«2) = <, E‘’ 

00517 

341* 


CALL PRlNTVri7 r HEADCR,60,U0) 

00520 

342* 


FCACEk^.I ) -”1 nooo - 

00521 

343* 


CALL FRlNTV</5, HEADER, 70 ,252) 

00522 

344* 


HEADER n J ) - fc >.0t V.V.t " 

CXI 523 

345* 


CALL PR I NTV<i 5, HEADER, 70 ,393) 

00524 

346* 


hEACEA^ 1 ) =*30000- 

00525 

347* 


CALL FRlNTVuS, HEADER, 70, 535) 

00526 

346* 


HEADER «tl ) — * 40000 - 

00527 

349* 


CALL FR 1 NTV«i5 , HEADER , 70 , 670 ) 

00530 

350* 


HEADER^ l ) = - 3LXXV.1- 

00531 

351* 


CALL FR I NTV<iO, HEADER ,70,81 8) 

00532 

352* 


HEADER til > = - 60000" 

00533 

353* 


CALL PRINT YdS ( HEADER , 70,950) 

00534 

354* 


CALL CHSIZV*3,5) 

00535 

355* 


CALL R 1 TE2 Val 00 , 1 000 , 1 OOO , 90 , 1 , 60 , 1 t FCSSAG , NL AST ) 

00536 

3 56* 


THCT A = RA D* ANGL E ) 

00537 

357* 


LIM-LIMIT-I 

00540 

35** 


CALL CFATHS 

00541 

159* 


LlwrT=i 

00542 

3 CO* 


DO 20 J=1,LIM 

0054 5 

361* 


LIHIT=LIHJT+1 

UU545 

*tv"* 


CAll RADHTC 

nil «i7 


?r, 

rr-+tr Tmi F" 

OG55t 



CO TO 100 

00551 

305* 

C 


00551 

3t>0* 

C 

i m‘ ut resfcwsc functions as percentages at midpoints cf each 

00551 

367* 

c 

interval 

00551 

3G8* 

c 


00552 

3C9* 

200 

CCNTINLE 

00553 

370* 


LFACER«1 ) -HEACl 

00554 

371* 


LEACrRa2)=HCAC2 


i 
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00555 

372* 


00 3 56 

373* 


00360 

374* 


00363 

375* 

30 

00565 

376* 


00 5 TO 

377* 


00374 

376* 


cogoo 

379* 


00601 

360* 

40 

00603 

361* 


00603 

302* 

C 

00603 

303* 

c 

00603 

3G4* 

c 

00003 

305* 

c 

00604 

306* 

41*0 

00607 

307* 


00611 

300* 


00614 

389* 


00615 

390* 

500 

00616 

391* 


00617 

392* 


cr »620 

393* 


00620 

394* 

c 

00620 

395* 

c 

00620 

396* 

c 

00621 

397* 

coo 

00627 

390* 


oocsa 

399* 

*00 

00630 

4 OCT* 

c 

00630 

401* 

c 

00630 

402* 

c 

00631 

403* 


□0633 

404* 


00635 

405* 


00643 

406* 


00651 

407* 


00654 

408* 


00062 

409* 


00662 

410* 

c 

00662 

411* 

c 

00662 

412* 

c 

00670 

413* 


0O6T3 

414* 


00675 

415* 


D0676 

416* 

815 

00676 

417* 

c 

00676 

410* 

c 

00676 

419* 

c 

00677 

420* 


C2Z77 

*21 


00677 

422- 

r 

OOG77 

423* 

C 

00677 

424* 

c 

00701 

425* 


00701 

426* 


00701 

427* 

c 

00701 

420* 

c 


LEADERS) =HEAC3 
WU1E«6,19> 

do 30 i=i,e3 

RESPCNit > = 0, 

DO 40 1*1 .LlHll 
REACTS, 15> INTERV , FUNCT 1 
W<nea6,16) INTERV , FUNCT I 
J=«INTCRV-12.5)/2S. 

REST<NrtJ>=FUNCTl 
CO TO ISO 

I^fUT THE TEMPERATURE RAMSC OF TFC BLACK-OX Y INTENSITY/ 

temperature scale in cent i or ace 

REAOl 5 1 1 5 ) TRANCE 
WRlTEaS.16) 

Lfi !TEa6 , 1 6) TRANCE 
CO TO 100 
CONTINUE 
CALL PLTtOiO) 

CALL PL TERM 
STOP 

IISPLTT DISPLAY LABELS 

REAT>i5,26)o^eSSAG<*I) ,1=1 ,10) 

CO TO ISO 
CONTINUE 

READ DEFAULT PROFILES ANC INTERPOLATE MISSING DATA 

IF<iSTEFS2.GT.O> CALL STEFCR 
WRITE*6,35> 

READ* 5 , 33 , E>C= 1 CO) PRESS P ,QZ&&r , TEPFkE , NUMI Or 
WUTE*6,33> PRESSr,02O'£P,TE^RE J Ha^I0P 

DO 620 J-2 .LIMIT 

REACTS ,33 ,EM3=100> PRESUH ,0 ZC*CS , TE>TGR , HUMIC 
WRITERS, 33) PRESUR ,GZ'INES , TEPFOR r HUM ID 

SCAN REQUESTED LEVELS FQk DEFAULT PRESSURE WINDOW 

DO 810 1=2 .LEVELS 

IF<iFRESS*I) -LT.PRESSP. AND. PRESSED *GE«FRESLI5) CO TO 015 

CO TO 010 

CONTINUE 

INTERPOLATE hissing tew-ekature 

I Erf TEMPER** I ) .CT*10 J J0.) TDEERf* I } = Tuvr^RCSSr ,f kE5LR i TthPRE , Tl^POu , 


INTERPOLATE MISSING SPECIFIC HUMIDITY 

IF*FLDa3 t l ,CEFAUL*I > ) .Ed, 1) HUM* I) =T ERFafRESSP » PkLSUR , HUMI CF .HUM I D 
♦ , PRESS* m 


rtw 

NEW 


interpolate MI5SINC 02 one 



TO 701 

4294 

c 


TO 70 3 

4304 


tF«FLCa2,l r CGFAUL*n> .CQ.l) 02CNE* 1 I -7EKP4FRCSSP f FRESUR ,OZCNEF « 

00703 

4314 


* OZONES, PRESSED ) 

00705 

4324 

010 

CCNllNUC 

TO 70 5 

4334 

c 


00705 

4344 

c 

FETCH 1CXT LEVEL CEF AULT PARAMETERS 

00705 

4354 

c 


TO* 71 *7 

4 364 


FRESST-zFRESUR 

00710 

4374 


OZCNEF^CSCNES 

0071 1 

43** 


TEHPRC=1EMPCR 

00712 

439* 


HUHJCr=HLM|C 

00713 

440* 

620 

CONTINUE 

0071 5 

4414 


CO TO 100 

00715 

4424 

c 


00715 

443* 

c 

lU'KJX TRACE CAS TRANSMISSION FUNCTIONS WHERE 

0071 5 

4444 

c 


00715 

445* 

c 

V^IICHT TRUNCATED TO AN INTEGER IS AN IMCEX TO THE UAVENL7CLR 

am 5 

4464 

c 

INTERVAL FCR VHICH THIS ANALYTIC FUNCTION APFLYS 

00715 

447* 

c 


con 3 

44*4 

c 

VEICHT FRACTION IS THE CrCCFETRIC ►CAN WEIGHTING KCIFIED DY THE 

0071 5 

4494 

c 

FLANCKIAN FUNCTION FOR THE WAVE -NO EER INTERVAL FRACTION 

00715 

4 504 

c 

OF THIS ANALYTIC FUNCTION 

00715 

4514 

c 


□0715 

4 52* 

c 

ADS-3RO - A THESE ARE PARAHCTERS APPEARING IN ANALYTIC 

OJ715 

4534 

c 

EXFCMj r B FUNCTIONS OF HE FORM A^iATrtj6-CENT*^C) tnEFttO 

00715 

4544 

c 

EXFCND = C FOR THE INTERVAL ANC WEIGHTING CESCRlEED A BONE 

00715 

4554 

c 


TO716 

456* 

1 1GCJ 

continue 

00717 

45/4 


W<lTEa6,3l) ►CADl,HEAC2 r lCAD3 

00724 

4564 


Do mo ;=i , limit 

00727 

4594 


iFetf .GT.1TO) 03 TO 100 

TO 731 

4 on* 


REACTS , 2ft , ETC =1 TO) WEIGHTS) , ACSORCtt 1) . EXPCMM I ) , EXFCTO* I) 

00737 

461* 


1 I=Vs£I GHTa I ) 

00740 

462* 


LElGH^veiGHTaI)-II 

00741 

4634 


BEGINS INTERIM) -12. 5 

CTJ74 2 

4644 


FINl S= INTER* 1 1) *12 , 5 

00743 

4CS4 


VftlTE«G,29) CEGlN^ims.AESlROil) jEXFCFCuI) |EXrCfO*l) , WEIGH 

00753 

466* 

1110 

CONTINUE 

0075 5 

467* 


CO TO 100 

00756 

466* 


function SPHERE**, L) 

0)756 

4694 

c 


00756 

4 m* 

c 

THIS IS A CONDITIONAL FUNCTION WHICH ►CCIFIES THE PRECIPITADLC 

00756 

471* 

c 

WATER AT SATELLITE ALTITUDES FOR THE SPHERICAL SHAPE CF THE 

00756 

472* 

c 

EARTH 

00756 

473* 

c 


0O761 

474* 


sphered. 

00762 

475* 


IF4.NDT.SHELL) RETURN 

007fc> 

4 i '6* 

c 


ivr»e-» 

4 '»'»* 

c 

rrct To tlr- Miir r*~ •rx.nr iik^i r rn^rmrir t#-*. t/> tut 

nr)76? 

4 7ft* 

c 


00764 

479* 


TE S T = 21 iS 9 ft 69 6 . /42O090 69 6 . ♦ HE I GH T<iL J ) 

00764 

4ft0* 

c 


00764 

4014 

c 

N> CONFUTATIONS ARE MACE WHEN THE FIELD OF VIEW IS ADOVE THE 

00764 

462* 

c 

HORIZON 

00764 

403* 

c 


00765 

4644 


IF«TKETA ,GT . ASlNctTEST) ) RETURNl 

00765 

405* 

c 




0076* 

466* 

C 

MICH IS THE AVERAGE HEIGHT CF THE WATER YAFOR LAYER 

00765 

467* 

C 


GG767 

466* 


HIGH=*KEIGM1«L> ♦HEtGH7*L-l I ) /2. 

00767 

»:69* 

C 


00 767 

490* 

c 

ADO THE AVERAGE HEIGHT CF THE LAYER TO TIC RADIUS CF THE EARTH 

00767 

49 1 * 

c 


00770 

492* 


FART 1 =M I GH4 200911 696. 

00770 

493* 

c 


00770 

494* 

c 

ML1_T 1 T L Y DY THE MEAN IhCEX CF REFRACTION AT T»C AVERAGE HEIGHT 

GO 7 70 

495* 

c 

CF THIS LAYER 

00770 

496* 

c 


OJ771 

497* 


FART2-PART1 *iREFRAXciL > +RCFRAX<iL-l ) ) /2, 

00771 

49** 

c 


00771 

499* 

c 

MULTIPLY E.Y THE l^EX CF REFRACT Km AT THC HEIGHT CF TVC SENSE* 

C7'»77l 

500* 

c 

AND BY TFC Site CF THE ANGLE CF CGSEnVAT ICN 

0.1771 

501* 

c 


00772 

502* 


FART3=FARTl*RCFRAX*L)*SlNiT*£TA> 

00773 

503* 


PART4= S*2RT«af= ART24FART?) -*FART 3*PART3) ) 

00773 

504* 

c 


0.1773 

505* 

c 

GOTLTTC THE h€AN FUNCTION ASSUMING T>C CE SERVING SENS.CR IS JUST 

00773 

506* 

c 

AD9VE TVC FCAN LAYER TOViUCH IT IS APPLIED 

□0773 

507* 

c 


00774 

50 * 0 * 


SPHERE = FART 2/FART 4 

00775 

5*09* 


RETURN 

00776 

5 ID* 


sucrclttinc cistly 

00776 

511* 

c 


00776 

512* 

c 

THIS SLCROUTtKE CIST-LAYS TIE MXIFIEC DATA BASE 

00776 

513* 

c 


OlOJl 

514* 


VRITTa6 f 1 1) 

01003 

515* 


OO ID l-t,LCVELS 

01006 

516* 

10 

W\1 TEtfQ, 12) PRESS* I ) r Mi>Ul) ,TEM ER*I) ,OZCtC*I) 

01015 

517* 

11 

FORMAT*- 1 PRESS SPECIFIC TE^F C VCt£ m ,/ t 

01G15 

516* 


* • HLHICITY SQU>ClNG- r / r 

01015 

519* 


* ■ ►© GM/KGM K GM/KGM N > 

01016 

520* 

12 

FORMAT «1X,F7.1 ,lX,El0,5|2X,F6.Z f lX i E10.5) 

01017 

521* 


L I MI T=LEVELS 

01020 

522* 


RETURN 

01021 

523* 


FUNCTION TERFaXl , X2 p Y1 , Y2, X) 

01021 

524* 

c 


01021 

525* 

c 

THIS FUNCTION PERFORMS LAGRANGIAN INTERFOLAT ICN 

01021 

326* 

c 


011124 

327* 


SLCFEA^*X-X2> /*X1 -X2) 

01025 

520* 


SLCFEE^X-Xl ) /aX2-Xl ) 

01026 

529* 


TERF = SLCF EA * Y 1 + SLCF EB* Y2 

01027 

530* 


RETURN 

01030 

531* 


sic^outim: ct-aths 

01030 

532* 

c 


c;07.s 

r.33* 

c 

tniS iLUxCLH If*. CO’TUlLii ANCLE DEPENDENT OPTICAL F ATH5 

oi rvsn 

*L-*jt* 

r 


01033 

535* 


ruin r^o. 

01034 

536* 


TWHCLE=D. 

01035 

537* 


OO 10 I-2.L1H1T 

01040 

536* 


OTCEN T* 1 ) "OTCLNT* 1-1) ♦uFRESS* I - 1 ) - CRESS* I ) ) **OZCfJEtf I ) *OZCNCaI -1 > ) * 

OJ 040 

539* 


* 237.960 

01G41 

340* 


WATERED rWATER* I -n ♦ aHU Vi I > *STMCREti$lO , I J *oFRESSa | -t ) -PRESS*! )}) / 

01041 

541* 


* 900. 

01042 

542* 


FRESUR^FRESS*! ) *FRESS*I-1 )>/2* 



Ill 043 

3434 


TEMPCR=aTEMrtRd I ) ♦ TEMPER* I -1 )) /2 . 

01044 

3444 


lF«l.GT.n 

01044 

5454 


I CHTa n =NC I GMT* l -1 ) *29 .3*TEHPCR*AL0G*PRE$Srt [ - 1 ) /PRESS* 1 1 ) 

01044 

3464 


***!.♦. 0006 1 vJHUKc I)J 43.2808 

01046 

3474 


ATCENl4l>;<iF*ESS'in-FRESS'iin«.2G0 

01047 

3464 


EUi:CE=EW*XE*«WATER*U -WATER*] -1 > ) *FRESUR 

01050 

3494 


TU ac= TW 1XE ♦ UWA T ER«i !> -UA TERu 1 - 1)) *TEM‘<R 

01031 

3304 


CPRESS*I) =*FRE5S*1 ) ♦PRESS* I > >/?. 

01052 

5514 


EFKCss^n =ew4:le/vatcr*i> 

01053 

3524 


ETEHTS^I) -TVHX E/WAT ERaU 

01054 

5534 

1G 

CCNTINUC 

01056 

3544 


RETURN 

01U57 

5554 


SUCRCVTISE STEFER 

01057 

5564 

c 


01057 

5574 

c 

THIS SLCROUTlNt EXP ATCS 0\ CONTRACTS TVC DATA TO TPC KCECCD STEP SIZE 

01057 

356* 

c 


□ 1062 

5594 


j=i 

01063 

5604 


LCD I NT=HI NrtFRESSrf 1 ) /STEP SZ . $1 ZPTM J 

01063 

561* 

c 


01063 

562* 

c 

TEST FOR ENC CF AVAILABLE DATA 

01063 

563* 

c 


01064 

564* 

5 

IFdj.6T.LCCl NT) RETURN 

□1064 

5654 

C 


01064 

5664 

c 

SET AVAILABLE CATA 

01064 

567* 

c 


01066 

568* 


PRESSP-FRESSaj) 

01067 

569* 


FRESUT=FRESSrtJ+n 

01070 

5704 


TErFRE=TEK'-ERdJ) 

01071 

5714 


TE>PCR=TElf'CRrt J f 1 ) 

01072 

5724 


OZO^F^OZataj) 

01073 

5734 


ozocrs^o/ci^d i ) 

01074 

574* 


HUMlDF=HUMdJ) 

01075 

5754 


HLMID=KJPMJ + 1) 

01075 

576* 

c 


01075 

5774 

c 

SET WANTED DATA 

01075 

5784 

c 


01076 

579* 


FREWAN^FRESSdl } -STEFSZ*J 

01076 

5604 

c 


01076 

581 * 

c 

IS WANTED DATA WITHIN AVAILABLE DATA 

01076 

5624 

c 


01077 

5834 


lFdFREVAN.LT .FRE55F . AFC .FREWAN .GE .PRESUR) GO TO 10 

OllDl 

5644 


GO TO 40 

01102 

565* 

IO 

CCNTINUE 

01102 

566* 

c 


01102 

567* 

c 

HCW KANY LEVELS CF WANTED DATA ARE WITHIN WlhODW 

01102 

566* 

c 


01103 

569* 


LCDSIZ=dPRESSF-PRESUR) /STEFSZ 

01103 

590* 

c 


Cl 103 

591* 

c 

KAKE SPACE If; AVAILABLE DATA TCR WANTED DATA 

01103 

592* 

c 


01104 

593* 


DO 20 K=LEVELS. J»~l 

01107 

394* 


FRESSdK+LCCSl2)=FRESSdX) 

OHIO 

5954 


TEHrERdK+LCCSlZ) -TEMPEfidK) 

01111 

59G* 


OZCNCaMLCC5lZ)=OZONE*K> 

0111? 

397* 

20 

HUHdK ♦ LCCS I Z) -HLFW K ) 

01114 

396* 


LE VEL S =L E VEL 5+ LCD S f Z 

01114 

599* 

c 




01114 

too* 

C 

load data cner wanted interval 

01114 

601 * 

C 


Gil 15 

602* 


DO 30 K,~ 1 >LCC SI Z 

01120 

603* 


J= J*t 

01121 

604 « 


TElFERa J) - TERFcir-RE 5ST- , FRESUR , 1 EH“RE , TfMT«t , F REWAN) 

01122 

603* 


02C*JE«J)^trRF*rRES5r .FRESUR .OZuNEP ,OZ'jf«'ES »F*REWAN) 

01123 

606* 


HU^WJJiTERFuFRESSF.PRESLR.HUMlCPiHUMlD.FREV^N) 

01124 

607* 


fRESS* J)=FRCWAN 

01125 

60S* 


FREWAfs^rRCWAN-STErSZ 

Oil 26 

600* 


CEFAUL'tJl -DEFAULT J“l) 

01127 

610* 

30 

CONTINUE 

01127 

611* 

C 


01127 

61 2* 

c 

RESET BOTTOM LEVEL CF AVAILABLE DATA 

01127 

613* 

c 


01131 

614* 


CO TO 5 

01132 

615* 

40 

CONTINUE 

01132 

616* 

c 


01132 

617* 

c 

CESTRi^r TOP LEVEL OF AVAILABLE DATA 

01132 

616* 

c 


01133 

619* 


DO 50 K-J .LEVELS 

01136 

620* 


FftESS'iKn ) -PRESS tK+2) 

01137 

621* 


TOPER*** 1 > =TT>r-ERaK+2) 

01140 

622* 


OZcrCaK ♦ 1 ) =OZ>C a* ♦ 2 ) 

01141 

623* 


HUMcUC*l)=HU*K+2> 

01142 

624* 


CGFAULaKM ) = DEFAULT 2) 

01143 

625* 

50 

CONTINUE 

□ 1143 

626* 


LEVELS: LEVELS- 1 

□1146 

62/* 


GO TO 5 

01147 

626* 


EM? 


DO CF CCM’ILAT ICN ' MD CIACNDGTICS. 


FRCP TFr$ . 

FUFPUR C/25D--OB/20-18 ’ 50 



001*4 

73* 

C 


0015* 

76* 


TRANS=TRANS*TRANSOuI) 

00156 

77* 


CO to 2150 

00156 

76* 

c 


00156 

76* 

c 

COFl/TE TRANSMISSION CUE TO WATER VAPOR, 

00156 

60* 

c 

- 

OCU37 

61* 

1*00 

CONTINUE 

00160 

6 ?* 


eCTA=l,16* 4l.36+,46*«FRATlC»**n**.5 

W161 

63* 


A=ABS:R«l>*W 

00162 

64* 


e=BETA**RRATtOtt«- .1 5) ) *A 

00163 

•3* 


Cb« 1 . ♦M .94CET A4^PRAT tOtta-1 . 1 5)) *A> ) **u- , 5) 

00164 

66* 


TRANS£CXFa~«&*C> > 

00164 

67* 

c 


00164 

66* 

c 

SEARCH AND INCLUDE C*C*C TRANSMISSION FOR THIS INTERVAL 

00164 

66* 

c 


00163 

90* 


IF*! .CC .76) TRANS" TRANS* TRANS^i I ) 

00167 

91* 

2130 

RETURN 

00170 

92* 


FUCTICN TRANSOM? ) 

00173 

63* 


TRASLM=0, 

00174 

94* 


TRANSOM 1 . 

00173 

95* 


DO 10 J=l,100 

00200 

96* 


I«€X=WEIGHT«*J> 

00200 

97* 

c 


00200 

96* 

c 

USE DATA WITHIN T>C CURRENT WAVEMJ-CER INTERVAL 

00200 

90* 

c 


00201 

ICO* 


if« i hcex.se. i) co to to 

00205 

toi* 


UErO^V£lGHTaJ)-INCeX 

00203 

102* 

e 


00203 

103* 

c 

OOFUTE MC€KALl>KO TRANSMlSSlCN FUNCTIONS FCR OZC*E USING TRACE 

00203 

104* 

c 

CAS COEFFICIENTS 

00203 

1P5* 

c 


00204 

106* 


B=ADSCRCM J ) **W5#CXrCHC« J > ) **FRAT IC**CXPC*£u J) > 

€0204 

107* 

c 


00204 

toe* 

c 

APPLT PLANUAN WEIGHTING FUNCTIONS TO TVE TRANSMISSION VALLES 

0U2O4 

109* 

c 


00203 

no* 


TRASUM=TRASLM^XF*-B) ICH 

00206 

m* 

10 

CONTINUE 

OGZlO 

112* 


IF*TRASSUH.CT.0. ) TRANSOM TRASUM 

€0212 

113* 


RETURN 

00213 

114* 


oc 


DC CF OOK=IL>TICN 


NO DIAGNOSTICS 



88 


FCft.S TRANS/SLANT, TRANS, 1RANS/SL ANT 
FCR 9*L-l'>*/2V»«1*'5» 


ruNCTlCN TRAUS ENTRY POINT CW3 r ,6 


a TOR ACC USED' CCCEM) 000510 DATA-O) 01X1105 CLANU CCHrCNtfZ) OW.7.7.) 

CC**»X CLOCKS* 

0005 21 001464 

CK XU 22 015530 


EXTERNAL REFERENCES -CLOCK, NATO 


0005 

KEXF61 

0006 

ALOGIO 

0007 

EXP 

OOIO 

fCRR3$ 


STORAGE ASSICNKENT -CLOCK, TYPE, RQ.ATIVE LOCATION, NATO 


0CX31 


000447 

1DL 

OGGI 


000.145 

IDOL 

0001 


0/0266 

12COL 

00.11 

000400 

17GG 

ODOl 


000353 

2150L 

oocu 


000146 

550L 

tXTOl 


000233 

600L 

oral 

R 

000007 

A 

00,13 

R 00/000 

aescr 

roo3 

R 

000620 

ad&:r 

0003 

R 

001320 

asscro 

0004 

R 

004 P.14 

ATCENT 

axvj 

R 

ixv/oio 

0 

a 000 R 0,7/013 

BETA 

0000 

R 

OOOGll 

c 

0004 

R 

011610 

CPRESS 

ooao 

R 

000012 

D 

0004 

R 

Q135GG DETAIL. 

O.X14 

R 007640 

DELTRA 

0004 

R 

010624 

EPRESJ 

0004 

R 

012574 

ETELFS 

0003 

R 

tXUCllll 

FXFOO 

0/03 

R 

root 4— 

Pvk-N 

0O/.3 

R 1*00*632 


COO: 3 

r. 

ni ** c* 

Exrcr:: 

0004 

R 

□02734 

ICIGKT 

CGC»4 

R 

003720 

HUH 

0000 

1 

0/0001 

l 

CGOJ 

I 000020 

ISCEX 

ccoo 


001X152 

IN IPS 

OUOO 


000056 

INJP$ 

OOC33 

R 

00031 0 

INTER 

oj:g 

I 

000017 

J 

0000 

R 000002 

ChE 

0004 

fi 

ry/issTO 

OTCEN1 

DCX14 

R 

014544 

Q7C»E 

DOOO 

R 

000005 

FRATIC 

0000 

R 

0/0004 

PR AT IO 

0004 

R 0/0000 PRESS 

0003 

R 

000454 

RESPCi 

0004 

R 

001750 

TEMPER 

GOOD 

R 

0000L73 

TRANS 

axw 

R 

0/0006 

TRANSC 

O.X14 

R 17/6654 

TkANSM 

0001 

R 

000471 

TRANS 

OOGO R 001X115 

TkANSO 

0000 

R 

00001 4 

TRANSft 

0000 

R 

0/0016 

TRASUJM 

0000 

R 000003 

TRATIO 

0004 

R 

000764 

WATER 

GOOD R 

000021 

VE1GM 

0003 

R 

00)644 

WEIGHT 













00101 

1* 


OOlOl 

2# 

C 

txuoi 

3* 

C 

0Q1D1 

4* 

c 

ooiai 

5* 

c 

OOlOl 

6* 

c 

OOlOl 

7* 

c 

0Q1D3 

e* 


00204 

s* 





00105 

114= 


oaios 

124 


00106 

13* 


00106 

14 * 


00107 

15* 


00107 

16* 


00107 

17* 



FUNCTION TKANSOF.W.LC.WD, LEVEL) 

THIS FUNCTION COMMUTES T>C PERCENT TRANSMISSION THROUGH WATER VAPOR 
IN TVe KANMR SUGGESTED CY DAVIS Arc VIE ZEE FOR THE WAVE NUK5ER RANGE 
25 TO 2150 1 /CM. IT IS SLFFLEPZNTEC BY MOSKALENKO TRANSMISSION 
FUNCTIONS FOR OZCKE 

REAL INTER 

FARAPCTER SI ZF TIL: 500 


I™ A«vAi C i 

CCfCt2fL / Zl/ ADSORB S I ZTuAJ , CXTCt&SlZTRA} « I NTCR^SI ZTRA) , 

* RESFCN-SI27RA> .ABSCKC-lG) , EXPCf4C*lC0 , 

* VEIGHT-SIZTRAJ , CXFOMCWS IZTRA J ,EXPC*0-$IZTRA> , 

* ABSCKO-SIZTRA) 

COMM3N/Z2/ PRESS-SlZPTH) , WATER-5 IZFTM ) , TEMPER- SlZPTH ) , 

* HEIGHT-SIZPTH) .HlAWSTZFTH) , ATCENT-SIZP7H) , 

* OTCENT-SIZFTH) , TRANSHaSlZPTH) ,CELTRA-5lZFTH) , 



00107 

164 


4 EFRESS'iSlZf TH) , CPRE$&«S1ZPTH> , ETEMrS^St ZP1 H> p 

00107 

194 


* CEFAULaSIZFTM) .OZCNC'iSl ZPTH) 

001111 

204 


TRANS=1 . 

CHHii 

214 


IFdW*LE*0.) GO TO 2150 

CXH 1 1 

22* 

c 


001 it 

23* 

c 

use wavenumber to compute an index 

CAM 1 1 

244 

c 


00113 

254 


IrdF-12. 51/25. 

00114 

264 


OC=l. 

00114 

274 

c 


00114 

264 

c 

OTUTE PRESSURE AND TEMPERATURE RATIOS AT THESE STAfCARC VALUES 

om 1 4 

294 

c 


oot 1 5 

304 


TRATIC^I * 

00116 

314 


FRATlQ^l . 

OOl 1 7 

324 


FRATIC=1 . 

01120 

334 


I F<i4_EVO_ * EO . 1 ) GO TO 100 

01122 

344 


TRAT lOcdETEf^ SdLEVTX) +ETE>FSdLEVEL- 1) > Z546.32 

00123 

354 


PRAT lOroEPRE$Sr4LEVEU ♦EFRESS^LEVEL" l ) > /2V-26.4 

01124 

364 


PRAT iC^tCFKESSoLEvTL ) ♦CFRESSoLEVEL- 1 ) > /?'*26.4 

00125 

374 

100 

CCKT1NIC 

00126 

364 


1 Fd I . GE . 32) GO TO 600 

01130 

394 


TRANSC=1. 

00131 

404 


IFdl .LT *22) CO TO 550 

00131 

414 

c 


00131 

4Z4 

c 

COf^UTE TRANSMISSION DUE TO CARDIN Cl OX ICE 

□0131 

434 

c 


00133 

44* 


A= . 4 -4 . 1 5 1 *ALOG 10*<PRAT IO -1 . G 

00134 

454 


AESORCa I -2 1 ) <VJC*4 . 67+ «PR AT I . 6) ) I 4FRAT l C*4A 

01135 

464 


C=TtvAT iCX^CXFCNCd 1-21) 

01136 

474 


C^*.4)*r<PRATlC<^6) ) 

00137 

464 


TRAMSC-EXPa-aC^.tr^C-^CNE) 4* . 5) KIG ) ) 

00137 

494 

c 


00137 

504 

c 

SEARCH ARC INCLUDE C^'tRE TRANSMISSION FC* THIS INTERVAL 

00137 

514 

c 


00140 

524 


TRARLSC'TkANSO*TkANSO-^I) 

□0140 

534 

c 


00140 

544 

c 

commute transmission cue to water vapor 

01140 

554 

c 


00141 

56* 

550 

CONTINUE 

DPI 42 

574 


ECTA=. 76+d. 58+d. 4 8*dFR AT ICK^2>))**. 5 

00143 

564 


A=ABSORd I ) <VA^TRAT lO:^CXPOmi > ) 

00144 

594 


6=1 . +a3 . 1 74CET A*dFR AT I 1. 1 > > *A > 

00145 

604 


C=E£TA*«xPRAT lOi=^i- . 1 >) *A 

00146 

614 


TRANSR=EXFd-C*dC**u- *5 ) ) ) 

00146 

62* 

c 


00146 

63* 

c 

MULTIPLY THE CARCCN-C IOXICE AfC WATER VAPOR TRANSMISSIONS 

00146 

64* 

c 


00147 

65* 


TRANS= TR ANSR * TR ANSC 

IX*13U 

€C* 


CO TO 2150 

rv-»« 


r 


nni 

6** 

c 

COMPUTE TRANSMISSION IN THE VtMEK VAPOR "WIJ33U“. 

00150 

69* 

c 


00151 

70* 

aoo 

CONTINUE 

00152 

714 


IFdI.GE.46) GO TO 1200 

00154 

72* 


TK •‘.NS=EXFd-<xADS3Fid I ) *W*FkAT I O) **ExFCNrt I ) ) 

00154 

734 

c 


00154 

74* 

c 

SEARCH AND INCLUCC OZONE TRANSMISSION FOR THIS INTERVAL 



out 54 

75# 

C 


00155 

76# 


TRANS'TRANSfclRANSOil > 

W1 56 

77# 


CO TO 2150 

001 56 

76# 

c 


IKil 56 

79# 

c 

COf UTE TRANSMISSION CUE TO WATER VAFOi. 

XXi 156 

80# 

c 


OOl 57 

81# 

1200 

CONTINUE 

00160 

62# 


£ETA=1 .16# <il ,38* . 481 'Mif'RATlOr . 23 ) #=>*5 

0O16I 

65* 


A=ABSO;nn*W 

CJO 1 62 

84# 


B-eeTA<-irRATio*^4-.i5) >«a 

0O1 65 

65# 


C-ti 1 .4«4.9^€ETAWtf’-RATlC<=*t-l.i5) )*A) )<=Wx~.5) 

0J164 

6 G* 


TRANS^CXF<*-o£ =?C) ) 

mi 64 

87* 

c 


CO! 64 

88# 

c 

SEARCH AiC InCLLCE OZ'ZbE. TkANSMI 55 ION FOR THIS INTERVAL 

mi 64 

69# 

c 


00165 

90# 


lFal.0C.78) TRAN5 =TRANS*TkANS:«I ) 

00167 

91# 

2150 

RETURN 

0O17O 

92# 


FUNCTION TRANSOM l) 

0G1T5 

93# 


TRASUK^O. 

00174 

94# 


TRANSOM 1 . 

00175 

95* 


DO 10 J- 1,100 

razoo 

90* 


XKCEX^VCICHTUJ) 

00200 

97# 

c 


00200 

98# 

c 

use DATA WITHIN TVG CURRENT WAYENUCER INTERVAL 

00200 

99# 

c 


00201 

too* 


IFalfCEX.M:. I> CO TO lO 

00205 

101 # 


WE I GH^ve I GMTtt J ) - IISCEX 

00205 

102 # 

c 

GCfF-Lm: HDSNALD^O TRANSMI SS I CN FUNCTIONS FCEc OZOE USING TRACE 

00203 

105# 

c 

C0205 

ir** 

c 

CAS COEFFICIENTS 

00203 

105# 

c 


00204 

106* 


6 ^ AGSCRO* J ) XrOMlM J ) ) AT I O^C XFCN X J ) ) 

00204 

107* 

c 

FLAVIAN WEIGHTING FUNCTIONS TO THE TRANSMISSION VALUES 

00204 

108* 

c 

002*14 

109* 

c 


00205 

110 * 


TRASUH^TRASLM+EXFct-B) 4A-EI GH 

0 O 2 O 6 

ill# 

JO 

CCNT IMJE 

00210 

112 * 


I FaTRASUH . GT . 0 . ) TRANSOM TRASUH 

00212 

113* 


ROURN 

00215 

114# 


EM5 


EM2 CF OCH^ILATICN' 


NO DIAGNOSTICS 



rc*.s cccr / slant »cco f o;o /SLANT 
FCft 94L-06/20-16' 50 -0,0 ) 


BLOCK DATA 

5TCAAGE USCC CCCL^l) WZ-'.- CATAnO) 000002 ELAhA C-lNMlN-*? 1 CttJP.HO 
CCMMIN CLCCKS* 

C XW3 Z1 0014G4 


STORAGE ASSIGNICNT -BLOCK, TYPE, RELATIVE LOCATION, KAMI) 


CV.703 ft CY.VA7.V.1 AESCft 
0003 R 000632 CXFCNC 

ra<3 r 000454 resfon 


0003 ft 000620 ACS1KC 
0003 001154 EXFOO 

0003 00064* WEIGHT 


0003 CH.il 320 AESOftO 

0000, 1 OOOOOO l 


0003 0311*10 EXrOMD 

0003 R 000310 INTER 


OCUCil 

14 


BLOCK DATA 




00102 

24 


real inter 




00103 

34 


PARAMTTER 5lZTKA=lOO 



00104 

4* 


CCMOt/Z 1 / 

ACSCft*SlZTRA) , EXFCNtSlZTRA) , INTERNS I ZTftA) , 

00104 

54 


4 

R£5PCNnSlZ 

TRAJ , A£SORC<xiro , EXrCNCajO) , 

00104 

64: 


4 

V£lCHT<tSlZTftA) ,EXtvMD4SI2TkA) .EXFCNl^SlZlftA) , 

00104 

74 


4 

AESOROaSIZTPA) 


00104 

64 

C 





00104 

94 

C 

PARAurrcRS 

APPEARING IN ANALYTIC 

EXPRESSIONS CF TRANSMISSION 

00104 

104 

c 

THitOUGM WATER VATO; IN 

TIC WAVt NUT-CEK RANGE 25 TO 530 1/CM. 

00104 

114 

c 





00105 

124 


CATAaAESrfttf I > , EXFONa I ) , 

,1=1 ,211/ 


00105 

134 


1 670., 

o., 

1350., 

o., 

00105 

144 


2 2450., 

o.. 

2500. , 

o., 

00105 

154 


3 2550. , 

0., 

2100. , 

0-, 

00105 

164 


4 1250. f 

0., 

1050., 

0., 

00105 

174 


5 955., 

.45, 

710. , 

1.15, 

00105 

164 


6 410., 

1.65, 

290., 

1.95, 

00105 

194 


7 265., 

2,6, 

140. , 

3., 

00105 

bo* 


6 53.5, 

2.45, 

37.5, 

2.55, 

00105 

214 


9 30.5, 

2.75, 

20.5, 

1.55, 

00105 

224 


1 11. , 

1.65, 

6.05, 

2.3, 

00105 

23 4 


2 3.7, 

1.55/ 



00105 

24 4 

c 





00105 

254 

c 

PAKAHTl EftS 

APPEARING IN ANALYTIC 

EXPRESSIONS CF TRANSMISSION 

00105 

264 

c 

THROUGH CAkDZN CIOXICE 

an; water 

VAPOR IN THE wave NUMBER RANGE 

O.'i 05 

27* 

c 

350 TO 000 

1 /CM* 



fining 

5**4 

r 





OOl to 

29* 


CATA*ABSCft*n ,EXPCN*I) ,1=22,31)/ 


oona 

304 


l 2.6, 

1.35, 

2.1 , 

■ 6, 

D0110 

314 


2 1.55, 

.0, 

1.1, 

.0, 

OOllO 

324 


3 .62, 

.0, 

.615, 

.a, 

001 io 

53* 


4 .47, 

.0, 

.37, 

.0, 

OOllO 

344 


5 .29, 

.0, 

.23, 

.0/ 

00113 

35* 


DATA-ADSOkC- I ) ,EXFONC«I> ,1 = 1 ,10)/ 


0003 ft OOT144 EXrON 
O.XVJ I 00000 1 J 



92 


CV.U13 

36* 


1 .00145, 4.7, 

.0365, 

4.1, 


00113 

37* 


1 .1*0, 

3.1, 

l*T t 

2.2, 


00113 

30* 


t 6.95, 

0., 

4.*, 

l.i 


OOl 1 3 

33* 


1 .33, 

3.0, 

.115, 

3.6, 


00113 

40* 


1 .0105 

. 4.1 , 

.000 96 

4.7/ 


00113 

414 

C 






00113 

42* 

C 

FARAJ^ITERS Af PCARlNC |N ANALYTIC EXF-RESSICNS 

CF TRANSMISSION 

OOl 1 3 

43* 

c 

THROUGH 

WATER VAfvi 

\ IN THE WAVE NUTTER RANGE 

8' Hi TO 1200 1/CM. 

00113 

44* 

c 






DO IX 6 

4 5* 


DATA* ADS- :RaH lEWKNuf ) ,1=32,471/ 


00116 

46* 


1 . 1 70 , 

.775, 

.135, .020, 

.115, .860, 

.105, .660, 

OOl 1 6 

47* 


2 .095, 

.665, 

.091 ,.885, 

.091 , .005, 

.091 , .865, 

00116 

43* 


3 .091 

.665, 

.091 , .865, 

.091,. 885, 

.091 ,.660, 

001 i s 

49* 


4 .095, 

.660, 

.105, .030, 

.115, .795, 

,125,. 760/ 

00116 

50* 

c 






00116 

51* 

c 

FARA MITERS act-car 1 

NG IN ANALYTIC EXFRESSlCNS 

CT TRANSMISSION 

00116 

52* 

c 

THROUGH 

WATER VAFOR IN T>t WAVE NUCCR RANGE 

1 200 TO 2150 1/CM, 

OOl 16 

53* 

c 






O'.nzi 

54* 


CATAaABSCRttl) ,1=48 

,0*)/ 



PJ121 

55* 


l .2a, 

.42, 

.75, 

1.5, 


00121 

56* 


2 3.1, 

6.4, 

13.5, 

30.. 


00121 

57* 


3 45.. 

TT9., 

120. , 

220.. 


00121 

56* 


4470. , 

370., 

295. , 

67., 


00121 

59* 


5110. , 

235., 

570. , 

405., 


0012! 

60* 


6320., 

220 ., 

135., 

67., 


00121 

61* 


7 52., 

30., 

18., 

U.. 


00121 

62* 


6 7.7, 

6 . , 

4.2, 

2.6, 


an 2 i 

63* 


9 1.5, 

.9* 

.6, 

• 4 1 



00121 64* I .ZB, .2/ 

00121 65* C 

00121 L W*VtNUtCD< TNTLkVAL MLCHOJNIS 

00121 G7* C 

00123 63* DATAdlNTtR-xI) ,1 = 1 ,47)/ 


00123 

69* 

1 

37.5, 

62,5, 

07.5, 

00123 

70* 

1 112.5, 

137.5, 

1GZ.5, 

107.5, 

00123 

71* 

1 212.5, 

237.5, 

262. 5, 

287.5, 

00123 

72* 

1 312.5, 

337.5, 

362.5, 

387.5, 

1X1123 

73* 

1 412.5, 

437.5, 

462.5, 

407.5, 

OOl 23 

74* 

1 512.5, 

537.5, 

562,5, 

507.5, 

□0123 

75* 

1 612.5, 

637.5, 

662.5, 

607.5, 

00123 

76* 

1 712.5, 

737.5, 

762.5, 

707.5, 

□0123 

77* 

1 812.5, 

037.5, 

862.5, 

887.5 , 

03123 

78* 

1 912.5, 

937.5, 

962.5, 

987.5, 

00123 

79* 

1 1012.5, 

1037.5, 

1062.5. 

1087.5, 

00123 

00* 

1 1112.5, 

1137.5, 

1162.5, 

1187.5/ 

00125 

01* 

DA TA<i INTER c*I ) 

,1=48,85)/ 


00125 

62* 

1 1212.5, 

1237.5, 

1262.5, 

1287.5, 

00125 

83* 

1 1312.5, 

1337.5, 

1362.5, 

1307.5, 

not 24 

Avi* 

* 1412.5, 

1 437.5, 

148*, «. 

1 407,4. 

00125 

85* 

1 1512.5, 

1537.5, 

1562.5, 

1587.5, 

00125 

86* 

1 1612.5, 

1637.5, 

1662.5, 

1687.5, 

00125 

87* 

1 1712.5, 

1737.5, 

1762.5, 

1707.5, 

00125 

88* 

1 1812.5, 

1837.5, 

1062.5, 

1007.5, 

00125 

89* 

1 1912.5, 

1937.5, 

1962.5, 

1987.5, 

00125 

99* 

1 2912.5, 

2037.5, 

2962.5, 

2087.5, 

DOS 2 5 

91* 

1 2112.5, 

2137.5/ 



09127 

92* 

CATArtKESFCN* J) , J=1 ( 05>/P5*9./ 


00131 

93* 

END 





END CT CCMPlLATlCtf 


ND ClACN'.STtCS 
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FOR <$ TFM1AB/SLANT , TENT AC, TLMTAB/SLAnT 
FOR 50 a'J f O> 


suuRccrr i ne temtab 
tadtcm 


ENTRY POINT OOO 166 
ENTRY Pot NT 000177 


STCRACC USED' CCCE^ 1 ) WV.iZl'j CATAo.0 ) 0G0343 BLANK CCMM:m2) OOOOOiO 
CCHH1N BLOCKS 1 
CiG03 21 O' 'i 1 4 64 


EXTERNAL REFERENCES «£L OCK, NAVE) 

0004 EXP 

0005 >CKR3$ 


STORAGE ASStWCNT ofcLCCK. , TYPE, RELATIVE LOCATION, NAPE) 


OGGI DJ0002 10L 

0001 OGGI 37 «0t 

OOOD R 000145 DECREE 
OCTO R l VOGOl HEATER 

0C7.V.1 : 000314 j 

OOOO R 0003 lb WA VEND 


0001 

000071 

133G 

0001 

000123 

SOL 

0003 

001 01 0 

EXPO-ID 

oooa 

I 0-7131 1 

I 

0009 

R 000*454 

RES*"N 

av.ij 

000644 

WEIGHT 


OGGI tXVX.75 1370 

0003 COSaXt AESDR 

0003 LOO 144 EXPCN 

0X0 000325 IN Jr I 

0X0 R 000316 SFCC8 


0001 000091 20 L 

0009 000620 aes:rc 

0003 000632 EXFCNC 

0000 R 0X000* I MEN 
0000 R 000912 SLCTEA 


0001 


000053 

30L 

0003 


001320 

AE50D“O 

0003 


00 1 1 54 

E XryJO 

0003 

R 

000310 

INTER 

0X0 

R 

177013 

SL'IFEE 


ooioi 

14 


SLCROUT I KC TEMTAEkiCETECT.TEM^DK) 

OGlGl 

24 

C 


OOIOI 

34 

c 

THIS ENTRY POINT RETURNS A TEMPERATURE IN DECREES KELVIN WHICH 

octal 

44 

c 

MATCHES A CETECTCR INTENSITY. IF ND MATCH 1 5 MADE A ZERO IS 

ooioi 

54 

c 

RETURACC . 

DCUGl 

€>< 

c 


00103 

74 


REAL INTEN 

00104 

£4 


REAL INTER 

00105 

94 


PARAMETER S1ZTRA-1DG 

ocune 

104 


CCHMCN/Zl/ AES^^IZTRA), EXFCT4iSl ZTRA) , I NTER^SI ZTRA) , 

00106 

1 14 


4 RESfO/4*SlZTRA) ,AESZRCalEO , EXFCNCalO) , 

00106 

124 


4 WClGHTaSlZTRA) ,EXrOOiSlZTRA) .EXPCPOiSlZTRA) , 

00106 

134 


4 ABSOKxiSlZTRA) 

00107 

144 


DIhCNSiCN HEATER ft 100) ,CEGRC EalGO) 

DJI 10 

154 


BGOCYaV,T) =<*£. 934 9E _ 1 34V443) /aEXFcml .43654V) /T> -1 .0) 

00110 

1 64 

c 


0U11U 

1/4 

L 

r*jn; i n<_ 

UU11U 

1*4 

c 


DJI 1 1 

194 


1 = 1 

00112 

204 

10 

1 = 1*1 

00113 

214 


IF* I .GT.10G) GO TO 30 

00115 

224 


IF*CCTCCT.GT.HEATCR«I-1) .AND. DETECT. LE-HEATERrl)) GO TO 20 

00117 

23* 


GO TO 10 

00117 

244 

c 




00117 

25* 

c 

RETURN THE SOUGHT TEMPERATURE 


<xn t 7 

264 

c 



OC«t 

27* 

20 

CCNTINUE 


001 21 

26* 


SLa-CA=*CETCCT-HEATER*m/«l lEATER«l-ll -HEATER*! ) > 


00122 

29* 


SLCPEB=*CETCCT-MEATER*l-l> >/aHCATER*t> -HEATOVtl -1) > 


0D1 23 

30* 


TEMPO* - SL<PEA’-^CEGRCE*l -1 ) ♦SLOPEC^CCORCEaU 

. 

W124 

31* 


RETURN 


DO 124 

32* 

c 



DOl 24 

33* 

c 

ERROR RETURN ZERO TEMPERATURE 


GO 124 

34* 

c 



1X1 125 

35* 

30 

TEMPCjR=0. 


00126 

36* 


RETURN 


00127 

37* 


ENTRY TACTEKtTEM<R»TRANGE> 


00127 

36* 

c 



OG127 

39* 

c 

THIS EMI RY FOINT GENERATES A TEMPERA! URE/lNTENS I TT 

CAL I NATION 

COl 27 

40* 

c 

TABLE 1 -V.1 CCGRCES WIDE IN STEPS <Jr C7t DEGREE. 


00127 

41* 

c 



exit 31 

42* 


DCGREE*1)=TEMFOVTkANM1/2- 


□0132 

43* 


CXD 40 1=1 ,11X1 


00132 

44* 

c 


THIS INSTKUCNT“S 

00132 

43* 

c 

INTEN IS THE TOTAL EUACXDCCY INTENSI TY ACCEPTED BY 

an 32 

46* 

c 

RESPONSE FUNCTION. 


1X1132 

47* 

c 



00135 

46* 


INTEN=0. 


00136 

49* 


CO 50 .1=1,65 


00141 

50* 


iFuRESPCNrt J ) .LE.ti.) GO TO 50 


10143 

31* 


WAVEfO^ INTER* J) 


00143 

52* 

c 


AtC TEMPERATURE. 

03143 

53* 

c 

CorOTE A CLACK- DICY INTENSITY AT THIS WAVDAM3CR 

00143 

54* 

c 



00144 

33* 


SFCt^t^CCY^ AVUO.U,Gt*lE<*1/ > 


0014 5 

36* 


I NTEN- I NTEN* STCRC-RE ST<*U J J 


OJ146 

57* 

30 

CONTINUE 


00146 

56* 

c 



00146 

59* 

c 

MULTIPLY BY TVC 25 1/CM INTERVAL AND CIVICE BY 2. 


00146 

60* 

c 



00146 

61* 

c 

CONSTRUCT A TEMPERATURE CETECTCR TABLE. 


00146 

62* 

c 



00150 

63* 


►CATER * I > = l NTENS2 5 . 


00151 

, 64* 


XF«I.Ed.lOO> GO TO 40 


00151 

65* 

c 



00151 

66* 

c 

increment temperature. 


00151 

67* 

c 



00153 

68* 


CEGREEetl ♦ 1 > =CEGKEEa I > ♦' TR ANGE/ 1 00 . 


00154 

69* 

40 

CONTINUE 


00156 

70* 


RETURN 


00157 

71* 


END 



r*un 4 yr cC*^ =i Tl* t TO w ' 


C I AGNDSTICS 



9b 


Fc*.s VAP RC S / 5L AMI , VAf-nE 5 , VAPRES ASLANT 

for 9 < 4 .*oa/ 2 '.'.-i**y:> ar,,i'.j 


FUNCl ION VAFRES CNTRY POINT 000066 

M CRACE USEC' CCCE-l) W®» CATAal'i) WOlIJl BLAf* CWfflWI W.VJ.V.V.I 


CXTtANAL REFERENCES <*£LOCK, NA*C> 

O.V.-.3 ALOCU) 

0004 fcxr-c$ 

OLX'*5 F£KR3$ 


STCRACC ASStCNS-ENT oCLCCH, TYPE. RELATIVE LCCATlCK. N^O 


cxvjto 000023 injp$ 

COUP R 000013 FARTS 


oxo r 0 : 0:02 parti 

0000 R 017X04 PARTS 


0.700 R 000003 FARTZ 
ro.Vj R O./VOG FAR 77 


COLV.T R 000007 PART 3 OB R CCO005 PARK 

OtV.Xi R 0001/31 TRATIO COOG R GOOOOj VAr*ES 


ooiai 

1* 


FUNCTION VAFRES^T) 

uoicu 

2* 

C 


OQIOI 

3* 

c 

THIS FUNCTION COMPUTES A VAf-CR 

001 Ol 

4* 

c 

It U&ES the "GCfT'ORATCH* 1 RELATION OPE 

0OIO1 


c 

►CTECRO CGICAL TALt-tS*. 

OUlOl 


c 


DO103 

7* 


THAT ICe 373. 16/T 

am 04 

s* 


PARTI -d 7.9 '.i29^*(TRATI<>1.)) 

00105 

9* 


F AR TZ = 5 . 02 £06 *ALCC- 1 0*T RATIO) 

00106 

ip* 


PARTS=H 3.49149*mTRATIOI .) ) 

ocnc.7 

li* 


PART4-*** 1 ! • /flli'i . *<FART 6) > -1 . ) ^ . 13?tfE 

am in 

12* 


PART 7= It ,344*»ll . -el ./TRATIO) ) 

pm ii 

13* 


P ART3 -*£ mIO . «<= ART7 > - 1 . > *1 . 3«U GE- 07) 

ootiz 

14* 


PART 5 sPART2-FART3*FART4 -PARTI 

COM3 

15* 


VAFRES^IO . **PART5 

DOHA 

IS* 


RETURN 

00115 

17* 


ETC 


DO C F 

COMPILATION* NO CIA&N^TICS. 



MAP , S RACMAP /SL AN! , TL0007/ SLANT 
MAP Ci’.<23-i'.<V2'/«-nV 5'J -uO t > 


1. 

Lie 

Liv*«r<UB 

2 • 

IN 

MAIN 

3, 

IN 

racm:c 

4, 

IN 

cccf 

3. 

IN 

VAPRCS 

6, 

IN 

TRAMS 

7. 

IN 

TEMTAG 


ADDRESS LIMITS COlG OO 025412 040000 073577 

STARTING ACCESS 021320 

W^lCS DECIMAL 10507 ICAWC 141*60 CCA^ 


SEGfCNT MAIN 


OOlOQQ 02S412 040000 073377 


nswtc$/For 

N^iuc^/tcr 

M^CS/FC#t 

M^CFS/FCR 

WnXMl/FCR 

f^ocv$/roR 

KFTVt/TCR 

NCNVT3/7C« 

nclos*/fcr 

nlclks/for 

r©scL*/rc« 

NLf=CA$/FOR 

rorons/FCR 

BEO-OR/UTL 

BBIN/LHL 

B6INDT/LHL 

BOCCM/LTTL 

NININ1/F0R 

NINFTJ/FOR 

NDTJN3/F Oft 
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